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FOREWORD 



The Commission on Engineering Education, 
through its Advisory Committee on The Use of Com- 
puters and Mathematical Techniques in Engineering 
Design, sponsored a conference at the Univetsily 
of Illinois at Chicago Circle, Chicago, Illinois, April 
21-28, 1966, dealing with the impact of computers on 
education in engineering design. The confermice 
steering committee was composed of: 

Professor Steven J. Fenves, Conference Chairman 
Dept of Civil Engineering 
University of Illinois 

Dr. Frankun H. Branin, Jr. 

Data Systems Division 
IBM Corporation 

Dr. Suluvan a. Campbell 
Xerox Corporation 

Professor Brice Carnahan 
D^;)t of Chemical and Metallurgical Engineering 
University of Michigan 

Professor Samuel E. Shapiro 
Department of Systems Engineering 
University of Illinois at Congress Circle 



The conference was intended to serve the faculty 
members of engineering schools who are responsible 
for the devdopment of courses and curricula in en- 
glneeriBg design. Its immediate objective was to 
explore all the facets that interact <ni dasigni dedgn 
education, and computer technology. The long-range 
objective was to set directions for engineering 
education and research both in design and analysis, 
wherever the latter is expected to change under the 
impact of computer concepts. 

These proceedings comprise most of the papers 
presented by the participants, as vrell as the discus- 
sions whi^ followed in some cases. The conference 
program and a list of those who attended it are in- 
cluded. 

The confermce was partially Ananced by a grant 
from the National Science Foundation to the G<nn- 
mission on Engineering Education. 
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OPENING SESSION 

SEOTION A 

CHAIRMAN: S. J. Fbnves, Prof, of Civil Engineering 

University of Illinois 



RECEDING RAGE BLANK- NOT Pii Mcn 

THE IMPACT OF COMPUTERS ON EDUCATION 

IN ENGINEERING DESIGN 

J. G. R. LICKUDER 
Thomas J. Watson Reseiurch Center 
International Businees Machines Corporation 
Yorktown, N. Y. 

(Summary) 



The keynote is the rate of advance of the tech- 
nology that provides the medium and the means — 
and constitutes in large part the context and the 
object— of engineering design. It is a high note. 
Several of the basic capabilities of informational 
technology, such as capacity for storage and speed 
of processing, have been doubling every two years. 
Over all, though undoubtedly not in every sector 
of tile technology, a high rate of change seems likely 
to persist throughout the university years of engi- 
neers who are now freshmen and well on into their 
postgraduate careers. Engineering education must 
do its best to prepare them to expect, to master, to 
employ, and, indeed, to contribute to the development 
of the more effective ways of thinking and working 
that the advanced technology will make possible. 

At the present time, we see a major wave of tech- 
nical advance about to break upon the field of engi- 
neering design. It may be only the first of a succes- 
sion of waves, or it may be the great main force 
that ynll shape the coming generation. All that is 
sure is that the wave is massive and that its impact 
will be great. 

The essential nature of the cresting wave is 
caught up by the current phrase, **close interaction.** 

^e process of engineering design involves inter- 
play between creative and routine activities, be- 
tween activities that are heuristic, such as formula- 
tion of hypotheses, selection of approaches, assess- 
ment of values, and activities that are algorithmic, 
such as determination of costs and calculation of 
stresses. In design by classical methods, those two 
kinds of activity were insofar as possible separated 
from each other, the designer attending to the 
heuristic part himself and assigning to draftsmen, 
clerks, and computers the part that he could lay 
out step-by-step in explicitly detailed procedure. 
That approadi was successful in inverse proportion 
to the challenge of the problem, for only when the 
design was confined to preformed diannels of solu- 
tion could anything like a dean s^aration of the 
algorithmic from the heuristic be achieved. It be- 
came obvious, indeed, that in truly creative engi- 
neering design the two kinds of activity guide and 
support each other in intimate interplay and that, 
if there is such a thing as a phase of creative in- 



spiration followed by a phase of clerical execution, 
the total effort involves hundreds of thousands of 
such phases. 

Thus creative engineering design requires close 
interaction between the heuristic factor and the 
algorithmic factor. Because people are so greatly 
superior to machines in respect of the former and 
machines are so greatly superior to people in respect 
of the latter, the requirement is essentially a re- 
quirement for close interaction betwemi men and 
computers. And that is predsely what is afforded 
by the new informational technology, which makes 
it possible: 

1. To bring the time scale of man-computer inter- 
action down to seconds, or even milliseconds, 
and thus to substitute for the old, ccarse con- 
catenation of human operations and machine 
operations a fine, tight interweaving in which 
even the smallest area contains both the man*s 
heuristic and the machine*s algorithmic 
threads. 

2. To rescale the computer*s contribution from 
uninterrupted execution of an ad hoc mono- 
litiiic program prepared to solve an entire 
problem or carry out an overall task to respon- 
sive execution of general purpose atomic pro- 
grams prepared to place a broad repertoire of 
fundammital information operations at the de- 
signer’s fingertips. 

3. To free the designer from the dilemma of hav- 
ing either to plunge forward one creative step 
to another with only the support of his own 
engineering estimates and artistic intuitions 
or to wait after each creative step for slow 
and laborious implementation and testing. 

The close-interaction approach to engineering de- 
sign depends upon i^ergic use of new hardware, 
new software, and new techniques. In this con- 
ference, we shall examine all three. 

In the field of hardwEire, the main areas are 
storage, processing, transmission, and input-output. 
In storage, the most dramatic advance is the advance 
up the scale of capacity. Stores capable of holding 
a million analog images or a trillion binary digits 
and affording access in a few seconds to any image 
or any bit are in advanced stages of development. 



o 







The trillion-bit digital atores seem especially signif- 
icant to me because they wiU make it possible to 
maintain vast libraries of engineering information 
in readiness for rapid access and for processing 
under program control. 

In processing, the advance is multiple access 
through time-sharing and multiprocessing, with 
memory protection, dynaniic relocation, and related 
features that make it possible for a team or com- 
munity of engineers to use extmisive (and expen- 
sive) facilities in a coordinated and efficient way. 
In transmission, the main relevant movement is the 
growth of digital transmission services and the 
melding of digital transmission wi^ digital compu- 
tation — ^which will let tiie user work in or near 
his office with information stored at diverse remote 
locations and processed wherever the computers are. 
In input-output, the critical area of advance is that 
of consoles for on-line man-computer interaction. 
Most of the basic capabilities are provided by one or 
another of the new-generation consoles. As consoles 
move from their positions of relative obscurity in 
traditional computing systems to become the main 
foci of interaction between the users and the facili- 
ties of on-line information systems, there is inten- 
sive development aimed at improving the basic 
functions of display and control, integrating them 
more effectively for convenience and rapidily of 
interaction, and making them widely available and 
affordable. 

The software requirements imposed by engineer- 
ing design and related multiple^iccess, interactive 
applications are more basic and more difficult, in 
my assessment, than are the corresponding be- 
ware requirements. The general requirement is 
the ^^software base,** the on-line library of widely 
useful subprograms ; the on-line languages and pro- 
gram-processing programs through which users 
will call and apply the subprogram from the library; 
the on-line programming systems that will facilitate 
preparation, testing, modification, documentation, 
and communication of new programs and subpro- 
grams, and on-line, controlled-access data banks. 
Problems raised by those general software require- 
ments are being attacked vigorously now, and the 
most basic requirements are gradually being met 
by research and development projects in universi- 
ties, the automobile and aerospace industries, the 



computer and communication industries, and several 
governmental and government-related non-profit 
organizations. Only two years ago, a large fraction 
of the progress in this area could be subsumed by 
citing **Sketchpad** and other programs at Masna- 
chusetts Institute of Technology in Cambridge and 
at the Lincoln Laboratory, the DAC-I Project at 
Gmeral Motors, and related work at Bolt Beranek 
and Newman, System Development Corporation, 
and very few other othw places. In the interim, the 
field has burgeoned. 

Undbr **new methods and techniques** in the 
phrase **new hardware, new software^ and new 
methods and techniques,** I meant to include new 
ways of using on-line facilities for engineering de- 
sign. It is now well understood that the main pro- 
mise of the on-line approach is not merely to 
automate drafting, but no one I know thinkn he 
yet sees clearly just how designers will use the new 
facilities ^r just what the process of design will 
be within the new context. The new ways of work- 
ing and thinking will evolve out of experimental and 
pilot-operational use. One of the main responsi- 
bilities of educators, therefore, is to incorporate such 
use of on-line facilities into the process of education 
in engineering design. Students should learn the 
new methods and techniques not so much by being 
taught them, as by contributing to their conception 
and development. 

What beyond on-line interactive design? Perhaps 
it is not too soon to think about that. My own best 
guess is that, when many engineers are working 
within the context of an on-line community, time 
will be a marked acceleration in the understanding 
of the heuristic aspects of design. More and more 
of the exploration of partly random, partly controlled 
variations on themes, for example, will be carried 
out by computer programs, and designers will tend 
increasingly to become the strategists, rather than 
the conductors, of explorations and experiments. 
Designers will function also as evaluators, of course, 
but perhaps even in that area they will become 
strategists — students of the processes through which 
**time** performs its ultimate evaluations. In any 
event, new possibilities and new opportunities will 
arise, and (from the somewhat non-central but none- 
theless interesting point of view from which I see 
it) the future of engineering design looks both chal- 
lenging and exciting. 
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SOME COMMENTS ON THE FUTURE PRACTICE 

OF ENGINEERING 

C. L. Miller 

Massachusetts Institute of Technology 
Cambridge, Mass. 



THE YOUNG ENGINEER— 1976 

The young engineers who will be 28 years old 
in 1976 will be entering college as freshmen this 
fall. While 1976 sounds like a long way off to you 
and me, age 28 sounds quite young to all of us. Cer- 
tainly we will want our 28-year-old engineers to be 
well prepared for their professional work. This 
means that the formal education we give these 
entering freshmen during the next four to eight 
Shears should reflect engineering as it will be prac- 
ticed in 1976, that far-off date. 

THE NEW TECHNOLOGY— 1976 

The new technology for engineering which is 
already available has been excellently described in 
the previous address by Dr. Licklider. With our 
concern for the 28-year-old engineer in 1976, we 
must project this new technology ahead for at least 
a decade of what is sure to be a period of accelerated 
progress. 

While no one can predict what the technology for 
engineering will be in 1976, it will surely be much 
more advanced than that available for exploitation 
today. One conservative approach to projecting 
ahead to 1976 is to say that the most advanced of 
the new technology which is just emerging in re- 
search environments in 1966 will be commonplace 
in design and production environments in 1976. 

THE PAST DECADE— ELECTRONIC 
CALCULATION 

The past decade might be characterized as the 
decade of electronic calculation in engineering. 
Since 1966, excellent progress has been made in the 
application of the computer for engineering calcu- 
lations and in the use of the computer as an aid in 
engineering analysis. In engineering schools the 
focus has been on the teaching, of elementary pro- 
gramming and an introduction to appropriate nu- 
merical methods. In addition, extensive use has 
been made of the computer in engineering schools 
for research calculations. During the past decade, 
interest in the potential role of the computer in 
engineering design has developed, and some progress 
has been made. Excdlent examples of such prog- 
ress will be presented during the balance of this 
program. 



THE NEXT DECADE— INFORMATION 
SYSTEMS 

As we look ahead to the next decade^ there is 
ample evidence that it will be characterized as the 
decade of engineering information systems. What 
has been done in the past for engineering calcula- 
tions will be done in the future for the whole area 
of engineering communications, decision-making, 
and systems design. We will be concerned with the 
computer as a component in large-ficale engineering 
information systems which will encompass the total 
spectrum of information acquisition, processing, 
storage, retrieval, and communication. We will 
be concerned not only with computer technology, 
but also with information technology in the broad 
sense of the word. 

THE NEW CAPABILITY— A NEW 
PRACTICE 

Information technology and the information sys- 
tems of 1976 will represent new kinds of capabilities 
for engineering practice, not just improvements on 
the old capabilities of past decades. To exploit the 
new capabilities will not be as simple as adding 
them to the old capabilities. It may take a rather 
complete restructuring of the practice of engineer- 
ing and wholly new approaches to the execution 
of the engineering process. 

It is reasonable to anticipate that the new tech- 
nology providing new capability will lead to a new 
practice of engineering. By a new practice is meant 
signiflcant changes in how engineering is ac- 
complished. The relative roles of men, machines, 
and organizations may be considerably different 
from that which we observe in engineering practice 
of today. Who does what, when, where^ and how 
will change. Engineering will still be an art in 
1976, but it will be quite a different art. 

THE NEW PRACTICE— A NEW MAN 

If we are to have a new practice of engineering 
in 1976 based on new capabilities, it follows that 
we may need a new kind of man to exploit these 
new capabilities and practice the new art By a 
new man is meant one who is educated and trained 
for the new art and who has the talents and compe- 
tence to use the new technology. 
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Engineering schools look ahead with respect to 
engineering science-oriented teaching and building 
on the science-oriented research. However, there 
is a tendency to look backward with respect to en- 
gineering practice-oriented teaching and to do little 
or no practice-oriented research. The art of engi- 
neering tends to be taught on a *i;his is the way it is 
done, boys** description of current and past practice. 
If the practice of engineering is going to change in 
the next decade, we have a responsibility to look 
ahead and prepare these young men for the future, 
rather than for the soon-to-be-obsolete methods of 
1966. 



KEY NOTES— KEY POINTS 
A kesmote address riiould presumably make a few 
key notes, or key points, to stimulate discussion dur- 
ing a conterence. It is not alwa3ni clear just what 
key points a keynote speaker is trying to make and, 
in order to bring mine into sharp focus, I will state 
them in summary fonn. My points are a combina- 
tion of speculation, opinion, and proposal. I trust 
that you will not entirriy agree with all those I am 
going to make because, if you do so, the keynote 
address will not have smwed its basic purpose— to 
stimulate discussion. To avoid the remote possibility 
that you might agree, I will deliberately state a 
somewhat extreme position which is sure to be some- 
what controversial. You should bear in mind tiiat I 
am speculating on engineering design practice as it 
might be in 1976 and not as it is in 1966. 



The ten key pohits I would like to make are as 
follows: 



1. Engineering design will be accomplished by 
an Site of gifted, brilliant; creative, highly 
Sailed, individual professional designers — 
hereafter called **mastor designers.** 

2. The master designers will be concerned with 
generating design concepts and with deciding 
between significant alternatives for broadly 
defined design problems. 

3. The master designers irill work directly and 
interactively with highly responsive design 
systems based on powerful and large-scale 
information systems, one component of which 
will be ultra high-speed computers, but the 
heart of which will be informatioh files and 
totally now kinds of software for working 
with these files. 

4. The engineering design process will be formal- 
ized and organized into a highly sophisticated 
strategy for the optimum use of wiaw and 
machine. The total design ss^stem will in- 
volve a highly sophisticated machine, a highly 
sophisticated man, and a highly sophisticated 
theory of design and decisiop-making. 






6. All aspects of engineering design which can 
be programmed will be programmed, includ- 
ing aspects of design which at this point in 
time are not wdl structured or wdl under- 
stood, resulting in a significantly different de- 
sign process from what we know today. 

6. Communications, in the broadest sense of the 
word, will be at least as important as proces- 
sing in the design system of 1976, and both 
will relegate calculation to a minor subset 
of total design activity. The conventional 
engineering drawing per se will cease to be 
the prime commiinicationB medium for engi- 
neering design decisions. Dynamic communi- 
cations will replace static communications. 

7. Design activity in an organization will be 
integrated with other decision-making and 
action-initiating activities, and \rill be con- 
ducted interactively with the total environ- 
ment and operations of the organization. 

8. The design systems for the late 70*s have 
not yet been fully designed, nor have the de- 
signers to use them been designed. The initia- 
tive and leadership for accomplishing such 
designs should be taken by the engineering 
schools. 

9. Engineering design education during the next 
five years should be highly research-oriented. 
The best design education we can give our 
students during the immediate years ahead 
is to make them junior partners with the 
facully in conducting design research and 
development. 

10. Future formal education for engineering de- 
sign practice will require a study program 
through the Fh. D., plus additional require- 
ments. The program should be highly sriec- 
tive and only those young men who are truly 
gifted for engineering design should survive. 
The design-oriented Ph. D. program for the 
future professional designer should emerge 
as the highest attainment in engineering edu- 
cation, involving all of the conventional re- 
quirements for demonstrated research capa- 
bility, plus demonstrated original design 
capability. 

THE PROTOTYPES 

While the ten key points represent an extreme 
view, they stem from extrapolation to 1976 of ex- 
perimental systems already under development in 
1966. In particular, they are heavily infiuenced by 
current work on the devriopment of the Integrated 
Civil Engineering System (ICES) underway at the 
MIT Civil Engineering Systems Laboratory and 
the large-seal^ time-sharing hardware/software 
devdopments underway in industry and the uni- 
versities, such as in MIT Project MAC. The 1976 
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yenions of these prototype systems, even if pro^ 
jected in a very conservative way, provide the basis 
for the ten key points. 

The protot^>es for the hardware/software sys- 
tems are under devdopment. The evidence that 
the technology will be available is quite clear. The 
evidence that a man will be available to use this 
technology is not as clear. That the engineering 
organizations will be available is even less clear. 
There is very little evidence that organizations are 
responding to the implications of such change by 
advance planning and preparation. However, it is 
anticipated that the prototype organizations will 
emerge in the next five years, and competition will 
take care of the balance in the follovdng five years. 

RATE OF CHANGE— RESISTANCE 

Change of any kind is always resisted by a sub- 
set of the people and organizations which are affect- 
ed. The subset includes those who are are already 
content but who fear that their relative iK>sition 
will be adversely affected. Engineering is no ex- 
ception. Many will adopt the attitude that the 
new capabilities are not needed. Historically, this 
has always been the case. (It wasn’t too long ago 
that the log-table people argued that desk calculators 
were not needed, and it was just a few years ago 
that many questioned the need for computers in 
engineering calculations — some still do!) 

There is no question but that the impending 
change in engineering practice will be resisted in a 
variety of subtle and not so subtle wasns. There is 
also no question but that there will be pressures 
for change. The only question is how the resistance 
will affect the rate of change. The rate of change 
is a major unknown at this time. It is quite possible 
that we may indeed fall way short of the situation 
described in the ten points by 1976 if the anticipated 
resistance to change is sufficiently strong. 

THE CHALLENGE 

Our engineering schools can have considerable 
infiuence on the rate of change in engineering prac- 
tice during the next decade. We can be an effective 
source of resistance and slow down the rate of 
change. Or we can accept the challenge of being 
an effective source of leadership. The choice is 
ours. This Conference on the Impact of Computers 
on Education in Engineering Design is particularly 
timdy as there will be an impact in the immediate 
smars ahead. We have the opportunitir now to ac- 
cept the challenge of being the designers of that 
impact and the instigators of change. The alterna- 
tive is to be followers and to face a future of 
responding to an impact which results from external 
diange. If we do not accept the challenge by con- 
scious decision, we decide on the alternative by de- 
fault. 



Perhaps the most intriguing aspect of the chal- 
lenge we face is that there are a variety of ways to 
pursue the challenge. Like all real design problems, 
there is no one solution. There are many solutions. 
Every school and every individual can participate in 
the search for good designs. We are sure to hear 
and see evidence of this during the conference. All 
of us have a role to play in designing the new prac- 
tice of engineering for the 28-year-old engineers 
of 1976. 



PAYNTER: Taking MUler^s fourth and tenth points 
I get the distinct impression that the person we 
are reaUy talking about bears a strong resem- 
blance to the person we have heard about in the 
past who was eaUed the creative mathematician, 
I believe that creativity in basic science and 
in mathematics is not as different even now as 
you seem to make it out to be, or hint that it may 
be, from creativity in engineering design. And I 
definitely feel that, in the period ten and twenty 
years from now, the man whom you will screen 
out at this level of the PhJD, and creative engi- 
neering design d la computer wUl be a man who 
simply has a tremendous, rational, creative gift. 
If you are talking about the strategy and the 
development of strategy, certainly no discipline 
has done more incisive, introspective looking at 
the strategy for creation of formal systems than 
mathematics itself, I would like to hear a 
response to this, 

MILLER: 7 don*t know whether I can give you a 
very eloquent response, I certainly toUl be wWing 
to admit rather candidly that, as you know, we 
are trying to practice a little of what Pve just 
been preaching in our own effort of taking a 
responsibly role and attempting to design and 
develop these new design systems. Some of the 
most able and valuable young men working on 
that project actually did their undergraduate work 
in mathematics. The young men who are making 
the most notable creative contribution to the de- 
sign of the ICES system are people who have a 
great flair for mathematics, a/nd very abstract 
mathematics as far as that goes, Vm unsure of 
how far what you say is correct in terms of deci- 
sion making. Even in 1976, udlh a very powerful 
design system of the type we are discussing; we 
are going to be working with incomplete and 
inaccurate information with problems for which 
there are resource eonstraimls, that is, limits on 
the available amounts of money and time. In this 
whole area of compromise and trade-offs. Pm sUU 
betting rather heavily on the man to play a major 
role in making many of the judgment deciMms, 
but he win be assisted in these by some very 
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aophistieoited, very fomuU strategies behind the 
scene represented in terms of the programming 
therein, which witt fadUtate his ai){Uty to make 
what may in the final analysis he judgment deeu- 
stone, Whafs going on inside wiU, in fact, he a 
representation of a degree of mathematical sophis- 
tication that one would associate with the best of 
creative imdhematies. What goes on outside when 
he turns around and walks a/way from that con- 
sole, / think, wHl still represent charaeteristies 
that we associate with the engineer today, at least 
our image of a good engineer, that is, his ability to 
make decisions in the fsuse of uncertainty and 
constraints and to do very well in the major 
decisions he makes, 

LICKLIDER: I would like to add a couple of per- 
haps divergent comments here. First, I largely 
agree with Paynter that creative mathematics and 
creative design, within the limits of what / know 
about it, have a lot in common. Second, / want 
to ^suggest that if there is any thought in your 
mind that designers will be raiher a drug on the 
market because only a few will be able to handle 
all the jobs, I would disagree strenuously be- 
cause it seems to me that the search for solutions 
in problem space, or the search for good and 
reliable synthesis, is an extremely difficult and 
low probability pursuit. As soon as we have the 
apmratus, there are going to be facilities for 
^ng a lot of that kind of thing. Then, I think, 
it win become important to do very much more 
than we now do, to carry planning which now, 
even at the international level, is about 1 1/2-or 
2-step affair, WeU, to carry design very much 
faHher, into working designs and into planning, 
will take all the inteUectual resources there are. 
One final thought is that if we ever do start to 
saturate what one might call modal design, that 
is, the finding of straightforward solutions to 
design problems, which I don*t think we have 
come nmr^ yet, then there will be a great premium 
on aohieving the odd-baU, but happily very good, 
designs tlud you ean*t get to along the standard 
channels, and that will readily take a lot of intel- 
lectual effort, 

MILLER: Let me follow up on LieMider^s comment 
by noting that at no place did / say we needed 
just a few of these master designers. By analogy, 
in the last ten years, with electronic computers 
we didn*t just chew up in a relatively few miUi- 
seeonds aU the ealcuUdions that we would nor- 
mally have done over many days. We simply 
went al^ad and did a lot more compuHng than 
otherwise would have been possible, / think this 
is also true of design. Certainly we wid do an 



enormous amount more of design, of a different 
kind and a higher level than we did before, I do 
feel that the effective design technician, while he 
will still exist and play a roUf-and there wiU be 
relatively large numbers of themr-wHL operate 
differently in terms of where he fits into the pic- 
ture, Bui fAo basio point is that the whole new 
technology is meaningless, or at least it does not 
achieve its full potential, unless we have the 
equivalent of the man I was labelling master de- 
signer to really master the technology that toiU 
be available, 

RESWICK: It seems to me aU the trouble started 
back with the Tower of Babel when people be- 
came unable to talk to each other, Licklider 
t^ked about the problem of design by teams, that 
is, by bringing various kinds of people together 
to achieve solutions. He also has had experience 
with group dynamics, I am thinking of what 
happens when an architect and an engineer get 
together to work cooperatively, to do creative 
design or innovation. It seems to me that lAck- 
li^ was saying that the computer offers a solu- 
tion to the Tower of Babel, that everyone wiU 
have to talk computer language and therefore will 
have to talk in the same language, and therefore 
wiU be able to understand each other, and that out 
of this can come optimum creative solutions. And 
. yet, what actually does happen when such a group 
of people meets to do a creative job; and how does 
the best kind of idea come out of this group? 
So many things happen which are way beyond 
the information or the knowledge content of the 
PMple, What is the information content of crea- 
tivity in a group, and when can a computer system 
ever reach the capability of transferring all of 
the things that must happen between people in 
sueh a situation to reach something like an opti- 
mum decision? 

LICKLIDER; I certainly have no adequate answer 
or solution, but I would suggest a couple of things. 
First, I wouldn*t hold that all the people would 
have to speak the computer's language, I think 
the trend is for the computer to speak the lan- 
guages of the various disciplines. The second 
point would be that much of the miounderstandf- 
ing between people from different disciplines is 
misunderstanding based on the fact that much of 
theyr talk is merely naming, and names are static. 
With the aid of the computer, one could substitute 
for the mere passage back and forth of static 
messages, the demonstration of behaviors. If you 
see the behavior, even if iffs just pictures of be- 
havior on a scope face, Ws better than passing the 
words book and forth. 
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1 . INTRODUCTION 
Over the past decade, the computer has been play- 
ing an ever-increasing part in nearly all of our 
activities at the Bell Telephone Laboratories. Since 
these activities span a wide range of technical areas 
— ^from basic ph3rsical and mathematical research 
to the design of pushbuttons for telephone sets-— 
it would be impossible to make a complete accounting 
of all the ways in which we have found the computer 
to be a useful tool. We will take two examples of 
practical design problems in the co mm unications 
industiy and give a brief history of how computer- 
aided techniques have developed over the years in 
each of these. By so doing we will be able to give a 
flavor of the manner in which the computer has 
changed, and is continuing to change, the way we 
'’.o our work. The two examples themselves serve 
to illustrate the variety of our problems. The first 
— ^that of a communication transformer design — 
illustrates some of the problems we face in selection 
and design of components, while the second — ^that 
of designing large systems — shows what we need 
to do in order to assemble hundreds of thousands 
of components into working systems. 

2 . COMMUNICATION TRANSFORMER DESIGN 

Quite naturally communications ^uipment is full 

of all sorts of broadband transformers, each of 
which must be designed to meet the particular cir- 
cuit requirements dictated by its environment. The 
history of the devdopment of computer techniques 
in this area is particularly enlightening for it shows 
how the computer has been used at many levels, 
ranging from specific one-shot calculations to a 
highly complex set of programs which virtually 
automates the entire process from design concept 
to the generation of detailed manufacturing infor- 
mation. 

The word design as we here apply it to trans- 
formers is really a misnomer in that it encompasses 
a good deal more than physical realization itself. 
Also involved in the overall problem are the pro- 
cesses of : 

Search : Can an existing transformer handle this 
job? 

Synthesis: From the given electrical circuit re- 
quirements, what are the parameters of the 
equivalent circuit of the transformer? 

Design: What are the phj^ical quantities re- 



quired to realize these equivalent circuit pa- 
rameters? 

Characterization: What is the total performance 
of the transformer we have either selected or 
ssmthesized? 

Specification: What is the detailed information 
required for manufacture of this transformer? 

Before the availability of the computer, these 
basic functions were carried out manually by refer- 
ence to catalogs, rough slide-rule calculations, labora- 
tory tests and drafting. We will now attempt to 
show how, over the years, increased use of the 
computer has gradually automated each of these 
functions. 

Analysis 

The first use of a computer in the design process 
is tsrically that of analysis: How will a given 
device or ssrstem behave under various conditions 
of operation? The first step in performing such 
analysis is to develop a mathematical model which 
suitably represents the ph3rsical situation. This 
requires a combination of intensive theoretical and 
experimental studies. In some cases — as in linear 
lumped-parameter filters — ^the approximation is very 
close, and one is able to make calculations with 
confidence that they will be accurate within a small 
fraction of a percent In others — as in modeling the 
behavior of a lypical modem central ofiice — ^we are 
unable to develop a model simple enough to handle 
on a computer which is, in turn, sufiiciently realistic 
to provide more than fairly gross approximations to 
the s3Tstem*s real behavior. In our example of 
designing broadband transformers, the model is in 
the form of an equivalent circuit made up of linear 
lumped-parameter dements. 

Next one must develop a method of relating the 
values of the parameters of the mathematical modd 
to the physical characteristics of the device or sys- 
tem being studied. Similar relationships between 
the behavior of the model and that of the phjrsical 
s3Tstem are needed in order to be able to interpret 
the results of the study. In our example the phy- 
sical characteristics of the transformer are core 
size and permeability, wire size, turns ratio, insu- 
lation thickness and several others. The parameters 
of the mathematical modd are the values of the 
elements in the equivalent circuit, and in this 
case the behavior of both the model and th e device 
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is either the transient response or the frequency 
characteristic curves, depending upon the trans- 
former’s application. Since our emphasis here is 
more on the design process itself rather than on 
broadband transformers, we will not go into detail 
as to the specific equivalent circuit which was used.^ 

In ihe first step of introducing the computer to 
the problem of design, we program it to accept phy- 
sical parameters of a given transformer, calculate 
the values of the dements in the equivalent circuit, 
and then i)erform standard network calculations on 
that drcuit. This allows an engineer to select a 
set of physical parameters, examine the resulting 
performance^ and zero in on a satisfactory design by 
repeated runs, varsring the phsrsical parameters on a 
cut-and-try basis. 

Synthesis 

One would initially suppose that the problem of 
sjmthesis would simply be a matter of reversing the 
problem of anabasis, that is, solving the equations 
for the physical parameters in terms of the desired 
characteristics. Howevmr, this is generally not the 
case; there are only a small number of phsrsical 
parameters which can be altered and one usually 
has a larger number of external characteristic val- 
ues to be satisfied. Thus a good deal more investiga- 
tion, both experimental and theoretical, is needed 
in order to determine what specific external meas- 
urements are of greatest significance, so that the 
inverse solution can be made possible at all. We 
are generally faced with a set of non-linear equations 
which defies solution even with a computer at hand. 
If we are fortunate, some further analysis will pro- 
duce a few gross relationships whidi allow at 
least an approximate solution to these equations. 
Now we are faced witii approximate answers to a 
mathematical model which in itself is an approxi- 
mation of the actual device. The next step is utili- 
zation of the actual device. The next step in utiliza- 
tion of the computer is to write a program which ac- 
cepts desired external diara^ristics as input and 
uses the rough approximations to estimate the ' 
physical parameters, providing input directly to 
the analysis program. 

This then allows tiie engineer to zero in on his 
final design with fewer computer runs. It also 
shows us the start of a system of progr am^ rather 
than of a single one-shot effort. The next three 
stages of development in our transformer design 
example are somewhat indep^dent— one aids the 
designer in his catalog search (after all we may 

^For a treatment of the broadband transformer design 
problem in greater detail, see "Ck>mmnnieation Transformer 
Engineering by Computer,** by C. M. BaUey, P. J. Kasper 
and D. Eats, in the Proceedings of the 1966 Electronics Com- 
ponents Conference. 



already have a transformer in stock which fills the 
bill), the second generates the detailed manufactur- 
ing information required by the factory, and the 
third allows his cut-and-try process to take place 
more conve^iiently and far more rapidly. 

Search 

Once a suitable model of our transformer exists, 
then a specific transformer can be exactly described 
by a set of numbers. If we have recorded all pre- 
viously designed transformers in computer memory 
(typically on a magnetic tape “catalog”), we are 
able to match the parameters of the desir^ trans- 
former against those for which designs already 
exist. Should any fall within suitable tolerance, 
then the analysis program is brought into play 
and the designer can see how well any existing 
design fits his requirements. If there are any that 
come close to fitting the requirements, they are 
so indicated and their complete characteristic curves 
are presented to allow ihe engineer to see if one 
of them fits the bill, (itherwise, we go through 
the synthesis phase (and, of course, add this new 
transformer and its characteristics to the catalog) . 

Output 

Once a design has been completed, the engineer 
has (as a direct result of the process described 
above) a complete functional characterization of 
the transformer. But the computer now “knows” 
^1 of the physical quantities required for actual 
manufacture. Therefore, a part of our growing 
^ign system is the production of manufacturing 
information in a format which is directly acceptable 
to the factory. 

Let us see what the computer has provided us 
with so far. Our program system provides the 
engineer with a significant saving in time, cost 
and accuracy. The S3rstem, through its search 
routine, will minimize the number of new designs. 
If a new design is necessary, the frequency band 
of interest is centered, thereby providing greater 
possibility for multiple-use transformers. This cen- 
tering procedure will tend to minimize overlap in 
designs; hence, a tendency toward fewer different 
designs. By iUling out the transformer design re- 
quest, the engineer can have witiiin a day’s time 
complete information on any transformer which 
meets his needs and/or design information, if it is 
possible to produce such a transformer. 

Although we do indeed have a system of pro- 
grams, it was designed in such a fashion tiiat part 
of it may be used independently. For example 
both the si^ithesis-design and the analysis portions 
have found extensive use as individual programs. 

Direct Interaction 

About the time that our transformer sdect«nd- 
design systmn became operational, a completely new 
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facility was made available to us at the Laboratories 
on an experimental basis. The present talk is not 
intended as a description of that facilily, for other 
talks in this session will go more deeply into 
similar developments. Briefly stated, this equip- 
ment, known as 6RAPHIC-I, provides a graphical 
input/output capability whereby a user can com- 
municate with a large computer via a CRT display 
and a light pen. Physically it consists of a CRT, 
a k^board, a light pen, and a smaller computer 
(a PDP6) with an added memory, plus hardware 
which allows data transmission to and from the 
7094. It differs from other similar devices (sudi 
as SKETCHPAD and the General Motors DAC 
Ssrstem) primarily in that the local computer and 
its added memory provide a good deal of stand- 
alone capability. Thus communication to and from 
the 7094 tends to take place relatively infrequently 
but with fairly large chunks of data handled in 
one interaction. GRAPHIC-I is really an experi- 
mental vehicle leading to the design of GRAPHIG- 
II, which win be used to provide graphical com- 
munication with our GE 645 S 3 rstem about a 3 rear 
from now. Its major limitation is one based on 
7094 software — ^we didn't feel it worth while on 
our particular time scale to provide capability for 
interruption of the 7094 on other than a between- 
job basis. For this reason, one might have to wait 
up to a few minutes to capture the 7094, but from 
then until the transaction is completed, he com- 
mands the 7094*s full attention. 

So much for the equipment. How has this 
changed our ouijook on design in general, mid on 
broadband transformer design in particular? Let's 
review the basic method of design we have used 
all along the way. Essentially it has been to guess 
what the proper ph 3 nsical configuration would be, 
anal 3 Tze this and then make another guess. The 
computer's role has been partiy clerical and partly 
calculational. It has been clerical in that the pro- 
gram first looks in a catalog to see if a new design 
is really needed. It has also been derical in that 
it "r^embers" past guesses and their consequences 
and in that it presents the finally agreed-on design 
in the proper form for manufacture. The calcula- 
tional roles of the computer have been to perform 
the analysis and to assist the engineer in making 
better-educated guesses. Now what the GRAPHIC 
facility provides is greater simplicity of use (the 
engineer no longer has to punch cards and assemble 
mysterious decks in the proper order)* and an 
order of magnitude in compression of the time 
scale. A complex design might require a dozen 



* At present he has to learn other mysterlons things lUf* 
poshing «aight buttons,*’ waving the light pens and the like. 
Onr ultimate aim is to make the system as simple as possible 
for the user. 



or so attempts — each depending heavily on the 
results of former trials. Using typical batch proc- 
essing procedures, this might stretch the overall 
time between orisdnation of the needs and realiza- 
tion of a satisfactory design to several days. Of 
course, during this time pmriod the engineer is 
presumably off doing something dse, but the turn- 
around time is still measured in dasrs. Having a 
^ectly interacting facility for entering data, view- 
ing results and entering new data can reduce this 
to a matter of minutes. 

While we are on the subject of eleclaronic com- 
ponents, we might profit by looking into another 
application of graphic consoles. A large part of 
what we generally refer to as design consists of 
looking through catalog information for a com- 
ponent which meets various combinations of re- 
quirements, including electrical and physical 
characteristics, costs, reliabilily, and numerous 
others. Although we have not 3 ret built up a file, 
we have developed a structure for such an applica- 
tion and have pla 3 red with it enough to be convinced 
that this will be of great value. What we are 
really doing is searching a huge file structured in 
the form of a tree, and at each node asking the 
user which branch he wishes to take. Although 
the identical searching process is feasible using a 
t3T>ewriter as the input/output device, the CRT 
light pen is far faster in displaying the information 
pertinent to each node, and it is far easier for the 
user to point to something by a light pen than to 
tsrpe his request out. 

3. SYSTEM DESIGN 

In addition to designing components we are, of 
course, interested in gatiiering these components 
into various levels of sub-assemblies which are 
built up to form large systems. It is not at all un- 
usual today to discuss a system with nearly a million 
dectrical components and hundreds of thousands of 
interconnecting wires. In this area, I am sorry to 
say, we can't point to as neat and orderly a growth 
of automation as in my first example, although there 
have been a few notable approadies to complete 
design automation systems. 

When we look at system design, the first problem 
that faces us is the sheer bulk of the data, together 
with the complexity of the interrelationships among 
subsets of the entire data. For this reason, most 
design automation efforts aimed at large systems 
have first of all attadced the clerical, record-keeping 
problems rather than the activities involving some 
creativity. Automation of record-keeping is vir- 
tually mandatory here, not only from the point of 
view of economy, but also in order to provide the 
host of design engineers up-to-date and consistent 
information as to the current status of the system. 
Furthermore, one can ask questions about the ap- 
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proximate cost of the entire system or of a sub- 
system, or find out how many times and where a 
particular component is used, and receive answers 
economically and quickly. With manual methods 
of keeping records, either of these questions would 
require so much thumbing through stock lists that 
the answers would cost too much, and by the time 
they were generated they might well be out of date. 

It is appropriate here to point out one of the 
major problems we have encountered in using the 
computer to document large ssrstems. Especially in 
engineering development efforts, the nature of the 
information we wish to keep is continually changing. 
Furthermore, the various people involved with the 
job are interested in different subsets of the entire 
information. However, the data-processing soft- 
ware we have available to us today renders changing 
the format of a file very difficult, since it necessitates 
changing each individual program which works with 
the file. Thus we are faced with a combination 
of two problems: a good deal of our programming 
effort is wasted by trivial but time-consuming edit* 
ing to introduce changes in file or input/output 
formats; and many potentially useful ideas are 
abandoned because of tlie cost involved in imple- 
menting them. 

A second area in which automated design of 
large eorstems has suffered stems from a heavy de- 
pendence of most of them on the detailed nature 
of the equipment being designed. That is to say, 
things like the type of logic used and the size and 
terminal identification of sockets tend to become 
so deeply ingrained in these programs that we al- 
most have to write a complete new ssrstem of pro- 
grams to handle a new line of equipment. 

So much for the negative side. There are a 
dozen or so companies now designing and manufac- 
turing large electoonic ssnsitems, and each has its own 
computer-based record-keeping system. Each of 
these provides to one degree or another a voice of 
conscience for the designer, checking his design 
against various rules of fan-in, fan-out, impedance 
matching and the like. This, inddentaUy, has 
proved to be one of the most valuable of the ser- 
vices provided by the record-keeping type of pro- 
gram. I fed with a fair degree of confidence that 
the problems relating to software and standards 
will soon be solved to the point where we will have 
really fiexible engineering data-processing systems. 

What now about the design of large systems? 
There are some design problems which are common 
to nearly every large system, and some are unique 
to certain specific technologies. In this review we 
will look at a few of the common ones. 

Perhaps the most worked-over problem in large 
system design is that of placement. This might 
mean location of components on a printed circuit 
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board, location of boards on a frame, or even loca- 
tion of equipment frames in an office. Many at- 
tempts have been made to do placement entirely 
automatically. One technique developed at Sandia 
involves establishing artificial force fidds among the 
components. * Another developed at the Bell Lab- 
oratories involves a rough placement of dumps of 
components, followed by an iterative scheme of in- 
terchanging pairs of components and determining 
the resulting effect. I know of other techniques and 
I am sure there are many more of which I have no 
knowledge. Of one thing I am sure, however: 
that is, that none have been found to be entirely 
satisfactory in a practical system design situation. 

Let us observe that the previous paragraph dis- 
cussed a problem in automated design without men- 
tioning any criteria by which a given design can 
be evaluate. This is really the heart of the matter. 
The most easily stated criterion is minimum total 
wire length, and even in this case no completely 
automatic techniques have been found to solve prac- 
tical problems. But the problem is far more serious 
in that most design engineers are unable to specify 
all of their criteria at all clearly. In some cases 
crosstalk is important, so that there exists a limit 
on the distance over which a pair of wires can run 
parallel to each other. In other cases, capacitive 
loading or delay are critical, and here one wishes 
to minimize the length of the longest lead in a cer- 
tain class. Finally, esthetic considerations come 
into play, primarily as an aid to the maintenance 
man. If we have a structure (say a shift register) 
which has a definite electrical pattern repeated many 
times, it is desirable to have the corresponding phy- 
sical pattern repeated similarly. 

For two reasons, then, most of us have given 
up on fully automatic placement, and have turned 
to use of graphic consoles operated by the design 
engineer himself. Here the computer is able to 
keep accurate records and remember, and cause to 
be repeated, specific patterns; it can also evaluate a 
given layout with respect to any of a number of 
different criteria. We leave to the judgment of the 
man, then, those criteria he cannot really explain 
explicitly. 

Another closely related problem is one we refer 
to as treeing and routing of wires. Given the 
electrical schematic and the position of each com- 
ponent, we know which terminals must be connected 
together. In this instance, completely automated 
methods have been develop^ for determining the 
proper order of connecting the terminals to give 
minimum wire length, and also for finding a route 



* ''ACEL** Automated Circuit Card Etching Layout,” 

C. J. Fisk, D. D. Isett, Proceedings of the. 1964 SHARE 
Design Automation Workshop. 




for each wire such that no rules related to bunching 
of wires are violated. 

Now let's look at the early stages of designing 
a large system. To my knowledge, each effort here 
has been of a one-shot nature. It is not unusual to 
see two kinds of stimulation used in the early stages 
of development of a project. The first of these is 
a fairly gross overall representation of the sjrstem 
to allow the designers to find out where the major 
flows of information will occur and which are likdy 
to be the most critical parts of the sjnrtem in deter- 
mining its ultimate performance. The second type 
of simulation is to take a small chunk of the system 
and simulate its behavior on a pulse-by-pulse basis. 
As far as I know, even our fastest computers are 
too slow to allow a complete, detailed simulation 
of a sensibly large system. Indeed, even if this were 
possible, it is hard to imagine what one would do 
with all of the numbers which would result. 

In the area of detailed logic design we again look 
to the graphic console to be of great use. The day 
is not far off when a designer can sit in front of a 
console, draw chunks of circuitry, asking the com- 
puter to repeat given chunks when a repetitive 
circuit comes up. The machine will then check 
the circuit against the design and manufacturing 
rules appropriate to the particular technology and 
tuck the information away for recall at any later 
time either by the designer, by one of his colleagues, 
or by programs used later on to generate the added 
information required to manufacture the system. 

4. CONCLUSION 

Computers have already become indispensable in 
the design of communications components and sys- 
tems. Today the uses are largely of a clerical na- 
turer—receiving information, massaging it gently, 
and printing it out in various forms. However, the 
problems are rapidly being solved to the point where 
we are looking to the design process itself. In some 
cases we are able to develop algorithms which solve 
a design problem completely. In most cases, on 
the other hand, we look to the development of pro- 
grams involving direct communication between the 
designer and the computer through on-line graphic 
consoles. 
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HART: I would like to have your opinion about 
knowing eomething about computing, hut not 
having any ccureee in it. Would you explain it a 
littlef 

BALDWIN : I think learning how to uee a eomputor 
is becoming easier and easier, I think that if a 
student were to spend a lot of time learning For- 
tran, Algol, Flugol, or what have you, and also 
the mechanism of coding, he would be wasting 
his time. There is a need for the high priests 
of computing who will supply the eventual soft- 
ware whereby the computer will talk peopWs 
languages; but a man who is going to be a designer 
needs to know very little about computers. What 
he needs to know is a lot more basic mathematics 
and physics, 

CONWAY: You, of course, are aware that almost 
all engineering students are now learning to solve 
problems with a computer and understand how to 
use it. My observation is that if they don*t learn it 
in school, once they get into industry it is very 
difficult for them to learn it, 

BALDWIN: I would like to address myself to two 
points. First, today you need to teach some pro- 
gramming because we have not reached the uto- 
pian situation where it is easy to use the com- 
puter, But, to me, the reason you teach pro- 
gramming is so that students can use the com- 
puters in their other courses to develop a good 
fundamental education- 

Now, my particular job is hiring people who 
will become these high priests and supposedly 
make the computer &isier to use, but w^re mak- 
ing it mare and more difficult. The last machine 
I really understood was the IBM 650, But whan 
I hire people, / da not care whether or not they 
have had any programming experience in schaol, 
and I run a programming department. What / 
want is good, fundamental mathematics and 
physics and a statement of interest in this par- 
ticular field, be&iuce people can pick up pro- 
gramming easily if they are going to be good at 
it. They can pick it up certainly in a year on the 
job, and welre hiring these people for 40 years, 
A second point is that schools have machines 
which are far behind the ones that are used in 
industry. So you get a fellow that says **Ah ha, 
I know how to program the Blubiac,** and this 
doesn*t do him any good at all because we*ve got 
something else, 

BRANIN : Vd like to make a comm&it on that point, 
1 feel that once an individual has learned to pro- 
gram in any machine language, or any langwige 
at all, he has mastered the fundamental principles. 
The translation to another language is nowhere 
near as difificult as learning the first process. 
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BALDWIN: Ma^he he leame to mistrust the eom^ 
puter. This is indeed a lesson that one must 
learn, 

LICKLIDER: Fd like to make m intermediate pro- 
posal because I agree with you that there shouldn*t 
he courses in computers which spend three weeks 
on binary notaHon and teach how to program 
in some particular machine language, I also 
agree that people can learn pretty quickly on 
the job, but I notice that the self-taught people 
miss something very important. They haven*t 
really learned to understand about formalisms 
or the syntax of computer languages, and they 
don*t have a set of the abstractions without which 
they cannot reaUy make any progress in the 
modem world of eompuiers. So I would argue 
that there is a course to be taught about com- 
puters; maybe it isn*t being taught anywhere yet, 
but there is some lore somewhere between com- 
puters and rnathemabics that would be awfully 
good to have, 

HART: It is very easy for a student to go through 
four years of engineering school and never solve 
a problem, / think one of the opportunities the 
computer provides is that the student, by getting 
involved in the sweat and tears of writing a 
program to wive a problem, develops the concept 
of the problemrsolving procedure, the idea of an 
algorithm, and / think exposure to this early in 
his educational career is extremely important, I 



don*t euppose that you must then continue with 
more haw-to-do-it courses, but give enough train- 
ing in programming so that a student can ap- 
l^owh a computer and solve a problem with some 
insight into what this involves, 

BALDWIN: One of my reactions is that he can 
learn the problemrsolving procedure, which is 
actually mostly problem formulation and not 
problem solving. Typically, you*re reformulating 
the problem as you solve it, whether you use a 
computer or not. But can*t he do this by standard 
analytical methods f You don*t think so, 

BRANIN: There is a strong point that Hart brings 
out, which is that the discipline of programming 
a problem forces you to loam it in more detail 
and more thoroughly than you would otherwise, 
so I would agree with him on that, 

SHANNON: With the introduction of our time- 
sharing system, we have moved in exactly that 
direction. We moved the teaching of computer 
programming into the freshman mathematics 
course, deleted a course in engineering, in **ln- 
troduction to Computation in Numerical Meth- 
ods,** and inserted a new course in **Field and 
Ehectricity and Magnetism,** This is the current 
direction in which we*re moving, 

BALDWIN: Well, that sounds good, / don*t be- 
lieve everything / say myself, / think certainly 
this is a very serious problem in education. 
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I. INTRODUCTION 

The intention of this paper is to discuss only one 
aspect of aeronautical design procedure which is 
influenced by the ^dstence of digital computers. 
This particular aspect is the analysis of structures. 
While limiting the technological scope of the dis- 
cussion, this topic provides an excellent illustration 
of the dramatic development of an analysis capa- 
bility to cope with a very wide range of problems. 
Such progress would have been impossible without 
the associated development of digital computation 
madiines. 

During the last decade the aeronautical industry 
has been mainly responsible for the introduction and 
establishment of finite methods. These methods 
have been denoted by various names depending 
upon whether they have inherently regarded forces 
or deformations as the unknowns of an analysis, 
or whether the method has been expressed soldy in 
terms of matrix algebra or has implied certain com- 
puting algorithms. These methods all possess one 
feature in conunon. That is, they depend upon the 
process of idealizing the structure into an assembly 
of finite elements which are analytically connected 
at nodes or boundaries. 

Examination of the literature will show that de- 
velopment of the finite element method is continuing 
vigorously both within industry and at universities. 
Effort is being directed along two main avenues 
of development: improvement of the element and 
improvement of the system that accepts the element. 
The former implies more precise representation of 
technological behavior and the latter implies more 
efficient computing procedures. 

The aircraft structural problem which can be 
said to deserve the strongest claim to having 
prompted the development of finite element methods 
is the multi-spar, multi-rib swept wing. This struc- 
ture possesses both complexity and a high degree of 
indeterminacy. A decade ago, such a problem was 
a diallenge to analysis ; today, finite element methods 
have made its solution commonplace. It is also 
evident from reports that much experience has been 
gained in applying such methods to other aircraft 
structural analysis problems. These applications 
of finite element methods reflect advantages over 
classical aircraft analysis procedures in that they 
reduce considerably the weight penalties associated 
with making assumptions in idealization and that 
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they are much more capable of producing refined 
descriptions of complex structural bdiavior. This 
paper will describe several typical applications and, 
where possible, show the corrriation between test 
results and analysis. 

It is believed that there will continue to be pro- 
gress made in developing and applying finite element 
methods of structural analysis. It is also believed 
that considerable progress will be observed in what 
might be regarded as a lateral direction, that is, 
the application of such methods to other technolo- 
gical areas. This trend can already be evidenced 
with the molecular vibration analysis by Kaufman 
and Kaufman, the seepage analysis by Zienldo- 
wicz» and temperature analysis by Visser“ and 
Hampton. To lend weight to this contention, 
two applications of relevance to aircraft design 
are included in this paper. 

Finally it can be said that, just as the finite ele- 
ment methods arc dependent upon the digital com- 
puter machines for their existence, so are they 
also dependent upon the appropriate education of 
engineers for their correct application, interpreta- 
tion and development. The overwhelming response 
to related courses being offered within the industry 
is a clear indication of the need for such education. 



II. FINITE ELEMENT METHODS 
It is not the intention of this paper to discuss 
finite element methods with any thoroughness, since 
this has been successfully accomplished elsewhere 
in the literature.^* The objective is to dwell 
upon applications and thus stimulate interest in the 
methods. However, some statements concerning 
concepts should prove beneficial within this context 
Before proceeding to a finite element analysis, an 
engineer must have at his command within the com- 
puter capability an array of mathematically defined 
structural elements. The force-deflection relation- 
ships of each element must be explicitly expressed. 
An array of typical linear isothermal elements is 
shown in Figure 1. The existence of such elements 
will enable an engineer to idealize a structure into 
an assembly of such elements. An idealization 
is shown conceptually in Figure 2 in which the plan 
view of an aircraft swept wing is covered by an 
analysis grid. The grid, in turn, implies the loca- 
tion of analysis nodes and finite elements. 

The forces and deflections of the elements are 
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FIG. 1 TYPICAL UNEAR ISOTHERMAL ELEMENTS 




FIG. 2 IDEAUZATION INTO FINITE ELEMENTS 

then expressed in a set of linear algebraic equations 
such that equilibrium and compatibility are satisfied 
at the nodes and boundaries of the structure. For 
a large complex structure, this leads to a large num- 
ber of simultaneous equations. In realistic prob- 
lems, many hundreds of equations can be involved. 
Thus it be^mes clear that for an engineer to merely 
discuss his problem he must resort to the compact 
notation of matrix algebra, but that to solve it, he 
must rely upon the availability of a digital computer. 

m. APPUCATIONS OF THE FINITE ELEBfENT 
METHODS 

1. Airplane Design Studies 

Extensive application of finite element methods 
occurred during the design of the Boeing Model 787 
Transport Airplane. A representative selection of 
the studies is shown in Figure 8 and their pertinent 
features will be described. In all cases, the com- 
putations were automatic and the engineering effort 
arose in connection with the interpretation of re- 
sults and preparation of data for idealization. It is 
in the two latter phases that engineering judgment 
must be exercised if an anabrsis is to be successful. 

Referring to Figure 8, Region 1 represents an 
analysis of the complex wing-fuselage intersection 
structure. This area, with its conglomeration of 
bulkheads, flexible frames, cut-outs, fuselage cross- 
section transitions and wing-fusdage interaction 



FIG. 3 BOEING MODEL 737 DESIGN STUDIES 




structure, contains a formidable challenge to classi- 
cal analysis procedures even when supported with 
empyricd data. Region 6 has similar features to 
those possessed by 1. Region 2, 4 and 5 are wing 
structure and inherently contain an analysis problem 
due to their being swept. For a two-spar structure 
this is the shear transfer of load from the front to 
the rear spar and the amount can not be determined 
by classical methods unless supported by empyrical 
data. Figure 4, a region on the Boeing Model 727 




FIG. 4 WING ROOT SPANWISE STRESSES 



main wing, shows the excellent agreement between 
stresses obtained by finite element analysis and 
stresses obtained during static tests on the airplane. 
Regions 8, 7 and 8 of Figure 8 represent cut-out 
analyses around doorways, windows and access 
panels respectively. 

One major consideration in such analyses are the 
boundary conditions for the regions of analysis. 
These are usually obtained from a prior gross an- 
alysis or by invoking Saint Venants Principle. Re- 
gion 9 involves ducting around the engines and 
requires consideration of thermal effects. Region 
10 represents an analysis that is smaU by compari- 
son with those previously mentioned and yet con- 
tains enough indeterminate elements to justify the. 
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use of a finite element method. The problem con- 
cerns interactions of floor beams and door sills. 

In order to give some meaning to what might 
currently be considered a Uirge analysis, one of the 
preliminary design studies performed for the Boeing 
version of the C6A Airplane is shown in Figure 6. 




This region of fuselage includes large access open- 
ings for loading vehicles, etc. The analysis grid is 
indicated in Figure 6 and the details that describe 
the size of the problem are included in the flgure. 




2. Non-Linear Material Behavior 

In assessing the fatigue and fail-safe characteris- 
tics of an airplane component, a preliminary phase 
involves the determination of stresses that are dis- 
tributed around an existent crack. A typical com- 
ponent to receive such attention is a panel comprised 
of cover plate and stiffeners. It is usually assumed 
that in the middle of the panel there exists a crack 
of finite length and zero width. In other words, 
the panel possesses a structural discontinuity addi- 
tional to those caused by its attachment to stiffeners. 
To solve the stress distribution, flnite element meth- 
ods are used. 

It is known that for signiflcant loadings on the 
panel, the concentrations will cause the stress within 
a certaifl region to exceed the elastic limit. Thus non- 
linear material behaviors have to be represented. 
These can easily be read into the computer in the 
form of a table that rdates stresses to strains for 



particular values. A typical curve is shown in Fig- 
ure 7. The solution, which follows the approach 




outlined by Argyris, uses an initial strain con- 
cept, and involves both iteration and inversion pro- 
cesses. 




An example of a simple problem is shown in 
analysis grid form in Figure 8. The plate is uni- 
form, loaded by a uniform stress along the edges 
paralld to the crack and, since it is doubly sym- 
metric, only one quarter need be analysed. The 
distribution of grid lines reflects an anticipation 
of the location of the stress concentrations. The 
grid lines deflne the positions of flange elements and 
the rectangular spaces imply uniform shear flelds. 
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In Figure 9 can be seen the results due to a cer- 
tain magnitude of applied stress. The stresses 
shown are those acting normal to the crack at the 
line co-axial with the crack (Y=0). Assuming 
first a fully elastic material, there can be a very 
high magnitude of stress near the crack tip. This 
corresponds to the singular point that would occur 
in an analsrtical solution. Allowing the material to 
have non-linear properties shows a reduction in 
peak magnitude^ yielding and general redistribu- 
tion. 

3. Engine Cowling Anabrsis 
The engine cowling structure was analysed by a 
finite element method and certain critical defiections 
obtained. The cowling was a tube constructed from 
titanium honeycomb panels and the tube was loaded 




by internal pressure. Analysis was made complex 
by virtue of the irregular ii^ape of the tube cross- 
section, the presence of plates (splitters) connect- 
ing opposite sides of the tube, and the need to 
represent the bending behavior of the honeycomb 
{Mmels. A view of the cowling is shown in Figure 




no. 11 COWLING ANALYSIS GRID AND PRESSURE LOADING 



10 and part of the analysis grid and the pressure 
loading is shofwn in Figure 11. 

The finite element method used to analyse the 
cowling structure included a triangular plate ele- 
ment that possessed bending and membrane prop- 
erties. This was used to represent the honeycomb 
panels in the idealization. The computing details 
of this analysis, which produced both stress and 
defiection distributions and an influence matrix, 
are as follows: 

Number of nodes =864 

Number of structural elements =401 
Computing time =46 minutes 

(IBM 7094 Mk II) 

Te^ were performed upon the cowling and 
deflection measurements obtained. A set of these 
results is shown in Figure 12 for one cross-section 
of the structure. It can be seen that , the 




IV. THERMAL AND ACOUSTICAL APPLICA- 
TIONS 

The concepts expressed previously with regard 
to structural analysis by finite element methods 
are valid also for analysis in other technologies. 
In the matter of analysing temperature distribu- 
tions, it can be shown that for linear conditions 
a direct analogy exists between the formulation of 
the structural and thermal equations. Thus one 
solution program could solve both classes of prob- 
lems, provided that the appropriate data was sup- 
plied. In the acoustical problem, a fluid element 
is treated as though it were a structural element 
and is therefore easily incorporated into the struc- 
tural solution. 

1. Thermal Analysis 

Using a variational approach, Visser “ has shown 
that a close analogy exists between the solution of 
conduction problms in thermal analysis and the 
solution of structural problems. Thus it is pos- 
sible to conceive of finite elements representing 
assumed temperature behavior and expressed math- 
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ematically by the corresponding heat flow-tempera- 
ture relationships. In matrix form tiiis is denoted : 
Where ®=0.T 

S=Golunm matrix of nodal heat flows 
T= Column matrix of nodal temperatures 
0= Square matrix of influence coefficients 
for .a thmmal element 

In the analogous structural equations, the mat- 
rix 0 would correspond to the stiffness matrix of 
a structural dement. 

Visser explained the 0 matrix for a flange and 
the 0 matrix for a triangular plate is derived in 
Appendix A. Using similar approaches, mathe- 
matical expressions for a. comprehensive array of 
idealized thermal elements can be constructed. This 
array could m^tch the ^^pes. of elenfents used in 
a structural analysis. . ; 

r If it is possible to idedize a thermal problem 
into a set of mathematically deflned finite dements, 
if i^^niaiiu for a solution procedure to assemble and 
. ’spiv^'fof the unknowns. The method used can be 
/Jf.;j^fical;''to i^at used in solving a; structural prob- 
* ^ disctissed in detail here. 

appl^tions of thermal conduction analysis 
* using 'the* .triangular plate element will be de- 
scribed. The flrst is a demonstration of Saint 
Venants Principle. Irregular distributions of 
temperature were imposed on the boundary of a 
plate that was ^^ng** in one direction compared 
to the other direction. The grid used is shown in 
Figure 13 and the boundary condition of the three 




remaining edges is such that there is zero heat 
flow across them. Interior nodes also have zero 
resultant heat flow. 

The results of imposing two sets of temperature 
conditions are shown in Figure 14. Plots of tem- 
perature distributions across the width of the 
plate are given at various distances (fl) from the 
application edge. It can be seen that the irregu- 
larities have smoothed out to give uniform average 
values within two width-distances of application. 



APPLIED 

POIMT 

30 



30 



I- .23 .3 .73 1.00 1.23 1.3 1.73 2.0 
AISOLUTE VALUES 

C- .23 ;3 .73 1.0 1.23 1.3 1.75 2.0 

JHII 

DIFFERENCE FROM AVERAGE VALUES (20) 




no. 14 SAINT VENANTS nUNCIPLE - TEMPERATURE DISTRIBUTIONS 

The S6v^x>nd application is shown in Figure 16 
and was performed to test corrdation with another 
metiiod. The exact analysis invdved a solution 
to the heat-conducti<m partial-differential equation 
and the matrix analysis used finite elements., . 
linear temperature distribution was impoi^ «^ng 
one edge, zero temperatures maintaii^ed j^|i^j|f ‘ tw 
others and zero heat flow along the remi^hg fdge* 
Plots of temperature distribution show ih^t 
relation between the two methods was" exodfent. 




2. Cavity Analysis 

The cavity illustrated in Figure 16 represents 
a problem in acoustics. The flgure shows a two- 
dimensional view of a rectangular air cavity 
bounded on flve sides by rigid surfaces and closed 
on the sixth side by a flexible panel. The problem 
is to determine the response of the combined panel 
and enclosed air volume to an impinging «u».ftii fftical 
environment. The representation of the pand in- 
volves st an da r d structural idealization procedures, 
and to represent the air system within the cavity 
a fluid finite element can be used. 

To test the quality of the representation, a flnite 
element analysis was performed for an air system 
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FIG. 16 CAVITY AHALTOS USING fLUlD ELEMENTS 

surrounded entirely by rigid surfaces. Boundary 
conditions such that motion normal to a sur- 
face was permitted. Submitting the assembled air 
system matrices to a vibrati<m • ahalysis yidded 
mode? and frequences that correspond to standing 
wav?| in the cadty. In^ Figure 17, results for the 
frequencies of the finite element approach are com- 




FIG. 17 CAVITY MODES AND FREQUENCIES 

pared with those obtained by classical analysis.^ 
The two sets agree very closely and mode shapes 
were virtually identical for the two methods. 

It can be seen that the ideas epcpressed here are 
capable of considerable extension, such as applica- 
tion to more complex configurations and to repre- 
- sentations of fiuids other than air. 
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FENVES: I have two Queetione. First, how does 
the engineer describe an ind^tertninate structure 
eoTwisting of 5,000 nodes or whatever the example 
was? Second, how can all this data be input 
with any relidbUity, considering the number of 
input cards that must be punched? 

BAMPTON: That is a very good question and / 
know one or two engineers with bloodshot eyes 
who would sympathize with it. For one thing, 
the program that resulted in this analysis involve 
ing 5,000 nodes ‘presumed a certain configuration 
so that the preparation of data was eapahle of 
being auto‘mated. Now this of course won*t gen~ 
erally be the situation. But there was a certain 
prescribed configuration built into the program, 
namely, that particular problem happened to be a 
fuselage problom and it was known that it would 
consist of grids defined by stringer lines a'nd 














frames, so this presented a lattice upon which 
the engineer could hang hie data. This eased 
the problem which I think you implied by your 
question, namely, how can any engineer or team 
of engineers handle so much data. 

FENVES: My next question is a little bit more 
fundamental and has something to do with the 
discussion we had earlier concerning the theoretic 
cal approach vs. the algorithmic approach. The 
procedu/re which you discussed was a very valid 
appUeation of computers, that is, you take a 
continuoue elliptic differential equation problem, 
chop it up into pieces so that you have only small 
pieces to concern yourself with, and then you use 
finite mathemaUcs to combine these pieces to- 
gether. Now you spend a tremendous amount of 
research on describing these particular elements. 
On the other hand, the technique, the algorithm, 
for interconnecting this very large number of 
elements so as to form semething that can be 
computed, you dismissed in one sentence. You 
simply said, **We obviously don*t use matrix alge- 
bra.** Why is there this disparity of emphasis 
and of prestige? There appears to be a great 
amsunt of prestige associated with deriving new 
finite elements, but all the work involved in try- 
ing to put together a couple of thousand finite 
elements, so that todays computers can get a 
solution rather than waiting tdl the seventh gerv- 
eration mashines are big enough and fast enough, 
is tossed away as a coding tricki 

HAMPTON : I don*t know what is the basis for your 
assumption. I*m not tossing aside the problems 
involved in the system which is accepting the 
elements, and I didn*t say that obviously we 
can*t manipulate matrices in a tone which im- 
plied that this was a trivial matter. I just hap- 
pened to have picked up a crutch, y<m might say, 
for discussianr-mamaly, the finite element. There 
is na doubt about it; there is a lot of work going 
on in examining the system which incorporates 
these finite elements in order to handle them in 
the most efficient manner. 

SCHMIT: I would like to ask if you cotdd teU us 
anything about what effect, if any, the computer 
is having in your organization, or in organizations 
you are familiar with, on arriving at the system 
you had when y<m set about making this idealiza- 
tion; in other words, you described an aft end 
of a fuselage and you had something that you 
were ready to analyze. Certainly a great deal 
of preliminary design or some sort of effort went 
into bringing you to that point. Are computers 
having any impact on this sort of thing— the 
preliminary design phase? What Vm really con- 
cerned with is this: are we faced with the prob- 



lem of optimizing what may not turn out to be 
a good idea, or attempting to get an efficient de- 
sign based on a preliminary design which is not 
necessardy the best of alternatives? Is there 
any activity in using computers in this prelimi- 
nary design stage in your organization? 

HAMPTON: Yes, there is, and the activity is di- 
rected partially toward structural representation, 
with which I*m most involved, and also toward 
what might be called integrated analysis. That 
is, you have a structures program, a weight pro- 
gram, and aerodynamics program, etc., and you 
attempt to integrate these in order to get an 
efficient preliminary design. This has only been 
possible on a limited basis because you*re faced 
with the task of preparing data. Thus, certain 
configurations have been assumed or buUt into 
the integrated design program and they do not 
reflect a structure that can be sophisticate as 
our structural program can produce for us. 

CHASE: rd like to answer the question, a little bit 
more. Actually there is a great deal of work 
going on in this area of the borne process. We 
have at Lockheed, for example, a bank of com- 
puter-system programs dedicated to ' the man- 
computer interaction that is being discussed to- 
day. In particular, one of the major areas of 
investigation involves the prelimiary design of 
aircraft in general, the surfaces, etc., in which 
structures of various sorts will be designed at 
some later stage. It is in its infancy. There*s 
a great deal of work being put into it. Definition 
of surfaces is a complex problem within itself, 
but it is being conducted in industry. Some of 
the work is done at MIT under Coons in the com- 
puter-aided design program. But the whole e<m- 
eept of design in its many phases and many 
aspects is being investigated and we^re whittling 
away at little pieces at a time. I think that you 
wiU see within a matter of two or three years 
quite a large accomplishment in which all these 
features will be put together in a somewhat more 
sophisticated package. 

PAYNTER: I*d like to propose a question today 
which I*U have a chance to ask again tomorrow, 
but which is particularly appropriate to the pres- 
ent speaker. So far we have not yet corns to grips 
with the problem of how can you evaluate whether 
a model is really any good or not. We saw here 
comparisons of a point-for-point type response 
between some test data and some model predic- 
tions. I*d like to raise a question in line with 
the comment just made in Fenves* second ques- 
tion: whether the following problem ought to 
be posed, that is, how does the presumed affect 
the decision process itself? In the context of 
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many of your oituatimo, ypu aro not going to 
design pUUoo Vfith variable ihiekneeeee at differ'- 
ent points on the airerdft shell, so that, if yon 
are restricted to using plates of uniform thick- 
ness, isn*t it true that as long as the model pre- 
dicts within one gangs thiekkess, yon would ac- 
cept it as an acceptable representaUonf 

BAIIPTON: / think the answer to that is in the 
affirmative, and I can only comment further that 
we have been investigating optimisation pro- 
grams using linear programming where such con- 
straints as standard plate gauges have been 
entered into the programming problem. 



SAUNPERS: / would like to inquire about the 
make-up of the teams Wha are doing this work. 
Are they pred/ominantiy of an engineering back- 
ground, or mathematics, or compute people, or 
fuel who are they f 

HAMPTON: Iffs a combination of engineers and 
mathematicians with the main emphasis on engi- 
neering, At Boeing, there is a distinct dividing 
Hne between mathematicians and engineers and 
they have to eomrnnnicate across that line, . Ocea- 
sionaUy, mathematiciaiw do cross that line, but 
the engineers are not allowed to do so. They co- 
operate pretty closely, and it is a combined effort 
that produces these kinds of formulations. 



APPLICATION OF COMPUTERS TO POWER 

SYSTEM ENGINEERING 

G. W. Stagg and a. H. Palmer 
American Electric Power Co., 

New York, N. Y. 



This paper will describe the application of com- 
puters in power system engineering and operation. 
It is essentially a case history of engineering com- 
puter applications in the American Electric Power 
Company. 

The American Electric Power Ssrstem provides 
electric service in a seven-state area. It has a 
generating capacity of over 8 million kilowatts, in- 
stalled in 15 plants with 38 generating units. Power 
is distributed over 81,000 miles of transmission and 
distributing facilities. The Ssrstem is connected 
to 21 neighboring power systems by 46 tie lines. 
In 1966 the peak load was slightly over 7 million 
kilowatts, of which approximately 99% was gen- 
erated by steam units requiring 16 million tons of 
coal. 

This particular case history demonstrates a more 
genera] operation than many others in that it in- 
volves electrical, civil, and mechanical engineering, 
as well as engineering information and project 
control applications. TTie achievements made in 
these areas will be described and some conclusions 
will be drawn concerning what the future portends. 

Initially mention should be made of the back- 
ground and resources which have been brought to 
bear on the problem. The word problem is used 
here in reference to provision to the practicing 
engineer of the computer capability of the digital 
computer, without requiring him to perform the 
task of program development. There exists a cen- 
tralized group which is responsible for the develop- 
ment of computer applications. 

For the past ten years, AEP has been requiring 
and training a staff of professional people. At pre- 
sent, the Engineering Analysis and Computer Divi- 
sion consists of forty engineers, analysts and 
programmers. In addition, a supporting group of 
ten people performs the computer operation, key- 
punching and computer services. 

Two major groups operate with the Division. 
The Engineering Computer Applications Section is 
concerned with the development of programs for 
use in electrical, civil, and mechanical engineering 
activities. A prime responsibility lies in the evalua- 
tion and development of techniques suitable for 
computer implementation. The applications are 
generally large-scale programs which consider the 
entire problem area and are intended for system 
design and simulation. 



The Engineering Information Processing Section 
develops programming systems for collection, pro- 
' cessing, recording and analysis of the company's 
operation infomuition in order to produce reports 
necessary for engineering control and evaluation, 
and also to satisfy regulatory r^uirements. 

In each of these groups a ^tinction is made 
between engineer and programmer-anal 3 rst. The 
programmer, or analyst, generally provides a high 
level of mathematical training, whereas the en- 
gineer provides an extensive background in the 
analytical aspects of the problem area. There are 
forty baccalaureate degrees and thirteen advanced 
degrees among the forty people in the two Sections. 
It should be stressed here that there exist varied 
educational backgrounds within the division. In- 
cluded are electrical, mechanical, civil, chemical 
and industrial engineers, mathematicians, physicists 
and business administration personnd. 

The staff has grown i>artially by obtaining per- 
sonnel from outside the organization, although of 
late some emphasis has been placed on recruiting 
from within, particularly in the Engineering Infor- 
mation Processing Section where experience within 
the organization is invaluable. 

Mention has been made of the resources which 
have been used. An important facet has been the 
availability of computing at reasonable cost. Al- 
though the units of computer usage have increased 
considerably, the cost through the years has grown 
at a much dower rate. 

A significant step in the continuing program of 
computer utilization for the solution of complex 
engineering problems occurred last year. On Nov. 
1, 1966, an IBM/360 Model 40, specifically intended 
for engineering, was installed in the Service Cor- 
poration Offices in New York. This sjrstem is 
scheduled to be replaced by a Model 60 in December, 
1966. It is expected that these machines will pro- 
vide an impetus to additional expansion in the 
computer application area. 

The efforts of this staff will be the concern of the 
remainder of this paper. Several broad subdivi- 
sions of application areas immediately arise. Al- 
though one of these is identified as design, all 
evidence a relationship to the design process to 
some extent. In the context which follows, design 
is to be considered the planned change of the AEP 
system. This is a continuing, complex, multi-step 
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process \vhich has as its turn the expansion and 
improvement of the systeuL 

The several areas which will be indicated by ex- 
ample are: 

1. Design. 

2. Simulation. 

3. Control. 

4. Information processing. 

DESIGN 

Among the most promising engineering applica- 
tions, from an economic viewpoint, are those in the 
design area. Programs of this tsrpe provide design 
assistance in connection with power-plant steel and 
concrete work and in location and design of trans- 
mission towers. In addition to material and design 
time-saving, the means are available for a more 
comprehensive analysis of each problem. 

The power-plant design function is supported by 
a system of programs developed to aid the design 
engineer in achieving an economical design within 
time limitations associated with this type of work. 
The Structural Steel Framing Program is used in 
the design of plants of braced frame construction. 
In the design of a plant of this type, the flooring 
system, made up of rolled sections and welded plate 
girders, is designed somewhat independently of the 
bracing and colunm systems. 

Data prepared for the flooring system is pro- 
cessed with or without specification of col umn re- 
quirements. As the floor design is devdoped, the 
reactions contributed to columns are automatically 
retained for column design. 'Ihe normal sequence 
of events on the computer provides the designer 
with all flooring member sizings and, when required, 
column sizes. 

An interesting extension to this program is the 
production of detail fabrication drawings. The 
information necessary for this step is all available 
within the design programs, and the technical and 
economic feasibility has been demonstrated of sub- 
sequently processing a magnetic tape file to produce 
drawings. The high-speed microfilm recorder dis- 
plays information on a cathode-ray tube which 
is film-recorded to 3 deld a drawing. 

The Foundation Slab Analysis Program currently 
being used considers the mat to be a rigid body 
developing a planar soil-bearing pressure distribu- 
tion, an assumption which is valid for the heavy 
slabs generally used in power-plant work. 

The utility of this program lies in the possibility 
of evaluating a number of alternative designs and 
considering the vast number of load combinations. 
Allhough this is a fairly simple approach to the 
general problem of mat foundations, considerable 
work has been done in the investigation of a pro- 
gram to design the mat based <m tfain-plate bending 




theory for a plate on an elastic subgrade. A third 
program is also under development which considers 
the same foundation slab on piles. 

In transmission line design, the computer enters 
the problem area at an early stage. Once a line 
plan is selected, field survey crews obtain terrain 
data. The field survey notebooks provide the data 
for the terrain model. It is hoped that eventually 
this data may be obtained even more directly by^ 
means of aerial photographs and subsequent auto- 
matic digitization. 

In addition to the terrain data, it is necessary 
that constraints be specified for line crossings, con- 
ductor clearances and restrictions on tower posi- 
tioning. The Transmission Tower Location and 
Selection Program uses a d 3 mamic programming 
technique to achieve an optimum economic ar- 
rangement of tower locations and tsrpes. Several 
additional supporting programs available for the 
transmission area are sag and tension analyses of 
conductors, both for the design and field-stringing 
conditions. 

Two salient factors should be noted. First, an 
integrated system of programs provides the capa- 
bility of operating in the entire range of the prob- 
lem area. Second, this system of programs, to be 
effective, must be inte^ated into the day-to-day 
work load. 

A number of other special-purpose programs are 
available for use in civil engineering functions; for 
example, analysis of backwater curve, gravity dam, 
reservoir capadiy and concrete stack. 

SIMULATION 

The planning and design of a large complex power 
system is virtually impossible without a computer. 
In the simulation of such a system it is necessary 
to develop a mathematical model which can be used 
for stud^ng changes in the real system. By the 
use of this modd, any number of ssrstem conditions 
can be simulated and their effects observed. These 
applications run the gamut from steady-state elec- 
trical network analysis to transient analysis. 

A steady-state model of the dectric power system 
is the load flow. This program evaluates station 
conditions and power flows and considers the volt- 
age-regulating capability of generators, s3mchronous 
condensers and tap-changing-under-load trans- 
formers. In addition, desired energy interchange 
between systems is considered. 

The system is described in terms of network 
connections, transmission line impedances, trans- 
former impedancfis and capacitor data. For a given 
study year, loads, generation, voltage levds and 
interconnection data are given. The study results 
include a detailed report on the status of the sys- 
tem during the study period. 
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This program permits the planning engineer to 
consider a wide variety of cases. The effect of 
future load, generation, transmission lines, trans- 
formers, interconnections and voltage control equip- 
ment may be evaluated. Furthermore, the outage 
of various components can be simulated. 

This iype of program can be used to analjrze vast 
systems. This particular version of the IBM/860 
can analyze systems with as many as 2,000 stations 
and 4,000 lines, and, with random storage, systems 
many times this size. Until the general availability 
of digital (^mputers this was not possible, nor were 
the solution techniques sufficiently well defined. 
Today, the design engineer has available the capa> 
cily to study in detail large interconnected systems, 
and thereby plan better future systems. 

The Transient Stability Analysis program simu- 
lates behavior of the electrical network during 
switching operations and fault conditions. The dif- 
ferential equations describing the machines are 
solved numerically to determine the system condi- 
tions, following the specified operations. Of parti- 
cular interest is the response of generators during 
the sequence of events. The aim in using this pro- 
gram is to provide data for use in the design of 
electrical facilities. 

The Switching and Lightning Surge program also 
evaluates transient currents and voltages. The time 
scale of interest here is of a much shorter duration 
than in the transient stability analysis. In effect, 
a voltage wave is simulated in order to determine 
design criteria for lightning arresters and switching 
equipment In addition, insulation requirements 
for electric facilities can be determined. 

In the mechanical engineering area a program 
is available to evaluate a steam power-plant <ycle. 
A one-line diagram defining the cycle is the starting 
point for a computer solution. The cycle is simu- 
lated by interrelating the various components of the 
(ycle. A subprogram utilizing the characteristics 
of a component as defined by input data can simulate 
performance of a segment of the cycle. Subse- 
quently, the remaining components are evaluated 
and eventually this iterative process converges. 
Results are provided for various points on the cycle. 
This program is used to design new cycles, evaluate 
current performance and to simulate modifications. 

CONTROL (Operational) 

The allocation of power-plant generation to effect 
optimum operating economy is an important prob- 
lem in an electric utility system. The complexity 
of the problem for a large, integrated power system 
induced American Electric Power Company to in- 
stall, in 1964, an automatic real-time control system 
to provide optimum operation of its power facilities. 
The facilities to provide on-line control consist of 



a 1,700-mile microwave system, analog and digital 
telemetering, analog equipment and a digital con- 
trol computer. 

The digital control system provides, on a real-time 
basis, economic loading schedules for 38 generating 
units and, in addition, performs power-system oper- 
ating studies and logging functions on a time-shared 
basis. The control system has proved to be effective 
in allocating generation for optimum economy in 
AEP System operation as well as affording other 
savings. 

The functions of the Digital Control System are: 

1. Economic dispatch under normal conditions. 

2. Study mode to permit investigation of 
energy interchange. 

3. Logging. 

INFORMATION PROCESSING 

The logging function associated with the Canton 
control computer is a part of the AEP information 
collection system. 

The MHW (megawatt hour) Data System estab- 
lishes a history file from the hourly collection of 
generation, interchange and load data. Hourly and 
daily reports are produced automatically from this 
file and are retained in conjunction with regulatory 
requirements. 

SUMMARY 

In summary it should be indicated what plans exist 
for continuing the efforts described. The primary 
objective is to give to the practicing engineer caught 
up in the swirl of details a more effective path to 
computing capability. This is not to be construed as 
a path to the raw computer, as much as a channel 
of communications with the programs which are 
appropriate to the special problem. This entails a 
three-step process: 

1. Development of integrated sets of computer 
applications associated with each problem area. 
These are well under way. 

2. Development of a data base upon which these 
programs draw. Some inroads have been 
made in this area. 

3. Development of a program to interrdate the 
problem program, the data base, the engineer 
and the computer. This presents a significant 
obstacle which at this point remains to be 
passed. 



PAYNTER: It looks to ms os if yoi&*ve outlifisd o> 
system with a great many fragmented programs. 
You spoke of integrating these in a system of pro~ 
grams hut how will Miller* s master designer fit 
into your scheme of things? 

PALMER: TYsU, 7 don*t see that this is at all in con- 
tradiction to Miller* s approach. / think this is 
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what we are heading toward, ^erhaph ih^ 
which we travel to gOb there It eomewh^ dif* 
ferent from hie, hut I think baeieaUy we are 
coming to that came point. 

LICKLIDER: Are there any parte of. the epectrum 
in power dietribution that are interfereneodree 
enough to coneider eendiing eignale over them, 
and, if eo, are there any thoughte about eending 
computer eignale over your power Uneef And 
if not, are the probleme legal or teehniealf 



PALMER: In my experience the only time we be- 
came involved with eomething of that eort woe 
eeverol yeare ago when the Civil Defence people 



wercthink^ in termeOf doing eomething of that 
eort, don^t think much work hoe 

bem dom, Pm noi^reaXly qualified to epeak on 
thie matter, 

OFP MIKE: Would you care to dal)- 
ordieuponmatf 

H think thie idea woe 

pretty muck the deeieUm to go to a 

microwave eyetem. But, in anewer to the ques- 
tion, it ie now being used by a number of com- 
paniee for control eyeteme, particularly in the 
far west where a powerline has been in use for 
twenty to thirty years. 



APPLICATIONS OF COMPUTERS IN HIGHWAY 

ENGINEERING PRACTICE 

Theodobe F. Mobf 
I llinois Division of Highways 
Springfidd, 111. 



There is an old sasring that has come down to us 
to the effect that **the cobbler’s children go bare- 
foot.” As in the case of many other enduring ex- 
pressions of traditional wisdom, this saying survives 
because it contains a large core of truth. 

In a certain sense this saying has been true of 
the techniques of the engineering profession. Engi- 
neers, as a profession, are a very creative group — 
perhaps the most creative of all the professional 
groups that there are. Yet this creativi^ is 
directed outward toward the design of products 
and remarkably little change has been made within 
the engineer’s own shop in the techniques that he 
uses during the creative processes. 

It always gives me the greatest pleasure to study 
the exhibit cases of the antique engineering, astro- 
nomical, and measuring intruments on display in 
the Adler Planetarium and the Museum of Science 
and Industry, here in Chicago, or in the Smithsonian 
Institution in Washington, D. C. Many, if not most, 
of these exhibited instruments are represmited in a 
modem form in our present-day engineering offices. 

If we seek differences between these antique 
instruments and their modem equivalents, we see 
less jewdlike elegance and perhaps somewl^t more 
utilitarian precision in the new ones. But, by and 
large, a modem engineer, draftsman, or surveyor 
could use without difficulty the instruments of his 
centuries-ago predecessor. 

Indeed, I often think that Copernicus, who lived 
about 450 years ago, would have far less difficulty 
in accustoming himself to the instruments and tech- 
niques in a modem state highway department de- 
sign office, than his wife — ^Mrs. Nicholas Kopemik — 
would have in preparing a batch of Polish Christmas 
cookies in my wife’s modem kitchen. 

Until recently, the most complex mathematical 
devices available to engineers were the desk adding 
madiines and desk calculators. Far from being 
modem, we may be surprised to learn that completely 
workable models were constructed by Blaise Pascal, 
who died at the age of 89 in 1662, about 800 years 
ago. To add a personal note, Philip Matthew Hahn, 
a south German Evangdical pastor (and my great- 
great-great-grandfather) designed and built a cal- 
culating machine whidi was probably the prototype 
of the modem CURTA computer over 200 years 
ago. 

Standing as it does in the mainstream of scientific 



and mathematical thought, the engineering profes- 
sion has advanced on the break-throughs of the 
God-gifted geniuses of the past— Aristotle, Archi- 
medes, Galileo, Napier, Newton; and more recently, 
Einstein, Fer^, and many o^ers. We may all 
stand a little taller and hold our dioulders a little 
straighter when we acknowledge this heritage of 
genius, as it speaks to us today across the echoing 
chasms of the centuries. 

WMle the bold advances in conceptualizations 
which we owe to these geniuses of the past and the 
present have normally led to advances in technology 
as we have come to comprehend their meanings, 
we may today be seeing something of an inversion 
of this normal process in the adoption of the elec- 
tronic computer in modem uses. 

The electronic computer, as it has been initially 
used, has been simply a far, far better mousetrap 
than we have had before. In its initial applications 
it has been like a |16 hamburger-— an enormoudy 
advanced device. its means we have done our 
conventional computations in an unconventionally 
rapid manner. 

However, as we begin to appreciate the pos- 
sibilities of the device, we find ourselves pressing 
against the limitations of our convention-bound 
tediniques. We are indeed looking for new con- 
ceptual approaches to exploit, in presently unknown 
W{^ Hie vast potentialities ^t we have discovered 
in the binary nature of an electrical pulse when used 
in combination with the rigors of yes-no logic. We 
are beginning to see, in real-time applications of 
computer control as well as in strictly mathematical 
processes, the intimations of a brave new world of 
conceptualization made possible by technical ad- 
vances. 

In the fidd of highway engineering it was my 
privilege to have been among the first to be exposed 
to the possibilities of the electronic computer. This 
was in December of the year 1966, and the place 
was the annual meeting of the American Associa- 
tion of State Highway Officials in New Orleans. 
Until then I had been dimly aware of the work 
going on at Harvard and a few other places. I 
had stood in the blinking, oven-like presence of 
Illiac at the University of Illinois. This was the 
“Tom Swift and His Magic Brain” period of popular 
journalism. As I remember, at tUs time a Univac 
was helping out on one of the popular television 
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prognana by finding, through mysterious mathe- 
maticsl means, a magical matching of matrimonial 
mates or something else equally preposterous. 

Then at Mi a mi in Janua r y,*^ 1956, another manu- 
facturer stimulated thinking among migineers by an 
exhibit at the American Road Beiders Convention. 
My own first serious introduction to this equipment 
was in a small discussion group, following the 1966 
Highway Research Board meeting, attended by only 
three or four highway engineers in the ofilce of the 
late H. A. Radzikowski, at whidi his two assistants, 
S. E. Ridge and L. R. Schureman, described the 
results of their preliminary investigations into elec- 
tronic computers during the past srear. 

The first electronic computer conference, under the 
auspices of the Electronics Committee of AASHO, 
was hdd in Chicago in early March, 1966, followed 
by meetings in Atlanta, Los Angeles, and Boston. 
Within about a year, nearly every hi^way depart* 
mmit either had an electronic computer on the fioor 
or on order. And within two years the various 
highway departments and consulting engineers in 
the highway fidd had over 1,000 programs in use 
applicaUe to the solution of highway engineering 
problems. In the years which have passed since 
then, the scope of applications has become so broad 
that inventories of numbers have become meaning- 
less. 

Without intending to venture into the jungle of 
superlatives as to which siiecialiy, wi thin the whole 
engineering profession, was the earliest, or the 
fastest or the most efficient in appljdng the com- 
puter to engineering problems, I would say that the 
highw^ group as a whole has made a rmnarkable 
^rd in this respect I would not be able to (nor 
if I coul^ would I wish to) bore you with a catalog 
of the highway engineering problems which have 
been aided by the dectronic computer. These are 
all wdl known, and, if I could find a problem not 
yet within the range of analysis by omnputer, I 
would only be putting the spotlight of challenge on 
the next problmn to be solved, in the words of the 
mountain climber, Mallory, *1>ecause it's there." 

Certainly, the dectronic computer, used as a 
calculating device, has profoundly affected highway 
engineering practice and all of those who are in 
highway engineering are aware of this change in 
varying degrees. The response to this change has 
by no means been uniform among individual en- 
gineers, as they show reactions ranging from 
tremendous enthusiasm through stolid Aiy>^ptfmce 
and all the way to baleful skepticism; but not indif- 
ference; No thinking person can be indifferent 
whmi his mental cage is being so vigorously diaken. 

How are dectronic computers being used in high- 
way departmmits ? Of course, as in any other large 
business organization, there are fiscal, statistical 



and accounting uses of Uie equipment. However, 
since ^s conference is concerned with engineering 
education, I would like to confine my remarks, to 
engineering uses rather than accounting uses, while 
always acknowledging thdr presence. 

Apart from conventional data processing, the 
uses in technical applicatimis are quite broad, rang- 
ing from information storage and retrieval through 
the fidds of traffic assignments, traverse computa- 
tions, structural computations, and soil mechanics 
if one might suppose there to be some kind of syste- 
matic or orderly continuum. 

There is also a certain futility in describing how 
elecixonic computers are being used in highway 
engineering, either technically or organizationally, 
because changes are taking place so rapidly. 
Changes in the tjrpes of problems susceptible of 
solution, in the scope of programs, in the software 
which permits better operational scheduling of work, 
also more efficient, so-called small computers and 
more practicable remote operations on large com- 
puters are all factors which modulate the momen- 
tary advantages of any particular mode of use. 

This is why a discussion of what any organization 
is doing must be understood to be a kind of amalgam 
of the outmoded methods and obsolete equipment 
that is actually in use, as modified by the intentions 
of use for the equipment on order, and liberally 
seasoned by ideas on techniques now only dimly 
perceived for the remote future. 

Generally the use of electronic computation in 
highway engineering today is not in itself a system 
in the sense that data processors use the term, nor 
in the light of present knowledge is it likely to 
become, broadly, a S3rstem in the foreseeable future. 
Some parts might be termed thus, perhaps, but in 
the main, computations, whether done by pencil 
or by electronic computer, are only a part of a 
larger system of design. 

(generally also, in the present state of the art, 
highway engineering problems are mainly of the 
anabrtical rather than of the design ^jrpe. In the 
analjrtical situation, the designer sets certain con- 
ditions and the computer is programmed to calcu- 
late certain results or to optimize certain aspects. 
The designer reviews these results and, if he is not 
satisfied, he will change the variables and repeat 
the analsrsis. 

This feature of a dialogue between the engineer- 
ing designer and the computer often argues for the 
closest association between the two and against 
the hig centralized computer with its theoretical 
efficiencies for data processing. Functionally, an 
el^tronic computer in an engineering analysis use 
might be compared simply with a tremendously 
powerful desk calculator available, if not imroedi- 
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atdy at the designer's elbow, at least in a menage 
where scheduling is not a formidable problem. 

Anotiher growing area of use of computers is in 
the real-time application^ traffic surveillance and 
control. For somd>tiiTO highway engineers have 
thought , of the traffic usage of an urban freeway 
as bdng analogous to any other continuously operat- 
ing process. In a chemical plant, for instance, a 
sensing and feedback control of readings of pres- 
sure temperature and composition has miabled 
process automation through the operation of the 
servo values. Similarly a sensing of speed, volume, 
and densily of traffic on a freeway will enable the 
optimization of vehicle throughput, or speed, or 
gross time of v^icles in the transpoi^tion corridor 
by controlling the traffic signals at the entrance and 
exit ramps of the freeway. 

We now have the hardware and the knowledge 
for traffic controls of tills kind and are, in fact, 
doing it to some extent on the freeways in the Chi- 
. cage ar^ andjn a somewhat different manner in 
* Defooit, in river tunnels in New York, in Seattle, 
and in other places. 

In Baltimore New York, Toronto, and dsewhere, 
entiio grids of traffic signals on suifface streets are 
controlled by electronic computers on a real-time 
basis. 

Further uses of process-control computers are 
seen in continuous sensing in the production of 
concrete, asphalt aggregate and other materials, 
but here the problem is less one of control than 
of knowing how to sense the variables that it is 
desirable to control. 

What do we expect from the engineering schools 
in educating the engineers whom we will be hiring 
as graduates this summer and from now on? At the 
bachdor's levd, not too much is expected beyond 
a general familiarization. What I am saying is 
tha^ given a fixed number of classroom hours, we 
would rather have the schools maki iig our prospec- 
tive employees better engineers than poor computer 
programmers. 

At tills level there should be several points of 
emphasis that are more important than instruc- 
tions in how to manage a problem in any particular 
computer language: 

1. First and foremost is the appreciation of the 
rigors of binary logic as appUed to the solution 
of a problem. Every day we can see an in- 
crease in the evaporation of the velvet fog of 
"engineering judgement** and we are applying 
informed and discreet choices in more and 
more situations. 

2. A disenchantment with the magic of ignorance 
concerning these fabulous devices and the 
admowledgement tiiat tiie limitations on com- 
puters are mainly the limitations of our God- 



given ability as humans to understand and 
analyze a problem or process. 

3. Instruction in a pe^gogical climate which 
acknowledges that yesterday's brave new world 
is already here and takes for granted the 
ever-widening application of computer uses. 

As to specific computer instruction for practicing 
engineers of our staff, we have been conducting 
schools of instruction and would prefer to continue 
this practice. . Our reason for tUs is that the in- 
struction must be closely related to the specific 
computer or computers whidi are being used by our 
engineers, and ^e program languages are those 
which are applicable to our programs, and we want 
to continue having it this way. 

Such schooling has been given to more than 500 
of our engineers in both the central office and in 
the district field offices. A further advantage of 
having tills instruction done by engineers of oijr- 
computer staff is tiiat it establishes their authority 
in the field and also aids in communications bet^6«^n 
those providing and those using the compute'r-cen- 
terse^ces. " 

Beyond the bachelor's level there woul^, ^f course, 
be specific instruction in one or more of the common 
program languages, and this is taken for granted. 
Here^ however, the emphasis should be on the pos- 
sibilities of tile compute* to lead the profession into 
higher levds of technical capabililly through better 
means of analsrsis and fewer approximations, or 
rule-of-thumb solutions. 

As I reach the conclusion of my remarks I am 
struck by tiie paradox of some of the positions I 
have taken. I have written with deepest respect of 
the capabilities of the equipment and with some 
disdain as to their physical presence, or the lan- 
guage in which we must communicate with them. I 
am more appalled tiian impressed that some of them 
will generate seven acres of paper covered with 
printing in a few minutes. I condude that I will 
have to explain this paradox to myself by seeing 
these devices as the man-made wings— -wings unim- 
portant in themselves— that will give added lift 
to tiie soaring spirit of man to enable him to move 
toward the fulfillment of those ends of creativity 
for the benefit of mankind for which he was made 
in the image of God. 



LICKUDER: Two or throe times there have been 
eommente about the UkeUhood or unlikelihood 
of ever pMng together the myriad problem- 
oriented programs in a fMd into a system. This 
was something that was approached this morning 
by Miller and it seems as though we ought to 
get a Uttle better agreement about the ground 
ruies, / think there is at least one concept that 
doesn*t suppose that you wiU ever puU sufficient 
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number of programs together to form a system 
that will itsdf do design, if you mean a system 
of programs and computer hardware. But there 
is an aUemative which says one should couple 
with the programs enough language implementa- 
tion to make it easy for a person skilled in the 
art of the particular field to shape up, with the 
help of these already programmed resources, most 
of the software he needs to solve his particular 
problem. Now this coupling of the concept of 
language to the concept of a system of programs 
seems to me to he extremely basic and Pm 
thinking it would be important for this group to 
approach closer and closer to, an understanding 
of whafs possible there, 

MORF: I think that there are about 12 reasons, and 
there might he IS or 20, why this is a result that 
I eee as being rather distant, I don*t eay that it 
iepenMhievable,andinfactthereareanurnberof 
•ejfbrte being made in this direction. One of these 
, . .ts being sponsored in part by the Bureau of Public 
•• Roadfi^^^^Jt is called TIES: Total Engineering 

Integrmed System, Another similar effort ie be- 
ing made at MIT, eaUed ICES, There are prob- 
lems here and they are very deep. One of the 
r^one ie that we don*t quite know how to score 
the ball game, Lefe eay, if you try to work a 
matrix on how to maximize the profits of a com- 
pany, you have a rather simple objective. But 
what are we trying to do in the highway field? 
Does anybody know what the objective of high- 
way adminietration or a road system ie? 

MILLER: We had some exposure, to TIES a year 
or 80 ago when it woe first being formulated 
as an ideal. Our efforts are now devoted entirely 
to ICES, and this ie different from TIES as ie 
shown by the absence of the word Total, The 
ICES system, by decision, makes no pretense to be 
a total system. But it makes at least a start 
towards the goal of integroMng what in the high- 
way field, as an example, have been UteraUy hun- 
dreds of pieeemeal programs. And it ie one single 
system, a whole set of languages or commands. 
The term integrated means that these eete of 
commands, or languagee, are so implemented that 
a bridge designer or roadway designer, a geo- 
metriee designer or foundation designer can all 
eroee-reference a common data base; meaning 
that any one designer who is concerned with the 
location of a highway can, via commands, carry 
out a loeaUone study which sets up files which can 
be cross-referenced by another man doing, lefe 
eay, geometric design of an interchange. Hie files 
in turn can be eroes-refereneed and utilized by 
the man designing the bridge, and so on. There 
is a whole set of commands for the various 




technical dieeiplinee which are integrated in the 
eenee that one can transfer data between these 
eub-systeme and carry out in one eontinuoue pro- 
eeae, using one master set of filee, a rather com- 
plete design job, I agree with you entirely that 
this matter of the objectivity, the criteria, or what 
it ie you*re trying to do, ie still open. The first 
operational versions of ICES wiU rely almost 
entirely on the designer himself, to provide 
these criteria. He will have a very powerful de- 
sign tool; how he will use it, though, ie still going 
to be up to him, 

MORF: 7 didn*t by any means mean to eay that this 
waen*t an eventuality. But I think that there are 
a lot of eonsiderdtione that have to be thought of 
before we have something that is available in the 
real world. 



BAUMANN: We have heard various things ^out 
making languages more useful for people, X)ne 
of the things that people do very well is use lan- 
guage redundantly; eomputere do thie nb^ak alii 
The real thing Miller is talking abmi 0^tegraU 
ing some of these special progranys ha^ug 
them done under one roof. What we really used 
to aceompUsh in terms of language is to somehow 
make it possible for redundant languages to be- 
come useful. An example I think we might weU 
look to ie in the field of elqetromcs. If you take 
a new television set down to your repairman, who 
may never have seen a model of this kind before, 
he can locate right away the 'power supply, the 
high voltage section, the picture tube, etc. In 
other words, there ie a lot more pattern in the 
eleetroniee eet-up than there ie in computer pro- 
grams as we now use them. We somehow need 
to put this pattern back into our programming 
language so someone else can recognize a standard 
element by just looking at it. Maybe we can go 
completely to a system where we can program 
computers topologically rather than by a listing. 
This seems to be one of the directions which will 
be very useful in integrating some of these 
separate items; that is, we*d be able to take a look 
at a specific computer program' and eee what in- 
puts and so on are required for another routine, 

MORF: If that was a question, the answer ie yes, 

PAYNTER: 7 would hope that this issue, which has 
been raised several times, is typical of some prob- 
lems that we might come to gripe with in the 
next few days. It ie your feeling, going back Is 
LicTdider*e charge to this group, that it is just 
a matter of time before we progress to a situa- 
tion in which we can expect in our computer die- 
eouree the same type of freedom that we almost 
demand of each other if we are working together 
on a problem? Or do you eenee that thsre may be 
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«ome fundamental Umitationi to what we will ever 
he able to attain in the way of an integrated deeign 
proeeduref 

MORF: Well, I have anewered a eirnUar queetion 
in one way. Now / toUl go ahead and anewer it in 
another. This is an organizational funetUm, 
There are people who determine alignmente, there 
are othere who determine the teneion in Unee, 
or thoee who determine the voltage dietrihution, 
and 80 on, and theee different functione are realty 
in varioue eompartmente. Perhaps one of the 
restraints on a totally integrated system is the 



organization and bureauerggy of a large industrial 
organization. This again woM argue that some- 
thing must be ehanged in the way of organization 
if you expect to put the puneh card in at the front 
door and have a final design eonie out at the other 
end. This woM involve cutting across quite a 
large number of organizational Unee, None of 
these are impossible but theifre not the kind of 
things youWe going to get right away and they do 
provide restraints, and they also provide an ex^ 
planation of the reason that the solutions in the 
earlier phases are undertaken one part at a time. 
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COMPUTERS AND ENGINEERING DESIGN 

IN INDUSTRY 

D. E. Habt 

Research Laboratories* General Motors 
Corporation 
Warren* Midi. 



I was asked to summarize industry's needs for 
future engineering graduates. Quoting from the 
advance abstract* ''Industry needs engineers who 
can approach new problem situations with intelli- 
gence* imagination* and confidence and who can 
arrive at satisfactory Mi'lutions in a reasonable 
period of time.^' * 

That statement was true ten years ago. 



It will still be true ten years from now. 

* T^e abstract also states that the computer is one 
of the many tools which the engineer should be able 
hse.. I am sure you would all agree that that is 
h^i^e s^temenl^ and* if that was all there was to 
file story*" ^ would be no need for me to say 
anything . further. 






01 all the tools available to the engineer* however* 
the computer undoubtedly has the greatest potential 
for increasing his productivity as a designer. Yet 
only a small fraction of the engineers now in indus- 
try<-^d only a slightly larger fraction of the 
engineers now entering industry from the univer- 
sities— «re making use of computers in design. 



This is the problem which we all face: How do 
we put the power of the computer into the hands 
of the engineer so that he can make effective use of 
it as a design tool? And* from industrsr's stand- 
point* how do we get the computer into the main 
stream of industrial design activity? 

What do I mean by "the main stream of design?" 
In the automobile industry* it is the design activity 
which starts with a blank sheet of paper and ends 
up with finished automobiles coming off an assembly 
line three years later. 



Who works in this "main stream?" Included 
among the workers are artists and clay modders 
in Styling* but the largest number are engineering 
designers in hundreds of drafting rooms through- 
out the Corporation who spend millio ns of man- 
hours each year designing GM's principal product — 
the automobile. 



I^s pretty obvious why these designers have not 
been able to use computes— the man and the ma- 
chine don't speak each other's language; The 
language of computers is numbers* while the lan- 
guage of the designer is griiphics— lines on a 
drawing. This language pr(fi>lem will be discussed 
in more detail later. 



l^t me elaborate on tiie concept of the "main 
stream" of design. In Figure 1* the main design 
stream is depicted as an arrow pointing to the 
right 

Outside of the main stream is a function labelled 
Researdi and Development It might also be called 
advanced design* or some similar name. 

The R&D function develops new technology which 
is fed into the main stream after it has been 
thoroughly tested* p.g.* Oldsmobile's front-wheel- 
drive Toronado. I%e if unction is also available on 
a consulting and troubleshooting basis. The three 
bars in the main stream illustrate that it takes 
three years to complete the design of a new car* but 
that a new design cyde must be initiated every 
year. 

This might be refened to as a steady-state design 
environment in contrast with the space-age design 
environment where the object is to invent a way 
to do something which has never been done b^ore* 
e.g.* get to the moon and back. 

My subf^uent remarks will be primarily directed 
toward this steady-state design mivironment— whidi 
I bdieve is reasonaUy typical of much non-defense 
industry. 

Where in this environment are computers not 
being usedr~and why? 

In the R&D area there has been a rapid increase 
in computer use as an analysis tool* but it is not 
being used as much as it should be. There is a 
phenomenon here which might be referred to as the 
Threshold Problem: an engineer will use a computer 
only if it's easier than some other way. Even if 
a program exists as a black box which converts 
data sheets into printouts* the engineer must still 
fill out the data sheet* get it key punched* and t hen 
wait a few hours for the results to come back from 
the computer. As an extreme case* this would be 
a ridiculous way to solve a single quadratic equa- 
tion * however* an engmeer might need to solve 
several hundred quadratics in a month— one at a 
time. If a program doesn't exist* the engineer can 
wait in line for a programmer or try to learn (or 
recall) progra^miing in a language like FORTRAN; 
however* FORTRAN is not good for the casual user 
since there are too many details that can be for- 
gotten. 

But something exciting is happening. I like to 
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call it conversational computing. We recently in- 
stalled a GE 266 computer system with remote tele- 
type terminals which operates on a time-sharing 
basis over telephone lines using the Dartmouth- 
developed BASIC ssrstem.' BASIC is easy to learn, 
and if the engineer does forget the rules he is re- 
minded of them at once, not four hours later. Term- 
inals are presently installed in test cells, drafting 
rooms, and engineers* * offices. 

An important point is that the threshold has been 
significantly lowered — ^the computer can now be 
used quickly to solve even one quadratic. 

Most important, a whole new group of engineers 
are now using computers— primarily because they 
can remain involved in the solution of their prob- 
lems on a minute-by-minute basis. I conclude that 
it*s got to be good if an engineer will put up with 
a telelype terminal as a means for communicating 
with a computer. 

Speaking of involvement, I recently received the 
Proceedings of the Third Conference on Engineer- 
ing Design Education on ''Authentic Involvement 
in Interdisciplinary Design." * This was a report on 
a series of workshops in which faculty members at- 
temped to guide their own students through un- 
structured design projects. In this report, I was 
unable to find reference to the use of computers 
on any of the design projects. 

My conclusion is that we indeed have not suc- 
ceeded in learning how to put the digital computer 
into the hands of the engineer as a design tool. 

Referring back to the main stream of design, 
you will recall that designers now working on draft- 
ing boards are unable to use computers because they 
don't speak the same language. It has become 
apparent that the solution is not to teach the de- 
signer the language of computers, but to teach the 
computer the language of the designer. 

I would like to point out that two designers often 
engage in a "conversation" about a design problem. 
During this conversation, they make reference to 
drawings, handbooks, and models which have been 
construct^ to help visualize the problem. They 
make sketches to help express their ideas. 

Such a conversation is essential if the job is to 
be completed on schedule — ^it couldn't be done if 
they had to send letters back and forth to each 
other. I'm emphasizing this point again because 
most problems are solved on a computer just that 



ip. T. Shannon, **Impact of Time-Sharing Computers on 
Education in Engineering Design,” Conference on the Impact 
of Computers on Education in Engineering Design, April 
21-28, 1966. Sponsored by the Commission on Engineering 
Education. 

* Proceedings, Third Conference on Engineering Design 

Education, "Authentic Involvement in Interdisciplinary De- 
sign,” July 12-18, 1966: Carnegie Institute of Technology. 



way— by sending it letters. Three or four hours 
after asking the computer a question, the answer 
comes back, and if the question wasn 't asked exactly 
right, the answer may be only, "WHAT?" 

There is another aspect of the design process: 
Design can be said to be about ten percent inspira* 
tion and ninety percent perspiration. But the plan- 
ning and the doing— or in this case, the designing 
and the drafting— are so tightly interwoven that 
the designer has had to be his own draftsman— he 
has to see what one line looks like before he can 
decide how to draw the next one. In a sense, he 
engages in a sort of conversation with his drafting 
board. 

The conclusion which we come to is this : "If the 
computer is going to help a designer, he has to be 
able to use it himself," and to do this we must 
make the computer appear to the designer to be 
a kind of high-speed mechanical draftsman which 
he can control on a minute-by-minute basis as he 
now controls his drafting board. 

In this way, the computer can do inuch of the 
routine drafting, and the engineer will have more 
time to dream up the real design which only a man, 
and not a machine, can do. 

In 1960, after two years of study and experimen- 
tation, the GM Research Laboratories initiated a 
major research and development project whose pur- 
pose was to learn how to put the power of the 
digital computer directly into the hands of these 
designers, and ultimately to move the computer into 
the main stream of GM's automobile design. 

This project is known as Design Augmented by 
Computers— DAC-I.’ Let me stress Design Aug- 
mented by Computers — ^not design by computers. 
This type of work is often referred to as computer- 
aided design* a term originally used by MIT. 

A key element in the DAC J syst^ is a graphic 
console which is attached directly to the Research 
7094 computer (Figure 2). It makes possible con- 
versation between the designer and the computer 
in the language of the designer which is, of course, 
graphics. 

This graphic console is a one-of-a-kind unit which 
was designed and built to our specifications by 
IBM. 

This computer can display drawings to the de- 
signer on the cathode-ray tube. He, in turn, can 
give directions to the computer by pointing at items 
in the display with the special pencil or by means 
of the keyboard. 

Inside the computer, and available on a moment's 
notice, are stored his drawings and his drafting 

* E. L. Jades, "A Laboratory for the Study of Graphical 
Man-Machine Communication,” Proceedings, 1964 Fall Joint 
Computer Conference, The American Federation of Informa- 
tion Processing Sodc^es. 
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tools — ^the computer equivalents of the T-square, 
triangle, ruler, compass, french curves, etc. 

All of the many computer programs— the soft- 
ware — which make the DAG-I system work were 
developed by the GM Research Computer Tedi- 
nology Depajiment. 

The DAC-I system went into operation in early 
1963, two and a half years after the project was 
started. This was a unique laboratory facility 
whose purpose, as I indicated before, was to learn 
how to put the power of the computer directly into 
the hands of the designer. 

But one important element was missing — ^the 
designers who were supposed to be able to use 
DAC-I. Computer technicians can design computer 
programs, but we can't design automobiles. 

Both Fisher Body and Styling agreed to assign 
designers to the DAC-I project. These were real 
designers, not computer people, and it wasn't Icng 
before they were making significant contributions 
based on their experience and knowledge of design. 
During their first year using DAC-I, they were able 
to design a complete trunk lid, both outer panel and 
inner panel, on an experimental basis. Figure 3 
shows the results of their efforts. This drawing 
was recorded directly on film by the computer using 
the Image Processor unit of the DAC-I system. 

In a written paper, it is not possible to describe 
the dynamic nature of the conversational interaction 
between designer and computer using a graphic 
console. For those who are interested, a 3-minute 
movie on DAC-I is available on loan. * 

Figures 4 and 5 are examples of the type of image 
which is displayed for the designer at the console. 
Figure 4 is a blow-up of a small portion of the 
inner panel of the trunk lid of a car. It shows the 
cutout hole and reinforcement where the courtesy 
light will be mounted. 

When the designer studied this display, he saw 
that the subassembly was not properly located rela- 
tive to the surrounding surfaces. In only a few 
minutes at the console he was able to supply all 
of the information necessary to redesign this sub- 
assembly 1^" to the left. The results are shown in 
Figure 5. 

With experimental applications such as this, we 
have shown that it is perfectly possible to place the 
power of the computer into the hands of the engi- 
neer so that he can use it as a design tool without 
becoming a computer programmer. I should point 
out that this was no small job— the DAC-I software 
runs into well over a million computer instructions. 
But with the DAC-I laboratory we have demon- 

^‘*Design Augmented by Computen: DAC-1,’’ available 

from the General Motors Corporation, Public Relations Staff, 
Film Library, General Motors Building, Detroit, Michigan, 
48202. 



strated how to make the computer available to a 
large group of engineers who previously had not 
been able to use it. 

Perhaps there is a parallel between what I have 
just described and the field of engineering education. 

If we change only a few words, Figure 1 pretty 
well characterizes the educational process. The 
arrow now becomes the main stream of education; 
the products engineers. The university, too, 
could be described as a steady-state environment; 
it takes four years to produce an engineer, but a 
new crop of students enters each year. Most uni- 
versity computer use has also been in the R&D 
area, and so far it has had very little impact on the 
main stream of education. 

I would like to suggest that the university's goal 
should be to learn how to place the power of the 
computer directly into the hands of the student and 
how to start moving the computer into the main 
stream of engineering education. 

PAYNTER: I gather that now weWe ready to 
coim hack to MiUer^a master d^signerf because 
he is obviously not the person Hart was speaking 
off certainly not the three-, four- or five^ear 
product, although perhaps the six-, seven- or 
eight-year product. I got the distinct impresswn, 
from Hart*s picture of two men talking to each 
other over the drafting board that these people 
were not necessarily master designers, but simply 
competent, typical engineers. Would you care to 
suggest how you think schools might make this 
particular group of people aware of what they 
could do with computer sf 

HART: I think this comes back to the last topic, 
which w(^ that part of the goal of the educational 
process is to impart to the student new concepts 
and ideas, and the goal of improved education, 
by rnoving the computer into the mainstream of 
engineering education, is to increase the rats of 
concepts per unit time that you*re ahle to stick 
into the headbone of the student. My initial re- 
marks were that we are interested in engincors 
who have intelligence, imagination and the ability 
to approach a new job with confidence (so far I 
havsn*t seen that the university imparts anything 
to the student but the intelligence and the rest 
of it he learns as he goes along). Part of the 
current attempt to get studonts involved in actual 
^sign problems during their university training 
is the wish to give them something more than 
tools to work with, and also some knowledge of 
how to apply these tools and what it means to 
apply these tools. I think this is very worthwhile. 

I would like to see the computer used more in this 
type of design workshop. 
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LOWE: It appears to me that today we have 
sensed, in the minds of many people here, the 
feeling that there is a real harrier between the 
computer designer and the computer user. How- 
ever, you said that the men who design the car 
body actually contributed to the design of the 
system, and that they came in as guinea pigs and 
in a very short time turned into contributors. I 
wonder if you could amplify this a little. 

HART: We could have spent ten years doing a 
systems study on how automobiles get designed. 
But we concluded that we weren*t going to be- 
come smart enough fast enough to learn hou) to 
design automobiles so that we could design an 
automobile design system. So our goal and ap- 
proach had to be entirely different. We observed 
one thing about the design process: whatever a 
man was designing, whether it was automobiles, 
spark plugs, or what have you, he used the same 
tools. He used a ruler, a french curve, a compass, 
and a straight edge. So our basic approach was 
to provide the designer with a set of reasonably 
primitive tools and means by which he could 
operate these tools in sequence in order to carry 
out a design process. As I indicated, we had to 
develop these tools in a kind of a vacuum, be- 
cause, when we approached the designers and 
asked them what sort of things they would need 
if they were going to use a computer to help them 
with design, they just raised their oyebrows at us 
and suggested that we quietly retire to our lab- 
oratory and leave them alone to do their jobs. So 
we had to da the best we could to buUd what we 
thought would be an adequate set of tools to get 
started. And the the designers came in as 
guinea pigs and attempted to use these tools to 
solve real design problems on an experimental 
basis. They rapidly discovered that some of the 
tools we provided weren*t good and some of the 
tools they needed we hadn*t provided. But an- 
other phenomenon also took place: as they com- 



menced to use the tools, they began to get some 
insight into the real capability which this new 
machine would provide for them. For example, 
if you work all day on a sheet of paper which is 
two-dimensional, you begin to think in terms of 
a two-dimensUmal world. We found that the de- 
signers who came in and worked with us did just 
that. They were working in a two-dimensional 
environment and it took them a while to adjust 
to the fact that, in working with the computer, 
they were really dealing with a three-dimensional 
medium, and they had to think in terms of work- 
ing in three dimension. So their concepts of how 
they ought to do design changed as the tools 
changed. These two phenomena interplayed back 
and forth in developing the laboratory to the 
state where it is today, 

QUESTION: Would not the same system, which 
could be the operational base upon which your 
engineers actually conduct design tvithin General 
Motors, be a very powerful tool for teaching a 
nsw generation of designers within the engineer- 
ing school? Suppose that we could afford to give 
the young men, before they came to you, two 
years of apprenticeship training under the guid- 
ance of some pretty good design people. Wouldn*t 
the DAC system be a wonderful way to teach de- 
sign in a university? 

HART: I think it would because ifs a rather pain- 
less way to learn descriptive geometry and the 
visualization and interactions of things in space. 

BRANIN : / just want to make a brief summary of 
a couple of points which seem to have emerged 
from our session today. The question of acces- 
sibility to the computer has been brought out. It*s 
really not accessibility to the computer per se, 
but to its facility. The other question was: how 
do we get engineers to use the computer? / 
think these points both have been very sharply 
highlighted by Hartfe discussion. 
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COMPUTERS IN ENGINEERING DESIGN 

Donald L. Katz 
The Univerai^ of Biichiffan 
Ann Arbor, Mich. 



At a meeting in 1962 of the Committee on Engi- 
neering Design of the Commission on Engineering 
Education, one task outlined for engineering edu- 
cators was the development of the computer as a 
design tool. Following these discussions, a pro- 
posal was written to (^tain support for the project 
which is being reviewed on this program. In the 
spring of 1964, a grant was received from the 
National Science Foundation for support. 

ENGINEERING DESIGN 

In the span of years encompassed by my contact 
with engineering education, engineering design has 
gone through three stages. The first stage was that 
of adopting the current technology along with modi- 
fications required by the new process under con- 
sideration. In 1980, as a student, I designed a 
petroleum separation plant utilizing information 
available on fired heaters, distillation columns, 
condensers, and treaters. Gradually there appeared 
a oonfiict be^een the devdoping engineering sci- 
ences and the art in form of technology. Accord- 
ingly, design problems of the second variety were 
developed. These problems required the creation 
of a process or system based on fundamental cal- 
culations for the various components using rigorous 
rate process equations, separation process rela- 
tiondiips, and othsr quantitative material. Such 
designs became rather involved and tedious. It was 
necessary to make a good choice of the parameters 
in the first place to permit completion of the design 
problem within the time limits permitted in the 
educational environment. 

With the advent of computer usage in the engi- 
neering sciences, it became rather clear that engi- 
neering design would soon be relegated to the 
computer — ^in good part because of the need for 
rapid calculations. It appeared that the universities 
should find a format by which design could be 
upgraded to include the computer techniques being 
adopted by industry. Also, the need for having 
design teachers conducting studies at the forefront 
of knowledge in this area was indicated. It was 
this view and the discussion with the Commission 
on Engineering Education’s Committee on Design 
that prompted the preparation of the proposal for 
oui; project 

The project organized in the fall of 1964 at The 
University of Michigan was modeled after an earlier 



iPir Jeet on Use of Computers in Engineering Education, 
lUnirersitj of Michigan (1969-68). 



Ford Foundation project ^ It had as major objec- 
tives: (1) the training of engineering design teach- 
ers in subjects related to computer-aided design, 
(2) the study of the role of the computer in engi- 
ne^ng design education, and (8) the generation 
of a substantial number of completely documented 
computer-oriented engineering design problems. 

In the fall of 1964, engineering design teachers 
with substantial computing experience were a^iked 
to apply for an intensive Mxuvvvedc program to be 
held at the University of Mkhigan during the mmi- 
mer of 1965. An Advisory Committee of the Com- 
mission on Engineering Education* assisted in 
selecting from among the applicants twenty-nine 
engineering design teachers from twenty-three 
different engineering schools (list attadied). 

Five weeks of the summer program were devoted 
to formal presentations on computation, numerical 
mathematics, modeling and simulation, mathemati- 
cal optimization techniques, problem-oriented lan- 
guages, man-machine interaction^, economics, and 
industrial design practiee. During the remaining 
four wedcs, each participrtmg professor solved and 
documented one or more individual design problems 
appropriate for inclusion in an engineering design 
course in his own field. 



During the fall of 1966, the project staff and one 
summer-program participant from^ each of five 
engineering disciplines edited or revised solutions 
to the participants’ design problems and prepared 
material for the final project report. This report 
consists of six volumes, one summarizing project 
activities and five oriented toward design in indivi- 
dual engineering disciplines as follows: 



Volume I. 
Volume II. 
Volume III. 
Volume IV. 
Volume V. 
Volume VI. 



Summary 

Chemical Engineering 
: Civil Engineering 
Electrical Engineering 
Industrial Engineering 
Mechanical Engineering 



Volume I includes an historical description of the 
project, details of the summer program, recom- 
mendations, review papers on the subjects covered 
during the summer program, and abstracts of 
forty-three of the design problems prepared by the 



>N. M. Newmark, Chairman (Univeraity of Michigan) 
S. E. Elmaghraby (Tale University) 

S. J. Fenves (Univerdty of Illinois) 

D. E. Hart (General Motors Corporation) 

A. Sehnlts (Comdl University) 
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summer program participants. Each of the five 
disciplinary reports includes a brief statement on 
the role of computer-aided design in the field and 
several completely documented computer problems. 
These reports result from the combined efforts of 
the project staff and the first author as principal 
editor. The example problems in some cases were 
modified after correspondence between editor and 
problem author. 

Some general conclusions and recommendations 
of the project staff are: 

1. Engineerin^r sdiools should continue to develop 
introductory computing courses to provide 
every engineering student with an adequate 
background in digital computation and instruc- 
tion in at least one procedure-oriented pro- 
gramming language. 

2. Computing work should be incorporated into 
the required upper-level analysis (engineering 
science) courses; before graduation, each 
student should solve msny engineering prob- 
lem of graded difficulty using a computer. 

8. If 'Stilts introductory programming course is 
taught early in the curriculum, as usually oc- 
curs and is recommended, then there should 
be a later treatment of more sophisticated 
mathematical topics available, to engineering 
undergraduates. The student will need: 

a. More advanced mathematics, in particular 
numerical mathematics and mathematical 
statistics; 

b. more careful instruction in the techniques 
of mathematical model b uilding and simu- 
lation; 

c. an understanding of the most important 
of the mathematical optimization tech- 
niques; 

d. exposure to a variety of computing sub- 
jects, including special purpose computing 
languages and the manipulation of graphi- 
cal information. 

4. The mathematical tools (numerical mathe- 
matics, mathematical model building, optimiza- 
tion techniques, etc.) needed for successful 
integration of computer work into engineering 
design courses are basically the same for all 
engineering disciplines, although there are 
some differences in emphasis among the en- 
gineering specialties. Therefore, it is not 
unreasonable to teach such material in an 
interdisciplinary course. At the same time, it 
should be recognized that problems are most 
meaningful when taken from the ai^ of 
interest of the student. Hence, problems used 
in an interdisciplinary course would have to 
be chosen very carefully. 

6. If we expect students to use the computer to 



its full potential for engineering design work 
and .wish to assign rather sophisticated prob- 
lems, special purpose problem-oriented lan- 
guages should be used when available. If this 
is not practicable, then at the very minimum 
a substantial library of useful design sub- 
routines should be available. 

6. Design problems usually have an infinite 
number of feasible solutions, and experience 
and intuition play a significant role in the 
solution process. Engineering design proce- 
dures are often very difScult to automate 
completely. In order to use the computer most 
creatively, the designer must be able to con- 
verse with the computer in a virtually conver- 
sational way. Turnaround times of days, hours, 
or even minutes are not generally acceptable. 
The interaction must be effectively instantane- 
ous. 

7. Time-shared computers will have a great 
impact on the place of computers iu engineer- 
ing education in general. Thei student will 
be able to debug his program and solye a 
proUem at one sitting. He should be able to 
solve reasonably difficult problems between 
class meetings. In addition, it will be practical 
for the engineering teacher to bring an on-line 
console, perhaps with graphical input and 
output, directly into the classroom. He will 
be able to display solutions for all to see as 
the computation is being done. Engineering 
faculty should be prepared to teach the use 
of such approaches within the next two or 
three years. 

8. With ^e introduction of time-shared hardware 
and the development of problem- or design- 
oriented software, the student designer for the 
first time will have the opportunity to investi- 
gate a large number of feasible solutions to 
his design problem. He will be able to gain 
some facets of engineering experience which 
cannot now be acquired outside the industrial 
environment. Hence, the next generation of 
computing systems should enable the student 
to improve his intuition about and under- 
standing of the nature of good design, which 
is, after all, the raison d'etre of engineering 
design courses. 

9. For financial reasons, it may not be possible 
or practical for a small school to acquire a 
large time-shared computing systmn for its 
exclusive use. The best approach for many 
schools will be to acquire remote terminals 
connected to a central processor shared with 
other university or industrial users. Even 
the larger schools will find the total cost of 
operation, though smaller on a unit computa- 
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tion basis, to be substantially larger than for 
the. present generation of. computing equip- 



ment. A very large increase in demands for 
computing services is almost certain to take 
place over the next few years. This will 
probably mean commitment of a growing 
share of the university’s budget to support the 
computer and its services. University adminis- 
trators should be planning now for support 
of these activities. 

10. The key to successful integration of computers 
into engineering design course work is the 
knowledgeable and interested teacher. A 
considerable educational effort will be required 
to train our engineering design teachers to 
use the powerful new computing systems and 
languages. 

Copies of the summary and individual disciplinary 
reports are being distributed to appropriate depart- 
ments at all accredited engineering schools. Hard- 
bound copies of the six-volume report are being 
sent to engineering libraries. Copies of individual 
report volumes will be furnished without charge to 
full time engineering faculty members. A limited 
number are available for sale to industrial organi- 
zations from the Publications Distribution Service 
of The Universi^ of Midiigan. 

During the 1966-67 school year, the project staff 
expects to conduct several two-day “information” 
sessions to stimulate interest in introducing com- 
puter work into engineering design courses. Rep- 
resentative engineering design teachers will be 
invited to attend local meetings to be held in New 
York, Boston, Atlanta, Dallas, Denver, San Fran- 
cisco, and Ann Arbor. In addition, an intensive 
two-week workshop will be conducted in Ann Arbor 
for about fifty engineering design teachers. 

Staff and Lecturers on Project 
Brice Carnahan, Associate Director in Charge 
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Dale Rudd, Universify of Wisconsin 
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Paul Shannon, Dartmouth College 
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PREPARING THE UNDERGRADUATE 
FOR COMPUTER-AIDED DESIGN 

Charles R. Mischke 
Iowa State University 
Ames, Iowa 



The principal effort in the preparation of the 
undergraduate for participation in computer-aided 
design is made before his last design course. The 
preparation by the student is not as involved, nor 
as extensive, as superficial mcamination might indi- 
cate. However, the woric to be done by the faculty 
seddng to offer computer-aided design is different 
to an extent that is greater than is generally ap- 
preciated. The purpose of this paper is to state the 
problem and to offer a technique for solution. 

DEFINITIONS 

An engineer plans and/or controls (i. e., designs) 
the interaction between matter, energy, men and 
money in order to accomplish a specified purpose. 
Design is the essence of the engineering function, 
although there are many engineering functions that 
are not design. 

To design is to specify precisely how a task is to 
be accomplished. To specify requires creativeness 
and understanding. ^ Clearly, design is not taught 
in school. It is learned, if, indeed, it is learned at all, 
in the fidd, with and under a master artisan. 
Schools can, however, teach attitudes, approaches, 
conceptual tools and an appreciation of the mor- 
phology of design. 

The basis of understanding is usually analysis, 
which is itself based upon mathematical and em- 
pirical models and experience, and is a useful means 
of evaluating alternatives. Selection of the best 
alternative requires: 

(1) ability to recognize a satisfactory solution; 

(2) abil% to assess the merit of alternatives and 
order them accordingly. 

Synthesis means the generation of a solution and 
it is arrived at by understanding the problem. The 
degree of merit of the synthesis varies directly 
with the ability of the student to understand. That 
is, poor understanding generates only weak solu- 
tions, whereas complete understanding may enable 
the student to derive a solution of high merit at 
the first attempt. 

Engineering advances in the field of analysis pro- 
mote understanding (and therefore the degance 



^Strangely enough, although the handmaidena of design 
are ereativeness and understanding, engineering schools do 
not select students on creative potential, nor do they as a 
rule make earnest constructive efforts to poliah or cull on 
creative bases. 



of the synthesis) until the ultinuite form of ssmthe- 
sis is possible. Although this state of affairs is a 
goal of engineering, its very attainment relegates 
such synthesis to the status of routine. 

The frontiers of engineering include the excite- 
ment (and demands) of design work in a domain 
wherein understanding is poor, but the need is great. 
Hence there is a high challenge in searching for 
meritorious solutions by better analysis techniques 
(one route) or by clever search techniques in the 
domain of satisfactory alternatives (another route). 

PERSPECnVES CONCERNING THE COBIPUTER 

The role of the computer is becoming more power- 
ful as the ability to simulate systems on it is devel- 
oped. The computer calculates with speed, precision 
and without fatigue. Are these abilities useful to 
the designer? Certainly. 

Should ability to exploit these capabilities be 
taught to undergraduates in engineering? The 
rapid growth of computer utilization in industry 
over the last ten years indicates **perhaps,** but the 
advent of time-sharing computer hardware can lead 
to the development of a computer utility, like an 
electrical power utility, which could place in every 
engineering office a remote access to the largest 
computers available. The use of large computers will 
become a common part of design life. The advent 
of time-sharing will demand that four-year engi- 
neering graduates have some experience working 
fruitfully with a computer and have a broad under- 
standing of current useful techniques. 

Let us consider wherein computer-aided design 
assists engineering education, and what skills are 
needed. The engineering designer often proceeds 
as follows: 

(1) Defines his problem in engineering terms. 

(2) Decides how he will recognize a satisfactory 
solution. 

(3) Determines how he will recognize a better 
(or best) alternative. 

(4) Generate alternatives. 

(6) Through analysis, retains suitable alterna- 
tives. 

(6) Chooses best alternative. 

(7) Ad» upon the decision in (6). 

Since the computer can oidy converse in a mathe- 
matical language^ the hovering of a computer in the 
background affects these seven steps as follows: 
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(1) The engineer must eventually prepare a 
mathematical statemmit of his problem for 
the computer. The pressure to state his 
problem in mathematical engineering terms is 
a wholesome influence. 

(2) A description of a satisfactory solution in 
mathematical terms is required by the com- 
puter; every constraint upon solutions 
should be involved here. The computer re- 
quirements pressure the designer to be 
orderly, explicit, free of ambigui^, quantita- 
tive and exact. 

(3) A figure of merit, (utility function, object 
function, etc.) is always required by engi- 
eers. It may be non-quantifiabl^ such as 
an experienced expert's value judgement (a 
tea-taster, for example) . The computer can- 
not function as a tea-t^ter unless the infor- 
mation and experience is placed in mathe* 
matical form. Again the pressure is upon 
the designer to make his utility function 
mathematical. (It also prevents double-talk, 
such as **I made it as cheaply and reliably 
as possible," while offering a mediocre alter- 
native.) A computer program would reveal 
explicitly the designer’s utility function. 
Again, the pressure to quantify judgements 
is wholesome. 

(4) In this generation of alternatives the com- 
puter works with a speed and accuraqr with 
which no man can compete. In this step the 
economic contribution of the computer is 
strong. 

(5) Analysis programming has a history as long 
as that of the computer, and again the speed 
and accuracy of the computer exceeds that of 
manual calculation. 

(6) Using the computer for this step removes 
much of the subjective prejudice in deci- 
sion-making ("I like a cam for this applica- 
tion rather than a four-bar linkage.**). Dis- 
enchantment with computer decisions tells 
the engineer either that his programming of 
step (8) is open to improvement, or that 
his prejudices are showing. 

(7) The computer even contributes to the imple- 
mentation of the decision. If the engineer 
needs clearance from a senior before he can 
proceed, Ibe coniplete computer program to- 
gether with the recommendation will allow 
the senior to quickly inspect the merit func- 
tion employed and consider the implications. 
The latitude available to the decision-maker 
is apparmit. A very flat figure of nierit sur- 
face in the domain of the extremum allows 
considerable variation in parameters at 

cost to accommodate non-quantifiable con- 



straints, or whims, or prejudices of the deci- 
sion-maker. 

If we concede that the computer can be helpful 
to the designer (and therefore to the student de- 
signer), what is the price of giving a suitable 
introduction to computer utilization in design to the 
four-year engineering student? 

First, let us be practical. If an engineering prob- 
lem had to be solved and programmed from scratch 
on each and every occasion, the computer would find 
little use in the design room. The overhead to plan 
and debug the computer program would be large, 
and only a small number of programs would be- 
come economically feasible. Any notion tha^ a prob- 
lem once successfully completed would be encoun- 
tered again should be suppressed. 

If the overhead of time and money associated 
with computer use is substantial, how can previously 
invested time and experience ^ salvaged so that 
the long-term cost per problem becomes reasonable? 
^e answer to this question is the same, surpris- 
ingly, whether you are talking about a human 
problem-solver or a mechanical one. A human is 
trained to provide conditioned responses to very 
elementary situations. A designer functions be- 
cause he is able to bring to bear upon a problem 
a whole bagful of conditioned responses, arranged 
and ordered by him through intellect, in such a 
way as to solve the complicated problem at hand. 
Just as the human is taught a repertoire of condi- 
tioned responses, so can the computer be instructed. 
It then becomes the role of the human to arrange 
the responses into a sequential strategy for a solu- 
tion. The computer capability corresponding to the 
learning and remembering of conditioned responses 
to elmentary problems is called the subroutine 
capability of the computer. 

A design room (or an educational department) 
requires a library of subroutines, available for 
quick access in any specified order by the designer 
(student) wishing to use the computer for design. 
The details on how to assemble and how to document 
such a subroutine library are not germane to my 
subject, but the existence of such a library is pre- 
sumed for the remainder of this paper. 

Given a problem in analysis, all that remains for 
the computer user to do is to plan a strategy of 
calling library subroutines in a specified order to 
initiate computation which will produce the required 
analysis in numerical terms. Our subject here is 
design by computer. The previous procedure is 
necessary to computer-aided design, but it is not 
sufilcient. 

What can a computer do? It can add, subtract, 
multiply, divide and perform other concerted actions 
which we simply call calculation. The fundamental 
and critical activity in design is decision-making. 
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Gan the computer make decisions 7 The computer 
can take direction and make logical decisions of the 
following kinds: 

a. It can take direction to the extent of trans- 
ferring control to a given location in the logic 
flow diagram of a program, regardless of the 
circumstances (the unconditional GO TO state- 
ment). 

b. It can make a quantitative comparison and, 
depending upon whether a mathematical ex- 
pression is positive, negative, or zero, it will 
branch in its logic and undertake a different 
course of action for each possibility (the IF 
statement). 

c. It can make a quantitative determination and 
undertake a number of different courses of 
action depending upon the value of a mathe- 
matical quantity (the computed GO TO state- 
ment). 

d. It can perform a calculation or a routine of 
calculation iteratively for any specified num- 
ber of times (the DO statement). 

e. It ca n read and write alphanumeric informa- 
tion. 

f. It can make decisions based on the truth or 
falsity of a statement and proceed on a logic 
path dependent upon that determination (the 
Boolean algebraic capability). 

g. It will remember (store in memory and recall) 
quantitative and alphabetic information, and 
it will seem to remember how to do a specified 
task in which it has had prior experience (the 
SUBROUTINE capability). 

h. It can draw, i. e., plot, two-dimensional figures, 
and (with some computers) draw on a cathode- 
ray tube to communicate with the designer 
(the graphical output capability) . 

We can that the computer has some decision- 
making capabilities. To the extent that the designer 
can reduce his decisions to combinations of the 
computer’s decision capabilities, then the computer 
can be trained to function as a capable apprentice. 
In some things, such as in speed and consistency, it 
will excel its master. 

The analysis capability of the computer with a 
library of subroutines serves the designer in step 
(6) of the design process. Step (2) is usually 
suggested by the constraints unearthed in step (1). 
Step (8), which is central to design, is often sup- 
pressed, or provided by a judgement. The computer 
demands that this step be quantitative, l^e de- 
signer must provide a quantitative relationship 
whose magnitude represents the merit of a satis- 
factory alternative. As soon as this is done, the 
problem remaining is to search over the domain of 
the satisfactory alternatives, and seek the one with 
the highest merit (largest value of the merit ftmc- 



tion). This is a procedure that can be effected with 
Uie computer and involves its decision- m a king 
capabilities. Step (4) requires creativeness and the 
computer (a model of consistency) has no creative 
talent. 

A TECHNIQUE FOR COMPUTER-AIDED DESIGN 

We have mentioned the elements of computer- 
aided design. A technique for making it practical 
consists of the following steps: 

a. The designer-student writes an EXECUTIVE 
PROGRAM to read in all necessary informa- 
tion and write the output documentation to his 
problenL This executive program calls the 
SEARCH subroutine. 

b. The SEARCH subroutine is a library program 
that will crawl over a merit hyper-surface and 
discover the location of the largest ordinate 
to the accuracy specified in the executive pro- 
gram. The search program calls the MERIT 
subroutine. 

c. The designer-student writes the MERIT sub- 
routine which generates alternatives and, if 
they are satisfactory, calculates their merit 
In so doing this routine calls DESIGN library 
subroutines. 

d. The DESIGN library subroutines are analysis 
routines that have been written, checked out 
found useful and are documented in a con- 
venient way to facilitate their use. (These are 
largly prewritten routines occasionally aug- 
ment by the designer-student.) 

In the beginning of an operation, as in a school 
seeking to initiate computer-aided design ^perience 
among its students, the design library will be modest 
and prewritten by the design faculty, but, as prob- 
lems are solved and new routines add^ to the 
library, an impressive design capability can be 
amassed. 

Now, given the design library and the search 
program, just what mathematical preparation must 
the students have in order to use the computer for 
design? 

The composition of merit functions is technically 
simple, but conceptually can be very complicated. 
This is not a computer subject; it is a design subject 
When lacking much sophistication, one can resort 
to the use of the simpler merit functions such as 
weight, cost, or geometrical attributes. 

Searching procedures need not be learned to any 
great extent because the search program is provided. 

The designer’s programming problems will be 
principally those of solving algebraic, transcenden- 
tal, differential and simultaneous equations, evalu- 
ating integrals and derivatives, curve fitting, and 
using some probabilistic ideas. 
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To a largo extent most of the places where these 
techniques are lu^ and applicable are within the 
design subroutines themselves, which are provided 
principally by tibe design subroutine library. To 
Mme extent the methods are used by the designer 
in constructing his merit subroutine. 

The University of Midiigan Project produced 43 
problems, solved and documented by individuals 
working without the benefit of a design subroutine 
library or a rather general search program. Tweniy- 
nine of the problems utilized optimization (search) 
techniques and were of varying complexity, cutting 
across most of the engineering disciplines. Never- 
theless, the list of optimization techniques and 
numerical methods is neither large in extent nor 
impressive in complexity. 

Thus, an instructor can choose several design 
problems constituting the entire computer-aided 
design experience wherein the level of sophistica- 
tion and involvement with specialized techniques is 
in keeping with a very modest preparation and 
ability among his students. Not much numerical 
analysis preparation is eiisential to show the student 
the potential of computer-aided design. 

Although some fine problems can be solved with 
a modest executive routine, how can some prepara- 
tion in mathematical methods suitable to computer 
impl^nentation be given to the undergraduate 
student? 



A SUBTLE PREPARATION 

Wher^ oh, where, in an already overcrowded 
curriculum, can this material be introduced without 
displacing other fundamental studies? 

This is a perplexing question. Perhaps it is time 
for the engineer to re-examine his view of mathe- 
nmtics and its use in the light of the computer’s 
role, present and future. The place to teach the 
materi^ needed for solving equations, integrals and 
derivatives is in the mathematical sequences re- 
quired of all engineering students and, to a lesser 
extent, in every engineering course in the curricu- 
lum. 

Let’s be specific. Is the time approaching when 
it would be more useful to view the entire spectrum 
of mathematics germane in engineering from a 
numerical analysis viewpoint? Cannot differentia- 
tion be taught from a finite difference and numerical 
analysis viewpoint? Cannot infinite series be 
riewed as partial sums with compensating trunca- 
tion errors? Cannot int^;ration be viewed as a 
summing up of rectangular (or trapezoidal) areas? 
Can’t sophomores study linear algebra, matrix 
ideas, and solution to simultaneous equations (alge- 
braic, differential, partial) together with compute 
implementation of same ? Under these circumstances 
students would not know a Bessel f un ction when 
they saw one, but do they know one now? 




In other words, cannot tibe basis for understiuid- 
ing tins m a tiiem atical matmial be brought into 
coincidence with tiie basis of computer techniques? 
To do so would place a double-edged sword in the 
hand of the student. If the mathematical sequence 
of some twenty-odd smnester credit hours taken by 
engineers is taught to a lesser extent from the tra- 
ditional viewpoint and to a greater extent from 
a numerical analysis vantage point, the student 
pattern of thought would be much more computer 
oriented, even if the computer were never men- 
tioned! (And his mathematics would be no less 
useful.) 

Given a new mathematical sequence and a simple 
programming experience in a language such as 
FORTRAN in the first year, and given continued 
interaction with the computer in the 
cowses, more experience with it thereafter is re- 
quired.^ In order to ke^ alive the student’s pro- 
gramming ability and his numerical analjrsis ideas, 
every engineering course should require at least 
one propamming problem each semester. I am not 
suggesting the use of just any old program to keep 
a curriculum committee happy. The program 
problem should be selected in order to the 
course better for its inclusion. Such opportunities 
abound. 

Under this arrangement, engineering students 
readiing their senior year and first comprehensive 
design course (not necessarily their first design 
course. There is something to be said for a design 
course ev^ year) will be able to use the computer 
as a genuine tool in their design problems. 

Design is a broad word and the meaning of the 
words design problem here is one in which the 
totality of the student’s prior experience is drawn 
upon, one that taxes their creative skills, and one 
that requires a search for optimality (say cost) or 

a trade-off optimum (say weight, cost and relia- 
bility). 



A FRINGE BENEFIT 

Probably the most important reason for intro- 
ducing computer programming early in the 
curriculum has to do with the handmaidens of de- 
sign— creativeness and understanding. Engineering 
students are attracted to enginee ring partly as a 
place to vent their creative abilities. The usual 
engineering curriculum stifies this natural enthu- 
siasm by demanding “conform now— create later” 
until by the senior year they are largly sheep. (If 
my colleagues doubt this, try the following experi- 
ment: Ask a group of seniors and a group of 
freshmen for an “off the top of the head” solution 
to a broad problem. You’ll be surprised where the 
creative answers originate.) 

The writing of a computer program is a creative 
act. Since there are many, many ways to skin a 
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cat, creativeness can flourish here^ and there is little 
cost associated with it. Mistakes are not pri^^ 
in cast iron, intricate machining or non-functioning 
equipment. The merits of alternative approadies 
can be demonstrated quickly. 

Another act of the designer which is submerged 
in engineering education is that of decision. Cour- 
age in decision-making comes from confldence in 
tools of decision and a prior history of success. A 
computer program is replete with decisions and the 
consequences follow immediately (within the turn- 
around time of the computer installation) . What 
better opportunity to sharpen the tools of decision 
and build confld^ce than in this aspect of computer 
programming? Thus, a simple course in program- 
ming can have secondary products of exercising cre- 
ative and decision-making skills. 

The changes indicated here will only come about 
through the initiative of the engineering educational 
community, courage to experiment; hardheaded 
evaluation, and considerable effort whose only sure 
reward is the excitemrat of doing and contributing. 
The dangers inherent in such changes are the same 
as in any change. Perhaps the most important 
necessary (but not sufficient) item is whole-hearted 
faculty support of whatever plan is implmnented. 
Nothing dulls a student’s enthusiasm more than 
to have had one professor assure him something is 
important^ only to have a subsequent professor 
show disinterest (or ignorance). 

SUMMARY 

The computer can perform the design function. 
An undergraduate student can be given a meaning- 
ful introduction to the computer’s capability in 
design if : 

1. He receives a programming course in his fk’st 
year. 

2. His mathematics sequence is enriched with 
ideas and approaches from a numerical ana- 
lysis viewpoint, and the computer is used as 
a teaching aid. 

3. Every engineering course requires one prob- 
lem in which the computer is used. 

4. The faculty develops a design subroutine 
library in advance. 

5. The faculty develops an efficient search sub- 
routine. 

6. Consideration is given to formulating figures 
of merit. 

7. A design experience is permitted in which the 
student writes executive programs and merit 
subroutines to sdve design problems. 

Once convinced of the potential of the computer 
as a design tool, the engineering graduate will be 
motivated to learn on his own all he can concerning 
methods and approaches useful in computer work. 






QUESTION: You keep talking about developiny a 
- ewnputw and then keeping 

it alive until you can use it. Why not develop 
it just before you need it'i 

MISCHKE: Ideally 1 1 think, a man ought to be ex- 
posed to design from the freshman year all the 
wa/y through. There*s a lot to be said for a 
design course every year in the man*s experience, 
and it can start out very trivially in terms of 
background required and give some ideas con- 
cerning optimization and simple problems (they 
might even be ones to show up the math depart- 
ment) but problems which are too hard to do 
by pushing a pencil. 

QUESTION : Are you then equating the use of the 
computer and the use of the design education ae 
being co-equals? 

MISCHKE: No. The computer is a tool. I think 
I made it clear that cooking up a merit function 
and figuring out how to generate alternatives are 
design activities. The machine does the rest. 

QUESTION: 7 wasn*t disagreeing with you, but I 
was trying to see whether you were saying that 
the purpose of introducing the computer into the 
engineering curriculum at all was for its use in 
design courses. 

MISCHKE: No, other people ought to use it, too. 
At my school the amount of use by the disciplines 
within M. E. is exceedingly small at the m/omcnt. 
We hope to jar them a little. 

CRANDALL: 7 know in our particurlar depart- 
ment in which people take four or, five courses, 
we give them five hours of work arid three hours 
worth of credit. Do you realize what a burdon 
this might put on a student as far as an academic 
load if every one of these classes required a com- 
puter program? 

MISCHKE: One program per course per semester? 
It doesn*t have to be very sophisticated to teach 
the ideas. In fact, that* s the challenge. 

(laughter) 

CRANDALL: No, Fm thinking that in our basic 
computer course we require approximately five or 
six problems in the one-quarter system. If you do 
this in every quarter from then on, they are in ef- 
fect taking an additional course load as I see it. 

MISCHKE: No comment. 

CARNAHAN: 7 think the new computing systems 
are going to make quite . a difference in what can 
be done in the way of assignments. That is, in 
the past it has taken a minimum of two weeks to 
do a problem with any complexity at all, bul with 
new systems I think it will certainly not be un- 
reasonable to ask a student to do a problem be- 
tween Monday and Wednesday if he can sit down 
at the console and get six runs and debug his pro- 
gram. 
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THE SIMULATION SPECTRUM 

Norbert Hauser 
Polytechnic Institute of Brooklyn 
Brooklyn, N. Y. 



This presentation is based on the Simulation 
Chapter of Reference (1), which has been distri- 
buted to Conference attendees. The following 
remarks confine themselves to a few aspects of 
model building and digital computer simulation. 

A process or system model is of practical value 
only if (a) it adequately represents the process, and 
(b) it is amenable to some form of mathematical or 
logical analysis leading to a better undei standing of 
the process. 

Currently available methods to implement (b) 
frequently require violations of (a). When there 
is a conflict between realism desired and power and 
elegance of approach, the choice appears to be 
between seeking (1) the exact solution to an ap- 
proximate problem, or (2) an approximate solution 
to the exact problem. The first choice (1) is per- 
fectly acceptable if the resulting distortion does 
not significantly affect the problem’s characteristics, 
or if the engineer is able to evaluate and counteract 
the distortion’s consequences. The second one (2) 
on the other hand implies the ability to duplicate 
an entire system (or to operate on it directly). 
Even if this were feasible, the time and expense 
required to design and execute controlled experi- 
ments necessary to isolate relevant variables makes 
this approach a last resort. In most cases the 
engineer looks for a point somewhere between these 
extremes. He is quite satisfied with a formulation 
which has abstracted the system’s essential ^ charac- 
teristics, permitting him to apply other than brute 
force methods. 

Simulation is one of the available methods. It 
requires no assumptions of linearity, differenti- 
ability, or any other form of mathematically good 
behavior. A set of mathematical and/or logical 
statements describing the exact (in a deterministic 
or stochastic sense) relationships among all compo- 
nents is needed. The simulation then produces the 
overall effect of these iteractions. It is important to 
realize that simulation need not be considered solely 
an alternative to analytical methods. Many times 
both approaches may be used concurrently. Indeed, 
one of the cardinal rules of efficient simulation 
design is not to simulate those parts of a model 
which can be treated analytically. 



^Essential with respect to a particular application. Thus 
a model designed to answer one set of questions about a 
system may be worthless when confronted with another set 
about the same system. 



MODEL VALIDATION 

Since it is impossible to prove the validity of a 
given model, the best that can be hoped for is that 
it will withstand a series of tests, performing 
reasonably and producing reasonable results. In 
analytic models the danger arises from oversimpli- 
fication introduced in order to make the model suit- 
able to mathematical analysis. In simulation, 
however, the danger frequently is of the opposite 
kind. The desire for realism may cause some im- 
port' nt features to be overwhelmed by a mass of 
realistic but trivial information. Assuming the 
primary purpose to be a gain of insight into a 
system’s cause-and-effect relationships, it is essen- 
tial that all material irrelevant to the particular 
investigation be left out. This introduces some 
circularity : if we knew exactly which variables are 
relevant, we probably would not need to simulate. 
Model construction thus becomes an iterative pro- 
cedure consisting alternately of (a) selecting varia- 
bles and specifying their effect on each other, and 
(b) running the model and observing its behavior. 
If this behavior is unsatisfactory, either some rele- 
vant variables were excluded, or the interactions 
were inadequately specified. When the behavior ap- 
pears to be reasonable, the sensitivity of all variables 
and parameters should be explored. Besides leading 
to the possible elimination of some variables, this 
will indicate the precision required of parameter 
estimates and other data used. 

Not until we have reached this point are we 
ready to collect data from the system, and to deter- 
mine how much time and effort should be expended 
in improving the accuracy and precision of this 
data. The mistake of building the model around 
available data, or of collecting the data before 
constructing the model is frequently made. The 
previous discussion should make the improvidence 
of such an approach apparent. 

EXPERIMENTAL DESIGN 

Since simulation as such provides no algorithm 
leading to improvements in succesive runs, effective 
search methods (machine or man-machine) have to 
be incorporated. This topic, including an extensive 
bibliography, is discussed in Reference (1). As in 
all experimental techniques, thought must be given 
to designing a method which will yield maximum 
information for a given amount of computation. If 
random events are simulated, estimates of the pre- 
cision of computed output statistics are needed. 



Alternately stated, how many runs are needed, or 
how long should a single run be, to achieve a given 
interval of precision at a stated level of confidence? 
To answer these questions, the simulator must have 
considerable knowledge of probability and statistics. 

COST 

The cost of programming, debugging, and run- 
ning a simulation are usually high. Furthermore, 
it invariably requires more time — planning, testing, 
executing — ^than originally anticipated. The avail- 
ability of simulation languages may reduce pro- 
gramming time and cost somewhat. 

Simulation does offer some outstanding advant- 
ages, however. It can be used effectively in training 
and demonstration. Sometimes, valuable insight 
can be gained by simulating a system which could 
also be treated analytically. This insight may lead 
to a better analytic formulation. like a good doc- 



tor, lawyer, or teacher, simulation may achieve its 
greatest success by making itself superfluous. 

SIMULATION LANGUAGE 
A host of specialized digital-computer simulation 
languages, intended to facilitate the construction 
of programs ranging from discrete queueing models 
to analog representations of differential equations, 
are currently in use, and are constantly being im- 
proved. The reader is referred to References (i) 
and (2). 
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OPTIMIZATION METHODS — A REVIEW 
AND SOME EXAMPLE APPLICATIONS 

Brice Carnahan 
University of Michigan 
Ann Arbor, Mich. 



INTRODUCTION 

During the 9-v;eek 1965 summer program of the 
Michigan Project on Computers and Mathematical 
Techniques in Engineering Design Education, 
several days were devoted to formal presentations 
on mathematical optimization methods. Principal 
emphasis was given to those optimum-seeking and 
mathematical-programming methods which are best 
suited for implementation on digital computers and 
which in addition do not require a mathematical 
background much beyond that of typical senior-level 
engineering students. 

During the hnal four weeks of the summer pro- 
gram, the twenty-nine perticipating engineering 
design teachers formulated, solved on the computer 
(IBM 7090) and documented one or more design 
problems suitable for inclusion in their engineer!'’" 
design courses. Most of them attempted to use one 
of the mathematical optimization methods to pro- 
duce a design that was “best” in some sense. 

These problems have been published in the five 
discipline-oriented volumes of the final project re- 
port^ which have been distributed for the first time 
at this conference. In addition to these volumes, a 
summary volume (Volume 1) will soon be available; 
it contains, among several others, a lengthy review 
paper on mathematical optimization techniques. The 
paper includes a bibliography of 65 pertinent refer- 
ences which should prove useful to those interested 
in introducing optimization methods into the design- 
course environment. 

The present paper is a considerably less detailed 
review of some of the better known mathematical 
optimization techniques and is deliberately brief and 
incomplete. No doubt there are exceptions to some 
of the generalizations which appear. I hope that 
these remarks will give the reader the flavor of the 
subject and lead him to study the source materials 
for details and rigorous mathematical developments. 
Some example applications of the methods will be 
described briefly in the latter part of the paper; 



iKatz, D. L., Carnahan, B., Douty, R. T., Hauser, N., Mc- 
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puters in Engineering Design Education, Volume I - Sum- 
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Engineering, Volume IV - Electrical Engineering, Volume 
V - Industrial Engineering, Volume VI - Mechanical Engin- 
eering. 



those interested in details should refer to the ap- 
propriate discipline-oriented report. 

ENGINEERING DESIGN AND THE 
OPTIMIZATION PROBLEM 

Virtually every engineering design evolves in 
an iterative or trial-and-error fashion. The design 
problem is often stated in an ambiguous and open- 
ended way. The designer’s first task is to decide 
what the problem is, and what the requirements or 
specifications for any proposed solution are. The 
designer then generates a design concept, usr.ally 
in the form of a rough configuration, for the object, 
system or process being designed. These phases of 
the design activity are most difficult and require the 
greatest ingenuity and engineering judgement. The 
designer next attempts to “model” the object, sys- 
tem or process. This model may assume a wide 
variety of forms, but those usually preferred are 
mathematical in character. The engineer may use 
many weapons from his arsenal of mathematical 
and scientific tools, from the most fundamental of 
nature’s laws to the most empirical of data correla- 
tions, from the most rigorous analytical mathemati- 
cal tools to the least formal cut-and-try methods, 
to assist in the analysis of the proposed design. 

The mathematical model normally contains several 
(say n) parameters which we will call the “design” 
variables, Xi, Xa, . . .x„. To simplify the notation, 
we will use the notion of a design “vector,” x= 
[Xi, Xa, . . . Xi)]. We may view each design variable 
as one dimension in the “space” of the design vari- 
ably and any particular set of values for the design 
variables (i. e,, any particular design vector) as a 
point in this space. It may happen that certain 
equality constraints 

gi (x) =0 

gr (X) =0 

and/ or inequality constraints 

gr+l (x) ^0 

g„(x) 

are imposed on the space so that only points in a 
portion of the total space may even be considered 
as acceptable candidates for solution of the design 
problem. Such a region is termed the “feasible” 
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space; points in this region are termed feasible solu- 
tions or feasible design alternatives. 

The equality constraints come from functional 
relationships among the n variables which must 
be strictly satisfied; these often follow from elemen- 
I tary conservation principles. If there is to be a 
design problem, i. e., if there is to be a choice in 
the values which at least some of the design vari- 
ables may assume, then r must be smaller than n. 
The inequality constraints usually come from some 
specified design limitations (maximum permissible 
stress, minimum allowable temperature, etc.) There 
is no upper limit on the value of m. 

Briefly stated, the optimization problem is to 
find that one design vector or point in the feasible 
solution space which corresponds to the “best” de- 
sign. In order to locate the best feasible solution 
of the usually infinite number of feasible solutions 
using a mathematical optimization technique, it is 
essential that a meaningful, computable, quantita- 
tive, single-valued function of the design variables be 
formulated to serve as a measuring stick for com- 
parison of feasible design alternatives. This func- 
tion is usually called the “merit,” “effectiveness,” 
“criterion,” “utility,” or “objective” function. 

One of the designer’s principal responsibilities is 
to construct a suitable objective function to be maxi- 
mized or minimized. The objective function might 
be quite simple and relatively easy to calculate (the 
total weight of a storage tank to hold a specified 
amount of some material), or quite complicated 
and difficult to calculate (the rate of return over a 
20-year period for a proposed oil refinery, taking in- 
to account likely variations in markets, raw materi- 
als, taxes, etc.). The objective function may be 
very illusive due to the presence of conflicting or 
dimensionally incompatible subobjectives; for ex- 
ample, one might be asked to design a plant of mini- 
mum cost which at the same time minimizes air 
and water pollution and is aesthetically appealing. 
In such cases, it usually will not be possible to find 
a quantitative objective function and hence to auto- 
mate the selection of better design alternatives 
using one of the mathematical optimization pro- 
cedures. In what follows, I assume that a suitable 
objective function has been formulated by the de- 
signer. 

The mathematical optimization methods can be 
classified in a variety of ways; one such classifica- 
tion is shown below in Table I. 

^ TABLE I 

Classification of Optimization Methods 

1. Optimization methods for unconstrained 
functions of one variable. 

2. Optimization methods for unconstrained 
functions of n variables (n > 1). 
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3. Optimization methods for constrained objec- 
tive functions. 

4. Optimization methods which exploit system 
or information flow structure. 

Optimization Methods for Unconstrained 
Functions of One Variable 

Some optimization methods for unconstrained 
functions of one variable are listed in Table II. One 
criterion for locating a local extreme point of an un- 
constrained, continuous and differential objective 
function f (x) is that its first derivative f ' (x) van- 
ish. The traditional tool for solution of such prob- 
lems is the differential calculus, but its applicability 
in design problems is severely limited by the com- 
plex and nonlinear nature of most objective func- 
tions. 

TABLE II 

Optimization Methods for Unconstrained Functions 
One Variable — f(x) 

I. The Calculus 
II. Numerical Mathematics 
A. Root-Finding Procedures for f'(x) 

1. Newiton’s Method , 

2. Fals(j Position Methods 

3. Half “Interval (Bisection) Method 

4. Methods for Polynomial Functions ^ 

a) Graeffe’s Method 

b) Bernoulli’s Method 

c) Iterative Factorization Methods 
ill. Searching Schemes 

A. Simultaneous Methods 

1. Exhaustive Search 

B. Sequential Methods 

1. Dichotomous Search 

2. Equal Interval Multi-Point Searches 

3. Golden Section Search 

4. Fibonacci Search 

5. Lattice Search 

Some tools from numerical mathematics may 
prove useful, particularly the numerical root-finding < 

procedures listed in Table II. Newton’s method J 

requires the computation of the second derivative ^ 

f” (x) as well as the first derivative f' (x). The 
false postion and half-interval methods require only i 

that f'(x) be continuous and computable. Should i 

the objective function be a polynomial, then one of 
the listed root-finding procedures for polynomials 
might also prove useful. 

In contrast with these numerical method.3, the ■ 

searching schemes listed under parts IIIA and IIIB ! 

of Table II require only that the objective function | 

f(x) be computable; for most of these procedures, j 

neither differentiability nor continuity of the objec- i 

tive function is essential. 
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In each case, the optimization problem is to locate 
the global extreme value of an objective function 
f(x) on some arbitrary starting interval 
such that final optimal x value found by the 
algorithm and the true optimal x value lie within 
a specified “interval of uncertainty," i. e., an interval 
of specified length known to contain the extreme 
value. To simplify comparison of the methods, we 
assume that an arbitrary starting interval [A, B] 
has been mapped onto the unit interval [0, 1] by a 
suitable linear transformation. Hence, the initial 
interval of uncertainty for each method will have 
length 1. Other nomenclature used in common 
throughout this section is : 

a = length of final interval of uncertainty 

N = total number of evaluations of f (x) required 
to reduce the length of the interval of un- 
certainty from 1.0 to <* . 

n = the number of iterations for the sequential 
searching strategies (see below). 

The one-dimensional searching schemes may be 
classified into two categories termed simultaneous 
and sequential. In the simultaneous schemes, all 
points ai; which the function is to be evaluated are 
selected a priori. One such approach termed “ex- 
haustive search" is illustrated in Figure 1. Here 
the function, which need not be diflferentiable or 
even continuous in some cases, is evaluated at equally 

spaced intervals, a /2, assumed to be small enough 
to catch all pertinent features of the function. The 

total number of function evaluations required is 




Clearly, this is not a very effective approach for 
small a, particularly when the function is difficult 
to evaluate. 




When the sequential methods listed in Table II 
are used, the arguments for which the objective 
function will be evaluated cannot be known a priori ; 
instead, the sequence of argument values depends 
uuon the already observed values of f(x). All 
these methods require that f(x) exhibit unimodal 
character, i.e., exhibit just one local extreme point 



on the interval of interest as shown in Figure 2. 
The function need be neither differentiable nor 
continuous. 




Figure 2: A Function Unimodal on the Interval 

[ 0 , 1 ]. 

These sequential methods all involve the computa- 
tion of the function at one or more arguments 
which, because of the unimodality assumption, per- 
mit reduction of the length of the interval of un- 
certainty. A second set of functional values is then 
established for the reduced interval and yet a smaller 
uncertainty interval found. The process is repeated 
until the interval length is reduced to the specified 

value a. All these schemes have the property that 
the number of function evaluations required to 
achieve a specified terminal interval of uncertainty 
is known a priori, although the particular arguments 
used cannot be. 

Clearly, the evaluation of the function at just one 
point on the interval yields insufficient information 
to reduce the interval of uncertainty. At least two 
function evaluations at different arguments in the 
interval must be made. The location of these argu- 
ments is arbitrary, but some choices are better than 
others. A few of these one-dimensional searching 
schemes are illustrated in Figures 3, 4 and 5. 

If on the first pass of the sequential interval- 
slicing process we evaluate the function at two argu- 
meiiits centered about the midpoint and separated by 
a small distance e, i. e., compute f(l/2-|-e/2) and 
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f(l/2-e/2 as shown in Figure 3, then the smaller 
sub-interval of length (l-c)/2 can be rejected on the 
basis of the functional values (the I’s shown near the 
curve) . In this case, the interval [ (1-e) /2,1] becomes 
the new interval of uncertainty; two new function 
evaluations are made at arguments centered about 
its midpoint (the 2's shown near the curve) and 
the interval subsequently reduced again based upon 
the^ new functional values. This process, which 
is termed the dichotomous search, reduces the in- 
terval of uncertainty approximately by half for 
each iteration of the procedure. N and a are re- 
late to the number of iterations n as follows: 

o = + (1 - 2E“ ) E 

N = 2n = 2 1“ = 2.891n (^) 

Other choices of argument locations for the evalu- 
ation of the f (x) on each iteration are illustrated in 
Figures 4 and 5. . In both cases, the selected points 
are equally spaced in the interval. Figure 4 shows 
the two-point equal interval method applied to some 
unimodal function. On the first pass, the function 
is evaluated at x=l/3 and x=2/3) the I’s near the 
curve). That two-thirds of the original interval 
which must contain the extreme function value is 
retained as the new interval of uncertainty for the 
second pass; for the function shown, this would be 
the interval [1/3, 1]. On the next pass, the function 
is again evaluated at 1/3 and 2/3 of the span of the 
current interval, in this case at x = 5/9 and x ='i/9 
(the 2’s near the curve). Again, the 2/3 of the 
interval known to contain the extreme function 
value is saved, in this case the interval [5/9, 1]. 




Figure 4: Two-Point Equal Interval Search. 

The procedure is repeated until the interval of 
uncertainty is reduced to the specified a . The 
parameters «, n and N are related by: 




N * - 4.95 In (a) 



A similar three-point equal interval search is 
shown in Figure 5 below. In this case, the function 
is first evaluated at x=l/4, x=l/2 and x=3/4 (the 
I's near the curve) . That half of the original inter- 




val centered about the argument of largest function 
value (x=3/4 in this case) is retained. In the 
second pass of the sequential scheme, the function 
is evaluated at x=5/8, x=3/4, and x=7/8. In 
this case, and in fact for all the odd-point equal 
interval methods, there will be one function value 
from the preceding iteration which will appear 
again and need not be recomputed [f(3/4) here]. 
The parameters a , n and N are related by : 




N = 1 - 2.89 In (a) 

It can be shown that the three-point equal inter- 
val search illustrated in Figure 5 is the most etficient 
(i. e., requires the fewest function evaluations N 
for a given a ) of all possible equal interval search 
procedures. 

Intuitively, one would guess that rejection of one- 
half of the uncertainty interval per iteration of a 
sequential interval-slicing procedure would be the 
best one could achieve on the average, given no 
advance information about f (x) except for its uni- 
modal character. And from the viewpoint of the 
“minmax” strategy, i. e., the strategy which mini- 
mizes the likelihood of the worst possible outcome, 
this is the case. However, to reduce the uncertainty 
interval by half, using either the dichotomous or 
three-point equal interval search, requires two func- 
tion evaluations. Are there other yet more efficient 
interval-slicing methods? Yes! Two such methods 
described below are called the golden section search 
and the Fibonacci search. 

The golden section search is illustrated in Figure 
6. In this strategy, two selected arguments are 
positioned symmetrically about the center of the 



interval of uncertainty in such a way that t, the 
ratio of the lengths of successive intervals of un- 
certainty, 





Pi = a 



“ ^ 1-1 + h-2' i > 2 . 



where 






remains constant and is approximately equal to 
0.618, the so-called “golden-section” number. Hence, 
each iteration of this sequential technique reduces 
the interval by the factor 0.618, which would 
not appear to be as effective as the inteiwal 
halving achieved by the methods described above. 
However, clc ler examination shows that, after the 
first pass, one of the two arguments chosen for the 
j’th iteration will coincide with one of the arguments 
from the (j-l)’th iteration (e.g., x=0.618 in 
Figure 6) . Thus, only one new function evaluation 
is required for all iterations after the initial one. 

The parameters « , n and N are related by : 

« = (0.618)“ 

N = n-f-l = l-2.081n(a). 





0.764 

Figure 6: Golden Section Search. 



It has been shown that the very best strategy for 
optimizing a unimodal function is the Fibonacci 
search, which is based on the Fibonacci number 
sequence Fi where 



The first few numbers of the Fibonacci sequence 
generated by this recursion formula are : 



1 
0 
1 

2 

3 

4 



F. 

1 

1 

2 

3 

5 



i F. 

5 8 

6 13 

7 21 

8 34 

9 55 



The arguments at which the objective function f (x) 
is to be evaluated at each iteration of this sequential 
strategy are determined by ratios of appropriate 
Fibonacci numbers. Given a starting interval of 
length F„, the Fibonacci search requires, at most, n 
function evaluations to reduce the interval of un- 
certainty to length 1. 

To simplify the explanation of the method some- 
what but without any real loss of generality, con- 
sider an objective function f(x) on the starting 
interval [0,13] shown in Figure 7. Suppose it is 
desired to reduce the length of the interval of un- 




certainty to 1. Then one examines the Fibonacci 
sequence, finds the F„ which corresponds to (or is 
next greater than) the length of the initial interval, 
13. Since F# = 13, no more than n =6 function 
evaluations will be required. The first two argu- 
ments are positioned by the ratios F4/F« and Fg/Fe, 
namely, at 6/13 and 8/13 of the span of the inter- 
val, or at X = 5 and x = 8 (the I’s shown near the 
curve). Since f(6) >f(8), the new interval of 
uncertainty for the second pass of the iterative 
scheme is [0,8] . The arguments in the new interval 
are positioned by the ratios F3/F5 and F4/F8, i.e., 
at 3/8 and 6/8 of the span of the interval, or at x=3 
and x = 6 (the 2’s shown near the curve) . In gen- 
eral, the ratios involved on the j’th iteration are 



me 
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\ given by and F„.j/F„.j,i. One of the 

I arguments for the j’th pass will always coincide with 

\ one of the arguments for the (j-l)’th pass. Hence, 

I . just one new function evaluation need be made for 

I each interation after the fi '^it. On the last (n=5)’th 

F pass, one of the arguments is located within a small 
I distance e of the other to halve the interval of un- 
certainty after the (n-2=4)’th pass. 

If the problem is reformulated in terms of the 
I starting interval of length 1 used to illustrate the 

j other searching strategies, the parameters a , n and 
N are related by: 

a = i/F„ 

N = n. 

For i > 5, the ratio Fi_i/Fi = 0.618; thus, the 
method is closely related to the golden section search, 
and usually no more than one additional function 
. evaluation can be saved by using the Fibonnaci 
f search rather than the golden section search. 

I A closely related search procedure called the 

“lattice search" can be used to find the maximum of 
a discrete valued function which is unimodal. 

' A comparison of the number of function evalua- 

' tions required to reduce an initial unit interval to 

I a specified value of a is shown in Table III. The 

sequential strategies are clearly superior to the 
exhaustive search. 



TABLE III 

Comparison of Number of Function Evaluations 



a 


Exhaus- 

tive 


3-Point 
Even Int. 


Dicho- 

tomous* 


Golden 

Section 


Fibo- 

nacci 


0.1 


19 


8 


7 


6 


6 


0.01 


199 


15 


14 


11 


11 


0.001 


1999 


21 


20 


16 


16 


0.0001 


19999 


28 


27 


21 


20 



“Depends on e. 

One might reasonably ask why so much emphasis 
has been placed on the one^dimensional uncon- 
strained strategies, when most engineering design 
problems will have many design variables and prob- 
ably constraints on the design variable values as 
well. The principal reason, which will become 
clearer later, is that many of the multidimensional 
strategies involve sequences of uni-directional 
searches, and these one-dimensional methods can 
be adapted very easily to handle the suboptimiza- 
tion problem. In addition, although these one- 
dimensional searching strategies have been devel- 
oped without reference to constraints on the space 
of the variable, the selection of a starting interval 
places an arbitrary constraint on both the largest 
and smallest argument for which the objective 
function will be computed. Hence, the methods 
are directly applicable to constrained one-dimen- 
sional problems if the starting interval is chosen 
properly. 



GEOMETRY OF MULTIDIMENSIONAL 
SURFACES 

Some of the geometrical properties of a two- 
dimensional function f(xi,X 2 ) =. f(x) are illustrated 
in Figures 8, 9 and 10. A point x in the two-dimen- 
sional space, i.e., in the X 1 -X 2 plane of Figure 8, 
yields an observed value for the objective function 
f(x) which is displayed in the third dimension. 




(0) (ii) 

Figure 8: Response Surface and Level Contours 

for a Function of Two Variables. 



perpendicular to the space of the variables. The 
locus of objective function values in the third dimen- 
sion is sometimes called the “response surface." The 
intersections of planes parallel to the Xi-Xa plane 
(planes of constant function value) and the response 
surface, when projected onto the problem space of 
the Xi -X 2 plane, produce curves of constant function 
value called “level contours." Assuming that f(x) 
is everywhere continuous, the level contours will 
form connected curves ; should the function be every- 
where differentiable as well, then the level contours 
will be smooth curves (see Figure 8b). 

If f(x) is continuous and differentiable, it can 
be expanded about some point, say a, (see Figure 
9a) in a Taylor’s series. If only the linear terms 




Figure 9: Tangent Plane, Contour Tangent and 

Gradient for a Function of Two Variables. 
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in the series are retained, then the equation of the 
tangent plane of the response surface at a will be 
generated as 

y(x) *= f (a) + (i) + (x - a2)f (a) 

If the intersection of the contour plane f (a) and 
y(x), the tangent plane at a, is projected onto the 
Xi-X 2 plane, a line tangent to the level contour at a 
called the “contour tangent” is generated. Its 
equation is 

^*1 ■ + (^2 " = 0 . 

The projection onto the Xi -Xa plane of that line in 
the tangent plane along which the function y(x) 
increase m< st rapidly per unit distance, is called 
the ^direction of steepest ascent” or simply the 
“gradient.” The gradient at a is orthogonal to 
the contour tangent at a and its direction cosines 
di and ta are given by: 



T 

(x -a) A(x - a) , 

is definite. Here A is the matrix of second partial 
derivatives 






A = 



XgXi 



^2^2 



‘Vl 



n n J 
S 



If A is positive definite, the extreme point is a 
minimum point; if negative definite, the extreme 
point is a maximum point. 



> '^2 = 



S= /(f^^(3))^ + If (S))^ 

The n-dimensional generalizations of the tangent 
plane and contour tangent are “hyperplanes,” the 
n-dimensional analogue of a plane in three space. 
While impossible to picture graphically, their mathe- 
matical formulations follow in a simple way as : 

y(x) = f(a) + 2 (x - aj^)f (a) 

and n 

= 0 

The direction cosines of the gradient, which remains 
a line even in n-dimensional space, are given by 

= ^x^ (a)/S , 

where 



S = 







A point a is a “stationary point” when all the 
first partial derivatives vanish, i.e , when fx, (a) = 
0, i = 1,2, . . . , n. An “extreme point” is a station- 
ary point at which the quadratic form 



The concept of unimodality carries over from one 
dimension into n-dimensional space in a natural 
way; unimodality implies that the response surface 
has just one extreme point in the region under con- 
sideration. Three kinds of unimodality are illus- 
trated in Figure 10. A function is said to be uni- 
modal if, given an extreme point a, there is some 
path from every other point x in the space along 
which the objective function values strictly rise 
or fall (see Figure 10a) . Figure 10b shows the level 
contours for some two-dimensional function. Given 
an extreme point a, it is not possible to find a 
strictly rising or falling path from point b to point a ; 
therefore the function is not unimodal. A somewhat 
more stringent form of unimodality called “strong 
unimodality” requires that the straight line path 
from every point x to the extreme point a be strictly 
rising or falling. This is illustrated in Figure 10c. 
A yet more restrictive form of unimodality, ^linear 
unimodality,” is illustrated in Figure lOd. In this 
case, all straight line paths between any two points 
in the space (e. g., cd) must yield objective 
function values which are unimodal. 

While a strictly unimodal function would have 
just one “local” extreme point equivalent to the one 
and only “global” extreme point, it is often possible 
to confine observation to a small region of interest 
in which the only concern is that the local behavior 
of the function be unimodal. Most of the optimiza- 
tion procedures to be described in the next section 
require that the function exhibit unimodal character, 
at least locally ; some require more restrictive geo- 
metric properties. Virtually all are more efficient 
when the function is strongly or linearly unimodal. 
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*Z X2 




{c) Strong Unimodality (d) Linear Unimodatity 

Figure 10: Unimodality for a Function of Two 

Variables. 



OPTIMIZATION METHODS FOR UNCON- 
STRAINED FUNCTIONS OF n VARIABLES 
(n > 1) 

Several optimization methods suitable for finding 
extreme values for objective functions of more 
than one variable are shown in Table IV. 

The calculus is the classical tool for solving 
extreme problems in n dimensions. Unfortunately, 
its usefulness in the engineering design area is 
rather limited because of the difficulty of evaluating 
high order derivatives and subsequently of solving 
systems of usually quite nonlinear simultaneous 
equations. 

TABLE IV 



B. Sequential Methods 

1. Methods Requiring Evaluation of 
f(x) only 

a. Imbedded Random Search 

b. Lattice Search 

c. Univariate Search 

d. Rotating Coordinate Method 

e. Direct Search 

2. Methods Requiring Evaluation of 

a. Contour Tangent Method 

b. Steepest Descent Methods 

c. Optimal Steepest Descent Methods 

d. Partan Methods 

Since one of the criteria for finding a stationary 
point a of a continuous and differential function is 
that all the first partial derivatives vanish there, 
it may be possible to use one of the numerical 
root-finding procedures to solve the generated sys- 
tem of nonlinear equations fxi (a) = 0,1 = 1,2, . . ,n 
for the unknown ai,a 2 ,. . . ,an- The Newton-Raphson 
procedure is a common iterative solution method for 
systems of nonlinear simultaneous equations ; it has 
excellent convergence properties (when it con- 
verges) and is probably the best method when appli- 
cable. Unfortunately, the method requires that 
all second partial derivatives fxiXj (x) be computa- 
ble, preferably from analytical expressions. Numeri- 
cal approximations of the required derivatives may 
be adequate in some cases. 

Fortunately, all of the simultaneous searching 
methods and som.e of the sequential searching meth- 
ods listed in Table IV only require evaluation of 
f(X). Hence, these methods are particularly useful 
where derivatives are difficult to evaluate, the usual 
case in engineering design work. The multidimen- 
sional methods will be discussed in the order shown 
in Table IV. 



Optimization Methods for Unconstrained Functions 
of n Variables (n>l) - f(x) 

I. The Calculus 

A. Root-Finding Procedures for (x) 
1. Newton-Raphson Method 

II. Numerical Mathematics 

III. Searching Schemes 

A. Simultaneous Methods 

1. Exhaustive Search 

2. Random Search 



The n-dimensional analogue of the one-dimension- 
al exhaustive search already discussed requires the 
evaluation of the objective function at all nodes of 
a grid structure overlaid in the region of interest. 
For an n-dimensional hypercube (the n-dimensional 
equivalent of a rectangular area in two space or of 
a cube in three space) where a 1 is the desired inter- 
val of uncertainly^ for variable Xi and a,:^Xi^bi, 
the total number of function evaluations required is 



f2(bi - a^) “I 

N«TT [ -1 J 

i-1 



6S 



o 



Note that even for a three-dimensional problem 
with a rather modest starting hypercube (say 
bi -aj =1) and rather large interval of uncertainty 

(say «, =0.001), the toUl number of function 
evaluations required is about 8 x 10^; from a 
practical standpoint this approach is virtually use- 
less. 



For the general n-dimensional problem, assume 
that the feasible region is an n-dimensional hyper- 
cube with sides of length di in variable Xi, and that 
the final desired interval of uncertainty in variable 

X| is <*;. Then the ratio of the volume of the final 
“hypercube of uncertainty” to the starting hyper- 
cube !s given by 



a 




«1 



If we now choose a random point 3T, i.e., a point in 
the space for which the n coordinates have been 
selected randomly, then a is also the probability that 
3f is the hypercube of uncertainty containing the 
extreme value of the function. The probability that 
a random x is not in this hypercube is 1 - a ; and the 
probability that not one of p such randomly gener- 
ated points will fall into the optimal hypercube is 
(l-a)**. Then the probability s that at least one 
of the p random x values will be in the optimal 
hypercube is 

s = l- (l-a)p. 
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AREA ■ 0.01 



Figure 11: (a) Random Search, (b) Imbedded 
Random Search. 



Assuming that the at are much smaller than the 
corresponding dS 

_ ln(l - s) 
a 



Some values of p for selected values of a and s are 
shown in Table V. An illustration of this “random 
search” approach for a two-dimensional problem 
is illustrated in Figure 11a. 



TABLE V 

p (No. of Random Points) 



a 


00 

o 

II 

GO 


s=0.9 


8=0.95 


8=0.99 


0.1 


16 


22 





44 


0.05 


32 


45 


59 


90 


0.025 


64 


91 


119 


182 


0.010 


161 


230 


299 


459 


0.005 


322 


460 


598 


919 



As described above, the number of function eval- 
uations required for the random search increases 
exponentially with the number of dimensions, hence 
it is no real improvement over the exhaustive strat- 
egy described earlier. 

A great reduction in the number of random points 
required with only a marginal decrease in the over- 
all probability of success can be effected by imbed- 
ding the random search to give it the sequential 
character illustrated in Figure 11b. Here the original 
rectangle of uncertainty with sides of length 1 
and area 1 may be reduced to a smaller rectangle 
of uncertainty with aides of length 0.316 and area 
0.1, using just 44 function evaluations with a prob- 
ability of success 0.99 (see Table V). If this 
subregion centered about the best point (based on 
the function value there) in the first pass is subse- 
quently isolated and an additional 44 random points 
are generated in it, then a second rectangle of un- 
certainty with sides of length 0.1 and area 0.01 can 
be generated with a probability of success of 0.99. 
The overall probability for success will be degraded 
somewhat to (0.99)= = 0.98, but thO total number 
of function evaluations required is only 88. To 
achieve a comparable likelihood of success using just 
one iteration of the procedure would require 392 
function evaluations. 

The imbedding process can be carried out indefi- 
nitely with some decrease in the overall s at each 
iteration. There are no restrictions on the character 
of the objective function, i. e., differentiability, con- 
tinuity and unimodality are not required. However, 
since the only criterion for choosing a given random 
point as the center of a new block for further im- 
bedding is the value of the objective function there, 
the probability holds strictly only when the function 
value at every point in the hypercube of uncertainty 
about the extreme point is larger than at every 
point outside this hypercube. 
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I 

‘ SEQUENTIAL METHODS REQUIRING 

^ EVALUATION OF f(x) ONLY 

In the following section, the lattice, univariate, 
rotating coordinate and direct search methods, 
I which do not require computation of derivatives of 
the objective function, will be described briefly. 

Lattice Method'. In this method a grid of selected 
coarseness is overlaid in the region of interest (see 
Figure 12). A node of the grid is selected as a 
starting point. The function is evaluated there 
and at the 3"-l ‘‘adjacent" node points where n is 
the number of dimensions. The point with the 
greatest (or least) functional value is then selected 
as the next starting point and the procedure re- 
peated, except where the functional values already 
computed need not be recomputed. 

The process is continued until the greatest (or 
least) functional value occurs at the central point 
(starting point) . In this case, the grid size is halved 
in each dimension and the search begun with the 
reduced grid mesh. The grid size reduction is re- 
peated, when necessary, until the inter-point dis- 
tance in each dimension is smaller than some speci- 
fled tolerance value. 



k 



't 

$ 

I 




Univariate Search: In this search only one vari- 
f able is changed at a time. Typically the variables 
are changed in order, i.e., Xi,Xo, . . . ,x", so that the 
function is minimized (or maximized) in the co- 
ordinate directions. Starting with some arbitrary 
point 3Co, successive intermediate points Xi,X 2 , . . . , 
3L are found by carrying out linear minimizations 

I in the directions 

k" 
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0“ 




o’ 


0 




X 
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»PX ■ 
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» • • • Pn-x ■ 
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• 




0 


0 




0 




_x. 



such that Xui is the point of minimum functional 
value along the one in direction p, passing through 
point 3c t. Thus 

^ X 

where | a i| is the distance between and 
Although one need not compute the gradient 

g(Xw,) = gv,j 

to find a u a.\ is determined by the relationship 

i.e.. Pi is tangent to the contour at 5ci„. After com- 
pleting one round, a second round is completed, etc., 
until all I All are smaller than some corresponding 
tolerance values 

The method can be defeated if the variables 
Xi, . . . ,x„ interact strongly (Figure 13,) such as with 
a ridge-like structure for which the normalized gra- 
dient vector on the ridge has all direction cosines 
approximately equal in magnitude. If there is no 
interaction at all, i. e., if the principal axes of the 
surface are in the coordinate directions, then the 
method will find the extreme value in the least num- 
ber of rounds. 




(a) Little interaction. (b) Ridge with strong 

interaction. 

Figure 13: Univariate Search. 



Rotating coordinate method: In an attempt to 
overcome the ridge orientation problem of the 
univariate search procedure, Rosenbrock * developed 
a variant of the univariate method which involves 
rotation of the axes of the search. 



2 H.H.Roienbrock, “An Automatic Method for Finding 
the LeMt Value of a Function," The Computer Journal, 
Vol. 3, p. 176 (1960). 
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Heie the procedure begins with a standard one- 
round univariate search starting at point Xo (arbi- 
tiary), i,e,, n successive linear minimizations are 
carried out in the coordinate directions 



^i+1 ” + "i Pi» i ■ 0» •••» n-1 



such that 



p.^ln=0 



where 



0 

0 




1 



i*th position 



0 



Next, a new set of searching directions is generated, 
such that 



Pn 



n o 






and Pn+i 



» P2n-1 



satisfy the orthogonality conditions 

pT p * 0 » l>j*^» n+lj 2 n—l 

1 - J 

and 



-T - - 0. 

Pi 8i+l 

As before, 



The coordinate rotation technique used here 
causes the vectors Pnk to align them^ves along the 
direction of a sharp ridge, and hence the method 
tends to be quite good for tracking along a ridge, 
particularly along a relatively straight one. The 
method can be defeated however, given a bad start- 
ing point ^ 0 - 




Figure 14: Rotating Coordinate Method. 



Direct search: The direct search of Hooke and 
Jeeves* is a “trial-and-error” technique which has 
been used successfully by many investigators. There 
appears to be no one set of rules which constitute 
“direct search.” All the approaches consist of a 
local “pattern search” followed by a global move 
called a “pattern move.” 

The exploratory search consists of a restricted 
univariate search, such that at most one step (usual- 
ly small) of length is taken in the j’th direction. 

1. Start with an arbitrary base point bj^. 

2. Add (or subtract) 6 ]^ to x^ at and 
evaluate the function there. Let 



xi+1 * xi+ a± Pi 



In general, each round of the method consists of 
n one-dimensional searches in n orthogonal direc- 
tions. In the first round the n directions used are 
the coordinate directions. In the second round the 
first search is conducted in the direction given by 
a line passing through points Xo andxn; the following 
n -1 search^ are in directions mutually orthogonal 
and orthogonal to the first one as well. In succes- 
sive rounds new sets of orthogonal vectors R, . . . , 
Psn-i are generated so that for round k -f 1 



Pnk = 



^k ^(k-1) 

J [xnk - ^n(k-l) I I 



k 



1 , 



2 , 



• • • • 



The additional n -1 search direction vectors are then 
generated using the Gram Schmitt method. 




^1 



0 

0 



0 

and evaluate the function at + 63 ^ if 
necessary at and at bj^ keeping 
the "best" point. Let the new "tempo- 
rary" point be called til. 



*R. Hooke and T. A. Jeevei, “Direct Search Solution of 
Numerical and Statistical Problems,” J. A. C. M., 8, 2, 212- 
229 ( 1961 ). 



66 






{ bjL + 6^ if f(bj|^+6^) < f(b-j^) 

hi = I 4 f(bi-6i) < f(b^) 

{_ 

{ b-j^ if f(b^) < ipin{f (b^+6^), fCb^-T^)} 

3 . Next add (or subtract) 62 bo X2 at t^^^^ where 

0 ■ 

«S2 
0 
0 

• 

0 



^2 



o 



retaining the best point. Call it t^ 2 » 

{ *^11 ■*■ ^(^11 + *^ 2 ) < ^(^ 11 ^ 

^12 "“ { 4l ” *^1 ^^hl ~ ^ 2 ^ ^ ^^^ 11 ^ 

{ bli if fUll) < min {f(tll +^2^* f(tii- 62 )> 

4. Continue this process for each of the n variables Xj^, X 2 , ••• Xj^, 

each time finding a new temporary point tj^j where 

tij = if fdi.j.i-ij) < f(ti,j_i) 

1^1, j-l “ < min{(ti^j_l+6j)f(ti^j_l+ Ij)} 



5 . 




0 

0 



0 

0 



j * th element 



Let the second base point b£ be equal to t^^n* Note that the first 
exploratory search around the starting point b^^ requires at least 
n + 1 but no more than 2n+l function evaluations. 

Next, make the first pattern move by moving in the direction b 2 ”b 2 ^ 
along the line from bj^ to b 2 a distance equal to twice the distance 
from bj^ to b 2 » Call this point temporary point t2Q» i»e. 
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t2o - + 2(E2-^i) • 2^2 - 

6. Now conduct an exploratory pattern search about the point 72o as was 
done around the original base point bj^ (see 2 , 3 , and 4 above.) Call 
the terminal point of the search 7^ and let 

^3 ■ hn 

7 . Make the second pattern move in the direction b, - b, to locate a third 
temporary point ^30 

^30 “ ^2 2(b2-b2) = 

8. This process of conducting exploratory searches around the points 

^ko determining the k'th base point as b^+i"* followed by 

a pattern move to the point according to tj^^^ ^ =* 2bj^^^ - bj^ 

continues so long as the function value at is smaller than at 

9 . Should > f(bi^+i), then the base point hy, is dropped from 

consideration, and a new search is started by designating 

^1 • ^k+1 

and returning to step 2, 

10. Should a local exploratory search near the newly designated bj 

fail to find a temporary point tj^* different frc-m bj , i.e., should 
^2 ^1 then the sizes of the search increments 6j must be reduced 

and a new exploratory search made by returning to step 2. 

11 . When the 6^ have been reduced sufficiently to be "small enough," 
i.e., when 

where Gj is some small number, then the search is completed and the 
last designated b^ becomes the direct search approximation to the 
extreme point. 
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The method appears to be very efficient for many 
minimization problems, particularly those where 
ridges are encountered. The direction of successive 
pattern moves tends to become aligned with the 
ridge. Experimentally in such problems, the num- 
ber of function evaluafjions required tends to in- 
crease directly with n (the number of dimensions) 
rather than with some higher power of n as would 
be expected. 

A variant of the procc dure described above com- 
putes the temporary points as follows : 

^ k + 1,0 “ ■ "^1 • 

Thus all pattern moves are along lines with direc- 
tionj, bk+i-bx which pass through b,. The step sizes 
in the pattern moves are thus amplified more rapidly 
than in the standard method. 

It should be noted that since the criterion for 
tennination is failure of a univariate pattern search 
using sufficiently small 8j , the method can fail in 
those cases where n univariate search would fail, 
i. e., for surfaces wiith sharp ridges having certain 
orientations with respect to the coordinate axes. 

SEQUENTIAL METHODS REQUIRING THE 
EVALUATION OF FIRST PARTIAL 
DERIVATIVES 

In the following section, methods requiring the 
computation of first partial derivatives or the gradi- 
ent of the objective -function f (x) will be discussed. 
Where possible, the derivatives should be computed 
from analytical expressions. Since this is usually 
not feasible when dealing with complicated objective 
functions (e. g., the cost of a plant) , one must usually 
compute the derivatives numerically by perturbing 
csach variable in turn and computing the value of the 
objective function there. To compute all n first 
partial derivatives numerically (n l^ing the dimen- 
sion of the space) will usually require n -f- 1 func- 
tion evaluations. 

Method of Contour Tangents: If the function to be 
minimized (or maximized) is strongly unimodal. 



then the contour tangent line in two space (a plane 
or hyperplane in higher dimensional space) sep- 
arates the space into two half-spaces (Figure 16). 
The functional extreme point will lie in one of these 
two half-spaces. One “blocking” approach, called 
the contour tangent method, begins with the arbi- 
trary selection of a starting point Xo> computation of 
the first partial derivatives there, and subsequent 
determination of the contour tangent. Using the di- 
rection of the gradient as the criterion, one half of 
the total space, the half on the “wrong” side of the 
contour tangent, can be rejected as not possibly con- 
taining the extreme point. A new point Xi located 
in the “center” of the saved half-space is selected 
next and the process repeated. With each pass of 
the method part of the space is rejected. When the 
bounded or “blocked” space known to contain the 
extreme point is small enough, the search is com- 
plete. At each point selected, Xi, it is necessary to 
compute the equation of the contour tangent : 

(x - x^) g(x^) = f(^^) 

There are two major difficulties encountered in 
applying the method, selecting the locations of the 
Xi, and maintaining a running record of the extent 
of the space known to contain the extreme point. 
Wilde ^ suggests four different points, the “mid- 
pc»int,” “minimax point,” “median,” and “centroid” 
for Xi+i. From a practical standpoint only the mid- 
point can be calculated easily. Even so, the number 
of constraining hyperplanes may become very large, 
leading to an extremely complicated algorithm even 
for the two-dimensional case. 




Figure 16 : Method of Contour Tangents. 



'‘D. J. Wilde, “Optimization by the Method of Contour Tan- 
gents,” AICHE Journal 9, 2. p. 186 (1963). 
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Method of Steepest Descent : The successive points 
of the search Xi+i are in the direction of the negative 
gradient at 

Xi,"gi =-g(xi) 



The scaling problem can be illustrated by the fol- 
lowing simple function 

f (x,y) = x^ 4 - 25y2 

which, with the transformation u = 5y, becomes 



X • e • , ® i Pi * y X , • • • 

— ""Si 

where p£ = and is arbitrary 

but such that < f(x^) 

(see Figure X7). Justification for moving in the 
negative gradient direction is that the greatest im- 
provement in f (x) per unit length of travel, at least 
for infinitesimal ai, is in the negative gradient direc- 
tion. 



*2 




f (x,u) = x^ + u^. 

In the x,y space the gradient direction is not nearly 
the optimal direction (the direction toward the 
minimum) while in the x,u space it is the optimal 
direction. It is wise to scale the variables to make 
the contours as nearly circular, spherical, or hyper- 
spherical as possible. Since one can't very well 
know what the shape of the contours is ahead of 
time, such a scaling is generally not possible. One 
alternative is to choose a starting hypercube of 
uncertainty and then use the length of the i'th edge 
as a normalizing factor for the i'th variable. Scaling 
difficulties are almost always encountered when the 
variables are dimensionally incompatible (pressure 
and temperature, for example). All the gradient 
dependent methods which follow also have scaling 
problems. 

Optimal Steepest Descent Method: This perhaps 
misnamed method is similar to the method of steep- 
est decent, except that the oc , are determined 
automatically : 

*i+l = + “i Pi 

Si 

where Pf = - — 

Isi 

- T - 

Pi 8i+X " ® 



There are three principal difficulties with the 
method : 

X. Criteria for choosing the a , are not clear. 

2. Scaling can cause problems. 

3. It may take a very large number of steps to 
find the optimum value. 

One suggestion made by Marquardt ® which allows 
a more-or-less automatic step-size adjustment is to 
use the formula 

“i “ ®i-X^^X cos^ 0) 

where 6 is the angle between pi and pi_i and ai and 
as are constants (suggested values are ai = 0.5, as = 
X.6). Note that this causes ai>ai_i when the gradi- 
ent direction is not changing much (when Xi and Xi_i 
are far from the optimum) . 



^Marquardt, Chem. Eng. Progress, Vol. 66, No. 6 (June, 
1969 ). 



Thus the gradient direction at Xi is followed until 
the gradient of f (x) is orthogonal to it (see Fig- 
ure X8). 




Figure X8 : Optimal Steepest Descent Method. 
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A linear minimization in the direction pi must 
be carried out starting at point Xi. Any of the one- 
dimensional searching strategies may be employed 
to find Xi+i, the point of minimum functional value 
in the negative gradient direction. The method 
usually requires fewer functional evaluations than 
the method of steepest descent. This is particularly 
true for those problems which require that the 
partial derivatives be computed numerically. 

PART AN Methods: Because the steepest des<> 

cent methods often exhibit slow convergence to the 
extreme point, oscillating severely from one step to 
the next, particularly when the objective function 
surface possesses deep, narrow valleys, so-called 
acceleration methods have been developed which 
usually speed up convergence greatly. The methods 
are motivated primarily by some geometric proper- 
ties of quadratic functions. Since most continuous 
and differentiable functions assume a near-quadratic 
character in the immediate vicinity of an extreme 
point, the methods can often be used for non-quad- 
ratic functions as well. 

Consider the simple two-dimensional quadratic 
function : 



2 -2 
f (x^ , x^ ) ^X2 ^X ^2 ^22 '^2 

where = xj^ - Xj^ , X£ = X£ - X£ and 



X = 



4 



X 



2 



is the center (extreme) 



point of f(x]^, X£) 




It can be shown relatively easily that if f(x) is 
quadratic, then the contour tangents at Xi and X 2 
are parallel. Conversely, it can be shown that, 
if one has two parallel lines CTi and CT 2 and 
locates the point of minimum function value on 
each (say Xi and Xg), then the functional minimum 
is on the line passing through Xi and X 2 (see Figure 
19). 

In two dimensions the PARTAN method usually 
takes on one of the following forms : 

Method (a) (see Figure 19) : 

1. Arbitrarily choose a line CTi in the space of 
the independent variables. 

2. Choose any other line CT 2 parallel to CTi. 

3. Find the “best” point on CTi using one of the 
one-dimensional search techniques. The best 
point (call it Xi) is the point for which the 
line CTi is tangent to a contour of the objec- 
tive function. 

4. Repeat step 3 for line CTg and call the best 
point Xo. 

5. Search in the direction of the line passing 
through Xi and X 2 using a one-dimensional 
search strategy. If the function is quadratic, 
the best point on this line will be the functional 
extreme point. 

Method (b) (see Figure 20) : 

1. This method begins with two passes of the 
optimal gradient method as shown in the fig- 
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ure. Let Xo be the initial point and Xj and 
Xg be the terminal points for the first two 
optimal gradient searches. 

X2 =. Xq + tto Po 

■ 1^1 
Po’^ 82 ° 

X3 = X2 + tt2 P2 

82 

”2 = - in^ 

- 1 - 
P2 83= ° 



2. A final linear minimization is now performed 
in the direction Xo Xg along the line passing 
through Xo and Xg,. The minimum point x* will 
be the extreme point if f(x) is a quadratic 
function. This follows because Pa, the direc- 
tion of the contour tangent at Xg, is orthogonal 
to Po and hence parallel to the contour tangent 
at Xo- 

OPTIMIZATION METHODS FOR CONSTRAINED 
OBJECTIVE FUNCTIONS 

Figure 21 illustrates a typical optimization prob- 
lem in two dimensions with one nonlinear con- 
straint. Note that the unconstrained functional ex- 
treme point lies outside the region of feasible solu- 
tions. 




Figure 21 : A Constrained Optimization Problem. 



The methods available to the engineering designer 
for dealing with such problems vary in complexity 
and power, principally with the nature of the objec- 
tive function and constraints. Table VI lists a few 
of the methods available. 

TABLE VI 

Optimization Methods for Constrained Objective 

Functions 

I. Linear Programming 

II. Nonlinear Programming 

A. Lagrange Multipliers 

B. Iterative Linear Programming (Simpli- 
cial Methods 

1. Cutting Plane Method 

C. Adaptations of Methods for Unconstrained 
Functions 

1. Gradient Projection Method 

2. Multiple Gradient Summation Method 

Linear Programming 

A linear programming problem is a constrained 
optimization problem for which the objective func- 
tion and all constraints are linear in the problem 
variables. Very powerful algorithms have been 
developed for solution of such problems. Linear 
programming has in the recent past been the most 
used of all optimization methods, particularly in 
such business-oriented applications as resource al- 
location. 

Many industrially important problems are 
amenable to solution by linear programming meth- 
ods. Three of the more important types are known 
by the labels “The Transportation Problem,” “The 
Activities Analysis Problem,” and “The Diet Prob- 
lem.” A brief description of these three typical 
linear programming problems is given below. 

a. The Transportation .Problem: A company 

wishes to ship a number of units of some item from 
its warehouses to its customers (for example, a 
rubber company shipping tires to its wholesale dis- 
tributors). Each customer has ordered a number 
of units of the item; each of the company’s ware- 
houses can supply up to a certain amount of its 
current inventory of that item. What shipping 
schedule should the company use to minimize the 
total transportation cost, assuming that the total in- 
ventory available from the warehouses is equal 
to the total of the orders of all customers? 

b. The Activities Analysis Problem: A manu- 

facturer has fixed amounts of several resources 
available for his use (these might be raw materials, 
labor, equipment, money, etc.). The resources can 
be combined in several different ways to make one 
or more of several different possible products. It 
takes a certain number of units of a particular re- 
source to produce one unit of a particular product. 
He knows how much product he can make for every 
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unit of any particular product he manufactures. 
How should he allocate available resources to maxi- 
mize the total profit? 

c. The Diet Problem: Given the nutrient con- 
tents of several foods and a minimum daily dietary 
requirement for each nutrient, minimize the total 
food cost while maintaining at least a minimum 
nutrition level. 

Each of these problems can be written in the form 
of a standard linear-programming model of the 
form 

n 

maximize S Cj X. 
subject to the constraints 



Xz 




Figure 22: Geometric Interpretation of a Simple 

Linear Programming Problem. 



X^ ^ 0 f J “ 1^2^.. 



n { S } 

jSl { > } ‘'i’ "l ^ ^ 

The objective is to maximize the linear function of 
the n variables Xj, subject to m constraining rela- 
tionships also linear in the Xj ; the constraints can 
assume the form of equalities or inequalities (^,— ) . 
In addition, each Xj in the feasible solution must be 
non-negative. 

In two dimensions the linear programming prob- 
lem has a very simple geometric interpretation. 
Suppose the problem is to maximize the function 
f (xijXa) = 3xi -f- 4xo subject to the linear con- 
straints 

2Xi -f- 5x2 — 10 
4xi -f- 3x2 ^ 12 

Xx^ 0 
Xa^ 0 

This is illustrated in Figure 22. 

Notice that the inequalities require that the solu- 
tion be in or on the boundaries of the cross-hatched 
region. This region is termed the “feasible” region 
and any (x„ Xa) in the region is said to be a feasible 
solution. The linear inequalities enclose a feasible 
region which has the form of a convex polygon where 
“convex” implies that all points on the straight line 
joining any two points in the set of feasible points 
are also in the set. This simply means that a linear 
combination of two feasible solutions is also a feasi- 
ble solution, provided the combination lies between 
the two solution points chosen. 

For the given problem, the objective function is 
the family of straight lines with slope -3/4 ; its value 
decreases as the lines approach the origin. Note 



that the maximum feasible value of the objective 
function occurs at a vertex of the polygon (15/7,- 
8/7) . In fact, it will always happen (even in higher 
dimensional space) that the objective function for 
a linear programming problem will assume its ex- 
treme value at a vertex of the space of feasible solu- 
tions. All linear programming algorithms employ 
the strategy of examining the objective function only 
at vertices of the region of feasible solutions. The 
algorithm moves from vertex to vertex in such a 
way that the objective function value at the suc- 
cessor vertex is at least as great as at the predeces- 
sor vertex. In a finite number of steps the optimum 
can be found. 

The algorithms required to solve a linear pro- 
gramming problem are too complicated to describe 
briefly. Detailed descriptions can be found in many 
texts, among the more readable being Gass,® Had- 
ley,'^ and Dantzig.® 

Nonlinear Programming 

In the nonlinear programming problem, the ob- 
jective function or one or more of the constraints 
assumes a nonlinear form; the constraint functions 
are normally defined over continuous ranges of the 
independent variables. When the objective function 
to be maximized is subject to equality constraints 
and both the objective function and the constraint 
function can be differentiated, it may be possible 
to use one of the techniques of variational mathe- 
matics such Lagrange multipliers. Unfortunately, 
in most engineering design problems, constraints are 
of the inequality form and the objective function 



®Saul I. Gass, Linear Programming Methods and Applica~ 
tions, 2nd ed., McGraw-Hill Book Company, 1964. 

7G. Hadley, Linear Programming, Addison-Wesley, Read- 
ing, Mass. 1962. 

SGeorge B. Dantzig, Linear Programming and Extensions, 
Princeton University Press, Princeton, New Jersey, 1963. 
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is seldom easily differentiable. Consequently, La- 
grange multipliers have only limited utility for the 
more difficult engineering design problems. 

So-called simplicial methods have met with a cer- 
tain success. These all transform the nonlinear 
I programming problem into an iterative sequence 

of linear programming problems. One of these ap- 
P^oaches, known as the cutting plane method, is 
described briefly below. Still other approaches to 
the nonlinear programming problem involve modi- 
fications of methods for unconstrained functions of 
many variables, such as those already discussed. 
Two which appear to have achieved a substantial 
acceptance are the gradient projection method, a 
modification of the method of steepest ascent, and 
the multiple gradient summation method which is a 
modification of the direct search of Hooke and Jeeves. 

Cutting Plane Method: In this method » the gen- 
eral nonlinear programming problem maximize f(x) 
subject to 

gi(x):^0, i = 1, . . . ,m 
X = [Xi,X 2 , . . . ,x„] 

is rewritten as a problem, with a linear objective 
function and nonlinear constraint by adding one 
new variable and one nev/ constraint, 

ffm+1 * f (x) 

Then the problem is : maximize x„+i subject to 
gi (x):^0 , i = 1, . . . ,m 

Km+i— 0 

In the cutting plane approach, each constraint 
function, g,(x), is replaced by a linear approxima- 
tion obtained by expanding it about some starting 
point, 3c®, in a Taylor series and dropping all terms 
of order two and above. The nonlinear problem is 
then solved by first solving : maximize Xn+i subject to 
gi(x®) + gi (xo) (x-x»):^0 
i = 1,2, . . .m-j-1 

The solution, x*, of this linear programming prob- 
lem (which may or may not fail to satisfy all of the 
, original constraints) can be used as the base for 

new Taylor series expansions of the constraints, 
creating a new linear programming problem to be 
solved. When the solutions x*^ and for two 
successive iterations differ by less than some ac- 
ceptable amount and x*k+i satisfies the original con- 
straints, a solution to the nonlinear problem has 
been found. 

METHODS WHICH EXPLOIT SYSTEM OR 
INFORMATION FLOW STRUCTURE 

^ The methods described earlier have an essential- 

ly sequential point-to-point character; all the design 
variables are selected simultaneously for a given 
point, i. e., if there are n design variables, then the 

I Kelley, Jr., “The Cutting Plane Method for Solving 

Convex Programs," J. SIAM, 8, pp. 703-712, 1960. 
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problem becomes one of searching for the optimum 
in n-dimensional space. It frequently happens that 
the structure of the system can be exploited to allow 
the optimization problem to be reformulated as a 
stage-wise decision process, i. e., formulated in such 
a way that computation of the optimum solution 
takes on a serial character involving at each stage 
or step in the process the variation of no more than 
a very few of the n design variables at a time. One 
powerful solution approach to such problems which 
may involve a very large number of design variables 
is called dynamic programming. 

The stages are usually numbered in reverse order 
to the direction of information or material flow (see 
Figure 23). 




Figure 23: Dynamic Programming. 



The input or state variable to the i’th stage is the 
output state variable for the (i-fl)’th stage x^j. 

(xi+i might be a vector of state variables Xi+i.) A 
decision or design variable di (or vector of decisions 
di) is specified for the i’th stage. For each stage * 
there are two functional relationships. The first, | 
denoted hi describes the output state variable for 
stage i, Xi in terms of the input state variable Xi+i 
and decision made di. The second, denoted pi, i 

describes the profit (negative if a cost) associated I 

with the operation of the stage. The object of the dy- f 
namic programming algorithm is to determii:e how 
to operate the system (i. e., what decisions di,d 2 , . . . , r 
dn should be made) to yield a maximum overall pro- 
fit. 

Operation of Stage * . h^ (x^ , , d ) 

1 i 1+1 i 

Stage Objective ^ 

Function (e.g. Profit) = giCx^+^jd^) 

Optimal Cumulative Object Function; 
max.(p^(x.^^ dp + 

The procedure, which is a bit too complicated 
to detail here, is based on Bellman’s principle of 
optimality, which states: ' 

“An optimal policy has the property that, 
whatever the initial state and initial decisions ! 
are, the remaining decisions must constitute an | 

optimal policy with regard to the state resulting 
from the first decision.” 
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As a consequence of this principle, it is possible 
to optimize a serial system using a tail-to-head 
strategy, first examining a one-stage process, then a 
two-stage process, etc. 

Provided the program can be appropriately formu- 
lated, an apparently n-dimensional problem can be 
solved as a sequence of m-dimensional problems 
v?here (m « n) ; in some cases the n-dimensional 
problem can be reduced to a series of n one-dimen- 
sional problems. The number of computations re- 
quired can be reduced by several orders of magnitude 
over that required by standard optimum seeking 
methods when the dynamic programming approach 
is applicable. Even so, for problems with a large 
number of state and decision variables, at each 
stage the computations problem can become quite 
formidable. 

During the past year or two, and inspired pri- 
marily by the success of the dynamic programming 
concept, there have been several studies involving the 
decomposition of macrosystems, i.e., large and com- 
plex systems containing many components or sub- 
systems. The object of the decomposition schemes is 
to allow effective suboptimization of individual sub- 
systems, using virtually any of the various mathema- 
tical optimization techniques, both simultaneous and 
sequential, to reduce a problem of hopelessly large 
dimensionality to one of computable proportions. A 
recent paper by Wilde is particularly well written. 

DESIGN PROBLEMS 

During the last four weeks of the nine-week sum- 
mer program of the Project on Computers in Engi- 
neering Design Education, each of the participating 
faculty members generated, solved, and documented 
one or more design problems which he felt suitable 
for use in his design course. Most of these problems 
are included in the disciplinary volumes of the final 
report of the project. Each problem is stated as it 
might be presented to a student in an engineering 
design course and detailed descriptions of the solu- 
tion methods, computer programs and computer re- 
sults are included with each problem. 

Most of the faculty participants attempted to 
use one or more of the mathematical optimization 
methods of the previous pages to produce a design 
that was “best” in some sense. Table VII lists the 
volume number and general field of interest for 
those problems involving some optimization stratei?y. 
The number of the problem used in the reports, its 
title, the type of optimization method, the nature 
of the objective function, and the design variables 
involved are listed under each. Those interested 
in details should refer to the appropriate volume of 
the project report. 

10 D. J. Wilde, “Strategies for Optimizing Macrosystems,” 
Chem. Eng. Progress, 61, No. 3, pp. 86-93, March 1966. 

11 Comjyuters in Engineering Design Education, 



TABLE VII 

Engineering Design Problems Using Optimization 

Methods 

(Note) : The format adopted is: 

a. Optimization method. 

b. Objective function. 

c. Design variables. 

Volume II — Chemical Engineering 

I. Economic Optimum Conditions for a Staged 

Separation Process 

a. Univariate Search 

b. Rate of return = 

(Income - Product Cost) - Capital Investment 
Capital Investment 

c. XD = Product purity, ME = Multiple of 
minimum refiux, FR= Fractional recovery 

4. Design of a Mixed-Suspension, Mixed-Pro- 
duct Crystallizer 

a. Exhaustive Search 

Income-Total Cost 

b. Rate of return = Capital In T e l B iSr 

c. /=.„ Average crystal size 

Volume III — Civil Engineering 

6. Optimal Design of Two-Way Reinforced Con- 

crete Slabs 

a. Exhaustive Search 

b. Either weight of slab or cost of slab 

c. t = Slab thickness 

7. Decision Making in Building Planning 

a. Linear Programming, Simplex Method 

b. Cost of apartment complex 

c. Number of apartments built for each 
price range and size 

8. Optimization of a Two-Span Cover-Plated 

Steel Beam 

a. Nonlinear Programming, Cutting Plane 
Method 

b. Weight of the beam 

c. bf = Flange width, tf = Flange thick- 
ness, dw = Web depth, tw = Web thick- 
ness, /, = Left span length, t> = Right 
span length, bp = Cover plate width, 
tp = Cover plate thickness, etc. 

Volume I — Summary 

Volume II — Chewiical Engineering 

Volume III — Civil Engineering 

Volume IV — Electrical Engineering 

Volume V — Industrial Engineering 

Volume VI — Mechanical Engineering 

Final Report, Project on Use of Computers and Mathematical 

Optimization Techniques in Engineering Desig^n, D. L. 

Katz, B. Carnahan, R. Douty, E. McMahon, N. Hauser, 

J. Zimmerman, and W. Seider, March 1966, University of 
Michigan (800 pages). 
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9. Design v,f Rectangular Tied Concrete CoV- 
umns 

a. Exhaustive Search 

b. Cost of columns 

c. Size and number of longitudinal bars, 

, size and ^spreading of lateral ties 

10. Optimal' Design of Structural Steel Framing 

for Tier-Type Buildings 

a. Dynamic Programming 

b. Cost of a tier-type building 

c. Joist spacing, joist span, number of beam 
panels, number of girder panels, col- 
umn specifications 

11. Minimum Weight Design of a Linear Elastic 

Column 

a. Univariate Search 

b. Volume of material 

c. Column radius, length, and wall thick- 
ness 

12. Optimal Design of a Timber Warehouse Floor 

a. Dynamic Programming 

b. Cost of a warehouse floor 

c. Joist spacing, floor material, joist mate- 
rial, beam material 

13. Optimization in Site Selection 

a. Constrained Gradient Search 

b. Cost of a static missile test stand in- 
cluding utility and road construction 

c. Distance of missile stand from water 
sources, drainage pond, barricades, and 
roads 

14. Design Force for Bracing in a Deep Vertical 

Cut 

a. Direct Search 

b. Bracing design force 

c. Horizontal position of logarithmic spiral 
cut 

15. Optimal Design of a Wide-Flange Beam 

Using Standard Rolled Shapes 

a. Exhaustive Search 

b. Weight of the beam 

c. Beam dimensions 

Volume IV — Electrical Engineering 

17. Design of a Ferromagnetic Circuit 

a. Dichotomous Search 

b. Mass of the circuit 

c. B = Width of leg 3, K = Width ratio- 
leg 2/leg 3, etc. 

20. Class C Amplifier Design 

a. Gradient Search 

b. A trade-off function including input and 
output signal powers, and plate efficiency 

c. Plate supply voltage, grid supply volt- 
age, minimum plate voltage 



Volume V — Industrial Engineering 

22. A Dynamic Programming Problem for the 
Assignment of Tolerances 

a. Dynamic Programming 

b. Cost per combination within overall 
tolerance specified 

c. Component (stage) variaiices 

25. Optimal Assignment of Machines to an Oper- 

ator 

a. Exhaustive Search 

b. Cost per piece manufactured 

c. Machine cycle time, operator idle time, 
machine idle time 

26. A Dynamic Programming Approach to an 

Inventory Problem 

a. Dynamic Programming 

b. Coot of maintaining inventory 

c. Schedule for purchasing 

Volume VI — Mechanical Engineering 

29. Rocket Nozzle Design 

a. Dichotomous Search 

b. Weight of nozzle wall 

c. Wall thickness, materials of construc- 
tion 

30. Gear-Ratio Determination by Computer-As- 

sisted Search 

a. Cartesian Walk 

b. Precision of gear-ratio determination 

c. Number of teeth 

31. Optimization of a Pneumatic Linear Actuator 

a. Direct Search 

b. Operating cost (compressed air usage) 
usage) 

c. 11 variables including working area, 
stroke length, spring rate, working pres- 
sure, opening area, etc. 

32. Ejection Actuator 

a. Monte Carlo Methods 

b. Reliability of actuator system 

c. Resistance force 

33. Determination of Optimum Pipeline Route 

a. Dynamic Programming 

b. Cost of a pipeline route 

c. Number of sites, pump, and pipe sizes 

34. Cam Design of Minimum Size Based on 

Bearing Stress Limitations 

a. Direct Search 

b. Contact stress 

c. Prime circle radius 

35. Optimum Reheat Pressure for a Steam Power 

Cycle 

a. Fibonacci Lattice Search 

b. Heat cycle efficiency 

c. Extraction pressure 
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